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Introduction 



The investigation of non perturbative aspects of quantum field theories or statistical mechanic 
models is a fundamental topic in recent theoretical developments. These studies would give a 
more complete view of the structure of such theories (e.g. phase diagram), and are essential 
in order to elucidate non perturbative aspects like quark confinement in QCD or chiral 
symmetry restoration at finite temperature. Non perturbative techniques are both analytical 
and numerical. The most used numerical methods are Monte Carlo simulations that are very 
popular due to their universal setting.^ On the other side concerning with non perturbative 
analytic approaches, the 1/N expansions [4] have a central role, that seems to increase in the 
last years. The general scheme of this expansion (see chapter 1) is to consider Hamiltonians 
with a large number N of fields (that goes to infinity) in such a way that the equations of 
motion are simplified and can be solved analytically. Proceeding in such a way one recovers 
solutions similar to those found using other schemes hke Mean Field, variational methods 
or self-consistent approximations which reduce the interacting system to a interaction free 
one, with a self consistently determined parameter. For instance by considering a vectorial 
0^ theory, one can write 

and can recover a free interaction with a mass term that must be fixed in a self consistent 
way, by computing {(p). However the 1/N expansion with respect to the previous methods 
has the advantage that can be improved in a systematic way by including fluctuations. 
What does one expect from such an expansion? Of course one would extract informations 
about the physical case (at finite N), so that one would detect for instance the presence of a 
phase transition studying the N — oo solution. This is not always possible. Indeed, due to 
the fact that the large N limit puts infinite degrees of freedom in a finite volume, one could 
observe several non analyticities in the partition function that are due to the taken limit 
instead of the thermodynamic limit V ^ oo. This is the problem of the pathologies of the 
large N limit and has been investigated in detail in [11] by comparing the large N expansion 
with the available exact solution in d — 1. In the case in which the large N expansion is 
not fictitious one would also try to get 1/N numerical corrections to the N — oo solution 
to compare with available data at finite N (for instance MC simulations). Another usually 
addressed advantage of the 1/N approaches is that the theory can be solved (in certain case) 



^Maybe the claim is too optimistic. There are severe penalties related to MC simulations; the most 
popular are critical slowing down [1] or the sign problem in fermions system (e.g. [2]). 
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in a generic dimension,^ while for instance conventional renormalization group equations 
usually can be solved only when the fixed point is very near the Gaussian point (like for 
instance in the e expansion, e = 4 — d). 

In this work we want to investigate several aspects of the 1/N expansion for models with 
a global 0{N) symmetry.^ As stressed in [49] approaching a 1/A'" expansion one has to take 
care of two points (that are also two constraints of the theory): 

i) The theory must be both Infra Red Finite and Renormalizable. 

a) The large N expansion is just a technique, with its own limitation. The extrapolation 
of results (also qualitative) to finite must be discussed carefully. 

Point ii) deals with the question (introduced above) of the (possible) pathologies introduced 
by the large limit [11]. On the other hand, Renormalizability of large A^ actions [point 
i)] is quite delicate as the action for the auxiliary fields is usually not local. A scheme to 
circumvent this difficulty is presented in [24] and [49]. In this work we mainly focus ourself 
on problems related to Infra Red singularities [point i)]. The fact that 1/A" expansion is 
not Infra Red finite is usually addressed as a failure of the large A^ technique. There are 
several examples in which this happens. In this work we have studied the following cases: 
(a) Heisenberg models with generic ferromagnetic near neighbour interaction (chap. 1, 3, 
4); (b) fermionic models (Yukawa and Gross- Neveu) at finite temperature (chap. 5); (c) a 
tricritical model with a six-th order interaction (chap. 6). First, we will show (recover) that 
all the models introduced above exhibit a Mean Field phase diagram for N = oo. The 
failure of 1/A^ expansion calls for point ii): do the N < oo models show a critical point or 
the N = oo critical point is an artefact [11]? Scaling arguments for model (b) [50], a Monte 
Carlo Simulation for model (a) [82], and exact estimates for model (a) [79], suggest that the 
phase transition is present also at finite A^ with a critical behavior that is Ising hke. These 
results suggest that the phase transition near an Infra Red singular point is present also in 
the physical sector A^ < oo so that (at least in the present cases) the large A^ limit does not 
produce artefacts. However a question remains open on the nature of the Critical Point (if 
present). In this work we will confirm the results of [50], [82] and [79] that for every A^ finite 
one observes an Ising phase transition. However in the large A" limit the critical zone in 
which Ising behavior is observed scales with a proper power of 1/A^ reducing to zero in the 
spherical limit (A^ = oo). This explains why only Mean Field is observed for A^ = oo. This 
claim will be supported by the study of the effective interaction of the critical mode^ that, 
for the models presented above, looks hke a weakly coupled u(p^ interaction (with u ^ 1/-^) 
which exhibits an interesting crossover limit (called Critical crossover Limit) that matches 



^However in the approaches we present, as soon as one include 1/N fluctuations, there are severals 
distinguo in the scheme proposed. Of course this is not unexpected (we refer to section 2.1.5 fore more 
details). 

^For a review of 1/N expansion for this kind of models we recommend [49]. 

'^This is due to the trlog(- • •) terms that appear in the action for the auxiliary fields (see for instance 
chapter 5 cq. 5.7), and it is due to the integration of the physical fields. 

^Infra Red divergences are related to the appearance of a null eigenvalue in the inverse of the propagator 
P~^, which eigenvector is usually named zero mode. 
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Ising and Mean Field criticality. The mechanism is the same of what observe in medium 
range models [36] [34], where now plays the same role of R {R being the range of spin 
interaction) . 

The conclusion of this work are that for models (a), (b) and (c) the problem of Infra Red 
divergences can be solved using a generalized expansion which explains the universal^ 
crossover between Mean Field and Ising (or Multicritical Points in chapter 4). We believe 
that the scheme proposed is quite general and can be applied to other models characterized 
by the presence of a zero mode7 However in principle it could happen that the weakly 
coupled interaction could be unstable (this is related to a negative mass term in the zero 
mode interaction); until now it is not clear to us if this is the sign that the N — oo critical 
point considered is unstable (we have developed the theory around the wrong point!) or 
the sign that now really the N = (yD critical point is an artefact of the large N limit [point 
ii)]. Maybe, it seems that the answer to the previous question is not general but model 
dependent. Interesting also the fact that the method proposed gives a general scheme to 
compute 1/N corrections of the critical parameters, like for instance the critical temperature 
Tc. For example, in the critical models this can be done by studying the mass counterterm 
(and the magnetic one if the theory is not symmetric) that for d < 4 can be obtained by 
a computation in which only a finite number of diagrams enter^ (see chap. 2). This is a 
very important point in order to check large N expansions with MC simulations. In sec. 
3.5 we present some numerical predictions of a model that has been investigated in [82] by 
a MC simulation for = 3. The fluctuations seem to match the N = oo with the = 3 
results but there are not enough precision for any claims. As future plan we would compute 
corrections to the critical temperature for the model simulated in [50] . In this case there is 
the advantage (instead of [82]) that there are simulations for several N {N — A, 12,24) so 
that (we believe) the check will give more precise answers. 

Plan of the work 

The aim of the first four chapters is to introduce the generalized large formalism in the 
most general case (i.e. without any symmetry in the effective action). In order to do that 
we introduce Hcisenberg models with general nearest neighbour ferromagnetic interaction 

n = -A^/3^iy(l + <ri -o-,). 

(ij) 

In chapter 1 we show as for a large class of interaction W a continuous phase transition 
is present at finite temperature with Mean Field criticality. By an accurate choice of the 
interaction W, one is able to detect also multicritical points. However as soon as one tries to 



^Indeed we will show that the crossover functions of all the models considered are the same apart two 
obvious normalization condition. Details on how to compare numerically models a) b) c) with field theory 
are given in app. B. 

^This claim is supported by the general considerations presented in sec. 6.4. 

^This is due to the fact that in a bare <^'* theory, if d < 4 only a finite number of diagram must be 
regularized. 
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include fluctuations infra red divergences appear near the (multi) critical point and the 
1/N expansion fails. This problem has been solved in chapter 3 and chapter 4 (respectively 
for the critical and for the multicritical case) by a careful study of the effective action of the 
critical mode^ that looks hke a weakly coupled (p^ (</?^") theory. This exhibits an interesting 
universal crossover between classical to non classical criticality that has been just studied 
in the past for instance for medium range models [27] [34]. The studies of weakly coupled 
interactions have been presented in chapter 2. The discussion presented there, instead 
of the medium range models or field theoretical results, consider also the case in which Z2 
symmetry is not present in the starting interaction but must be fixed by other renormalization 
conditions. The presentation presented in the first part of this work follows [56] and [57].^° 
The problem of IR divergences is a well known problem in expansion [49] [14] [15]. 
This has suggested us to search to apply the same scheme used to study finite temperature 
(multi) critical points in Heisenberg models to investigate other physical interesting models 
that share the problem oi 1/N infrared divergences. In particular in chapter 5 we have 
focused our attention on system with Nf fermions coupled to a bosonic field through a 
Yukawa interaction, 

Nf 

/=1 

in the large Nf limit. The main result of this chapter will be to prove the claim that 
the crossover function obtained in the study of the Heisenberg model for A/" = 1 can be 
used also for the fermionic model (they are universal). We have just pointed out that an 
important check for the scheme proposed would be the computation of the 1/A^ correction 
to the critical temperature, to compare with the available results obtained in [50] for several 
N. In chapter 6 we will consider a vectorial model cf) — {0i, • • - ^at} with a cf)^ interaction 

The phase diagram is rather well understood since a long time (e.g. [96]). For H ^ there 
are two line of second order phase transition (that usually are called wing lines) that meet 
together in a tri critical point for H = 0. For the N = 00 theory a Mean Field Behavior is 
observed along the wing lines while, due to symmetry arguments, one expect Ising behaviour 



^The appearance of a critical mode at the critical point is strictly related to infra red singularities. 

^°In [56] the critical points in d = 2 have been explicitly computed. The results have been then generalized 
in a straightforward way to the d > 2 case. However in the case in which the interaction is not symmetric in 
higher than two dimension the discussion is more involved (see sec. 2.1.5) although the conclusions remain 
unchanged. 

^^Similar consideration holds for the Gross Neveu model (see [24] sec. 31.9), that is a model with a four 
fermions interaction. 

^^However in this case, instead of the Heisenberg models, it was pointed out (using scaling arguments) time 
ago [50] the mechanism of the critical zone reduction for which the region size in which Ising fluctuations 
become relevant go to zero for Nf 00, explaining the reason for which only Mean Field is observed for 
Nf = 00. 
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for every [96]. The apparent paradox can be solved including 1/N corrections in the same 
way of what we have done in the Heisenberg model. 

As future application we plan to apply our multicritical scheme to the study of Yukawa 
models with nonzero chemical potential and finite temperature which is believed to have a 
phase diagram with a tricritical point. This is an interesting topic in order to investigate 
the phase diagram of QCD. 

All the systems we have described are characterized by a global 0{N) symmetry. In 
models with local symmetry the Mean Field behavior of the large N hmit is questionable 
and currently under investigation [43] (see also [58]). It could be interesting to use our 
approach in order to recover (or not) previous results [43]. 



Chapter 1 



Large- expansions 



In this chapter we introduce the general technique commonly used in the study of models 
with 0{N) global symmetry in the limit in which goes to infinity. The basic idea deals 
with the fact that the 0{N) quantities self- average for — > oo; in fact this can be seen as a 
consequence of central hmit theorem [3] . Let we consider for instance the standard 0^ theory, 
where (f) is a, N component vectorial field and 0^ = (0^)^ is the interacting term. In order to 
compute 0{N) symmetric observables one can guess that in the large N fimit 0^ self averages 
on a given expectation value a{x) with fiuctuations that scale as 1/^/N. Concerning the 
expectation values, one can then write the equation reported into the introduction 



where now also fluctuations are included. In this work we implement the large limit by 
introducing several auxiliary fields that are related to 0{N) invariant quantities of the model 
[49] [24]. This allows us to write an effective action in terms of the fields introduced. In 
literature there are also other methods that give the same results like for instance Hartree- 
Fock approximations [5]. Also in the point of view of stochastic quantization or in the Critical 
Dynamic [8] the Langevin equation becomes linear and self-consistent for large N replacing 
0^ with {(f)^). On the other hand a famous related model is the spherical model [6], that was 
shown to share the same critical behavior with N = oo spin systems [7]. However Hartree- 
Fock or self-consistent conditions give only the N = oo limit while the method of auxiliary 
fields gives the possibility to develop a systematic 1/A^ expansion. This is a fundamental 
necessity for our work which focuses on 1/N corrections. 

As a prototype for IR divergences appearing in the large N limit, in sec. 1.1 we will 
introduce a general class of spin models (1.2) on a lattice, which exhibits a phase transition 
at finite temperature and spin-spin correlation length. In sec. 1.1 we also describe the recent 
interest in studying Hamiltonians (1.2). In sec. 1.2, after having introduced the auxiliary 
fields we study the saddle point equations (1.7-1.8-1.9) and the gap-equation (1.10) that 
give the physics for A?" = oo. We study both models with critical and multicritical phase 
transitions. Multicritical points can be observed choosing very peculiar Hamiltonians (1.2) 
so that the set of equations (1.15) can be satisfied at the critical temperature f3c < oo. We 
believe that the study of multicritical points is interesting because several physical systems 
undergo such a transitions. For example in the study of chiral symmetry in QCD several 
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Large-N expansions 



simplified models have been considered. In chapter 5 we investigate a Yukawa model^ at 
finite temperature T that shows the same IR singularities of (1.2) near a critical point Tc > 0. 
We plan to study the model also with a finite chemical potential^ fi > which is believed 
to exhibits a tricritical point so that our present investigation of multicritical interactions 
(1.2) could be useful. In sec. 1.3 we parametrize the gap-equation near the (multi) critical 
point obtaining the scaling fields and the equation of state (that are Mean Field like). In 
sec. 1.4 in order to include 1/A^ fiuctuations we develop the theory around the saddle point 
solution giving the expression for the propagator and the vertices of the Hamiltonian for the 
auxiliary fields. We show as near the critical point the propagator is singular at zero external 
momentum. In sec. 1.5 we show that, because this singularity, the standard 1/N expansion 
breaks down. 

1.1 The model 

Motivated by condensed matter or field theory, a lot of work has been focused on the inves- 
tigation of the following Hamiltonian 

H^-Nf3j2<^i-<^J, (1-1) 

(ij) 

where CTj is an A/'-dimensional unit spin and the sum is extended over all lattice nearest 
neighbours. Indeed (1.1) is usuaUy used as the prototype of short-range interacting models 
with global 0{N) symmetry. In two dimension the model is disordered for all finite P [62] 
and it is described for /? — >■ cxd by the perturbative renormalization group [63], [64], [65]. The 
square-lattice model has been extensively studied numerically [66], [67], [68], [69], [70], [71], 
checking the perturbative predictions [72], [73], [74], [75] and the non-perturbative constants 
[76], [77], [78]. Recently more general models than (1.1) have received attentions 

i/ = -Ar^^iy(i + (Ti-(T,-). (1.2) 

<ij> 

Indeed using exact estimates it has been pointed out in [79] that if one consider a generic 
0{N) interaction (1.2) then for a large class of interactions W, a first-order phase transitions 
appear for /9 > in two dimension. ^'"^ In [82] a particular class of interactions depending on 
a parameter p is considered and studied by Monte Carlo simulations 

W{x) ~ xP. (1.3) 

^This model in the large N limit [24] behaves exactly as the Gross-Neveu model [9] [10]. 

^This is relevant for recent experimental progress in the physic of ultrarelativistic heavy-ion collisions. 

^It [80] similar considerations are given for a gauge theory in 2 -|- 1 dimensions. 

^It is interesting to observe that this phase transition is not present for (1.1) although that the two 
continuum limit of (1.1) and (1.2) formally coincide. This is related to the fact that, according with [62], at 
the critical point the spin-spin correlation length is finite. 
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The results show, according with [79], the existence of a first order phase transition for 
high enough p at finite temperature.^ Otherwise also a continuous phase transition for 
p ~ 5 is observed and numerical evidences suggest that it belongs to the Ising universality 
class. It is important to stress that according with [62] at the critical point (CP) the spin- 
spin correlation length remains finite while the critical parameter seems to be the energy 
E = CTi ■ (Tj, where i and j arc neighbour. The same qualitative results of [82] was obtained 
in [55] taking the spherical limit {N — > cxo), but instead of the Ising criticality a Mean Field 
behavior was observed. 

The apparent paradox has been solved in [56] where it has been shown that in the large N 
limit the width of the critical Ising zone (that is the region near and (3c where fiuctuations 
become important and Mean Field behavior is suppressed) reduces as a power of for 
N — > oo. This explains why only Mean Field is detected for N strictly infinite. In this 
work we present the results obtained in [56] considering also interaction W depending on 
Af tunable parameter pi. We show the existence of multicritical points (MCP) for every 
model (with jV odd) that can simultaneously satisfy J\f + 1 conditions (1.14) at the critical 
point Pic - ■ -PMci f3c- Then following [56] we show how 1/N fiuctuations don't destroy the 
N = oo critical points so that CP or MCP are present at finite N and cannot be considered 
an artefact of the large N hmit [11]. Otherwise the Mean Field behaviour is destroyed as 
soon as one takes N finite. We explain this in term of a generahzed 1/N expansion using 
a weakly coupled ip^'^^ theory that will be introduced in chapter 2 and applied to (1.2) in 
chapters 3 and 4. 



1.2 The large- limit 

The large- A?" limit for the model (1.2) has been discussed in detail in Ref. [55]. Proceeding in 
a standard way [49] we introduce three auxiliary fields A-^^, p-r^, and /i^ in order to linearize 
the dependence of the Hamiltonian on the spins cr and to eliminate the constraint^ cr^ — 1. 
The partition function becomes 

J Hidp^^^dX,^] Hidii^da,] e-^^, (1.4) 

X 

where 

= [A^^ + \x^,crx ■ CTx+fM - Kf^pxf, + 2W{px^,)] + ^ X] {Pxcrl - Px) ■ (1-5) 

X/i X 

We develop the auxiliary fields near their saddle-point values 



Kn ^ Oi + ^=Xx^^, Pxn^T + ^Pxi^, 11x^1 +-^Px- (1-6) 



^There are other numerical works with similar results for different Hamiltonians than (f .3), see e.g. [83]. 
However Hamiltonian (1.3) is peculiar with respect other models -like for instance RP^~^ models- because 

it does not have any symmetry. 

^This constraint is usually irrelevant in the large N limit (see [24] sec. 30.6). 
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In Ref. [55] a, r and 7 were explicitly given. They can be written as 

7 = a(4 + m^)/2, (1.7) 
a = 2iy'(T), (1.8) 

^ - -(-0)^^ + ^-^, (1.9) 
i?i(mg) 

^ W^'(r,p)' ^'•'''^ 
where the parameter mo is related to the spin-spin correlation length = i/iTio and 



with the integral extended over the first Brillouin zone. The corresponding free energy can 
be written as 

F = -pdW{T) + ilog/(mg) + h{ml), (1.12) 

where 

L{ml)= [ log{p' + ml). (1.13) 
Jp 

In Ref. [55] it was shown that generic models may show first-order transitions. This 

happens when, for given /?, there are several values of itLq that solve the gap equation (1.10). 
Here, we will be interested at the endpoint of the first-order transition line for which the 
following two relations holds 



0. ^<0. (1.14) 



drriQ ' 9(mo)^ 9(mo)^ 

In order to recover the two previous conditions (1-14) one can use a family of one parameter 
of interactions W{x) = W{p; x) in (1.2) and tune both p and f3 to the critical point where 
(1.14) are satisfied. For instance in [82] and [79] the Hamiltonian 1.3 was considered as 
prototype of interactions W that show phase transitions. In [81] similar consideration to 
what obtained in [55] was given for the mixed 0{N)-RP'^~^ model W{x) ~ x + px^. In this 
work we will consider the most generic one-parameter families of interactions that satisfy 
(1.14) for a critical point Pc, fn^^. More interesting we will consider also families of interaction 
that depend on Af odd tunable parameters W{x) = W{pi, ■ ■ -pf^^x) in oder to avoid the 
following multicritical conditions 

0, ^. < (1.15) 



for i = 1, • • • jV -I- 1. This permit us to claim that for the Hamiltonian (1.2) not only Ising 

behavior can be observed but also multicritical points (described by scalar (fp^^'^ theory in 
chap. 2). This in principle gives us the possibility to study large- physical system with 
multicritical point. In the next sub-section we show that solutions of (1.15) with M even are 
unstable; i.e. if (1-15) is realised for mo = moc with Af even, then another solution m always 
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(a) 



(b) 




m, 



0,1 



mn 




mn 



m, 



0,1 



Figure 1.1: The two possibilities for P{ml) in the presence of a saddle-point with J\f even. 



exists with F{m) < F{moc)- The previous conclusion remains unchanged for the point m 
that satisfy (1.14) or (1.15) with opposite sign in the dis-equalities. Then in sec. 5.3 we 
will recover the N = oo phase diagram for the model (1.2) obtaining the scaling fields as 
functions of pi and (3 and explaining the mean field behavior. 



1.2.1 Saddle-point solutions with J\f even 

In this section we show that solutions of eq. (1.14) with Af even are of no interest. Indeed, 



given mo satisfying eq. (1.14) with J\f even there exists mg ^ such that 



F{ml,) < F{ml), 



;i.i6) 



where F{m^) is the free energy (1.12). Thus, is not the relevant solution and simply 
represents a mctastable state. 

We start by observing that, if M is even, is a local maximum or minimum of the 
function /5(mQ). Since Pimi) vanishes for rriQ oo, diverges for rriQ — > 0, and is always 
positive under the assumption that W'{x) is positive for 1 < x < 2, the function I3{m1) 
must behave as in Fig. 1.1. If it is a local maximum, there exists ^ such that mg ^ < 



m, 



0' 



f3{fn1) and (3(771^) < f3{mQ^) for ttIq < < TtTq. Similar considerations can 
be made if mg is a local minimum (see Fig. 1.1). Then, consider the free energy (1.12), a 
straightforward calculation gives 



dm?, 



;i.i7) 



Consider first case (a) of Fig. 1.1. Using eq. (1.17) we can write 

F{ml) - F{ml,) = -2/5(m^) [W{T{ml)) - W{T{ml,))] 
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+2 r dmlW'{r{mlmml)^. (1.18) 



■■0,1 

Now W'{r) > 0, dr/dml < 0, f3{ml) < f3{ml) in the interval, so that 



/ 



t2 



dr f"^o (jr 

dmlW\T{mlmml)—, > (3{ml) / dmlW\T{ml))- 



dml 1^2 dm?. 

0,1 " ™0,1 '-' 

= /3(m^)[W^(r(m^)) - W^(r(m^,,))]. (1.19) 

If follows 

F{ml) - F{ml,) > 0, (1.20) 

as required. In the case fnl is a local minimum the analysis is identical. If the disequalities 
(1.14) (1.15) have opposite sign one can repeat the discussion reported it is straightforward 
to recover the previous conclusion strictly following the case (a) of Fig. 1.1. 



1.3 Equation of state and scaling fields 

In this section we want to describe the phase diagram of the model outlining the mean field 
behavior. This is done in sec. 1.3.1 developing the gap-equation [3(171^) near the critical point 
(1.14). In sec. 1.3.2 we generalize the discussion for a multi-critical point (1.15). 



1.3.1 J\f=l 

We wish now to parametrize the singular behavior for fi ^ fic and p — > Pc- Expanding the 
gap equation (1.10) near the critical point we obtain 

P-Pc^Y. ""-rnip - Pcriml - mlr. (1.21) 

nm 

Because of the definition of (3c we have ooo = 0. Moreover, eq. (1-14) implies that aoi = 
ao2 = 0, ao3 < 0. For p — pcwe see that ml has the leading behavior 

ml-ml^ (^^J^ ' , (1.22) 



while for /3 = /3c, we have 



^0 - ^oc ~ I ■ (1-23) 

Q03 / 



The nonanalytic behavior with exponent 1/3 is observed along any straight line approaching 
the critical point, except that satisfying j3 — (3c — aio(p — Pc) = 0. Therefore, the correct 
linear scaling field is 

Uh = (3 - Pc- aio{p-pc), (1.24) 
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and we have ml — ~ Uj^ whenever Uh 7^ 0. To include the case Uh = 0, we write the 
general scaling equation 

ml-ml = ul/'f{x), (1.25) 

where f{x) is a scaling function and x is a scaling variable to be determined. Then, we use 
again the gap equation (1.21). Keeping only the leading terms we obtain 

Uh = a2Q{p-Pcf + ciii{p- pc)u]l^ f{x) + ao3Uhf{xf, (1.26) 

so that 

aosfixf + a^iip - p,)ul'''^f{x) - a^oip - Pc)W -1 = (1.27) 

Since (p — PcYu^^ — [{p — Pc)u'^^^^]^u]/^ , the third term can be neglected. Thus, we may 
take (the prefactor has been introduced for later convenience) 

X = an{p - Pc)\uh\-^^^ (1.28) 

The function f{x) satisfies 

ao3fixf + xfix)-l=0. (1.29) 

Such an equation is exactly the mean-field equation for the magnetization. Indeed, if we 
consider the mean-field Hamiltonian 

t u 

n^-hM + -M^ + —M\ (1.30) 
2 24 

the stationarity condition gives 

u 

-h + tM + -M^ ^0, (1.31) 

which is solved by M = h^^'^f{t\h\~'^^^), where f{x) satisfies eq. (1.29) with^ aqs = u/6. 
It is thus clear that the scaling field (1.24) corresponds to h, while p corresponds to the 
temperature. Note that this identification is not unique, since only the line — is uniquely 
defined by the singular behavior. For instance, in the usual Ising case, we could define 
t' = t + ah without changing the scaling equation of state, since t|/i|^^/^ = t'|/i|^^/^ + ah^^^. 
Since the scaling limit is taken with ^ 0, t ^ at fixed we sec that t\h\^'^/'^ ^ 

In the Ising case, however, there is exact Z2 symmetry and thus the natural t 
variable is defined so that it is invariant under the symmetry. In our case we could define 
Ut — p — Pc + A{P — /3c) and fix A by requiring the leading correction on any line (except 
Uh — 0) to be of order u'^J^ instead of order u]/^, recovering in this way an approximate Z2 
symmetry. For our purposes this is irrelevant and thus we will use p — Pc as thermal scaling 
field. 

Eq. (1.25) gives the leading behavior. It is also possible to compute the subleading 
corrections. We obtain for the leading one 

ml - ml = uff{x) + ufg{x) + 0{uh), (1.32) 

"^Notice that the sign of aos impHes that the low temperature phase is obtained -as expected- for p > Pc- 
However the sign of 093 is fundamental in order to guarantees the stability of the effective theory for the 
zero mode in chapter 3. 
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with 

9[X) — 2~? r~0 ft \2\ ■ [i.OO) 

Finally, let us discuss the singular behavior of the energy. We have 

E=2W{t,p). (1.34) 

Such a function is regular in ml and p. Since p — Pc ~ a^l-U/ip^^ and ttIq — ~ '^h '^t the 
leading term is obtained by expanding the previous equation in powers of ml — m^^. Thus 



o2 _ r> 



ABi 



ufm + Oiuf). (1.35) 

mo=moc 



where W' = dW{x,p)/dx with x — t and the suffix c indicates that all quantities must be 
computed at the critical point. 

1.3.2 The general J\f odd case 

Following the same step of the previous section 1.3.1 we want now the expand the gap 
equation (1.10) near a multicritical point (1.15). First we present the J\f = 3 case because it 
requires to consider non-linear terms in the definitions of the scaling fields, the generalization 
is then straightforward. We have 

p-p^= ^ afx,^{p, -picr{P2-P2cy'{P3-P3cr'{mi - miy. (im) 

with aQoo = for < £ < 4 because of eq. (1-14). To simplify the notations we introduce a 
multi-index A= (ii,i2,'?3), rewrite the previous equation as 

= J] a5J)(p - p,)\ml - mU. (1.37) 

and define the linear scaling fields as 

Uh^P-Pc- J2 (^AiP-Pc)^, (1-38) 

A:[A]=1 

«?(P-Pc)^ (1.39) 

A:IA] = 1 

^2= <^AiP-Pc)^, (1-40) 

A:lA] = l 

^3= Y ^AiP-Pc)^, (1-41) 

A:IA] = 1 

where [A] = a indicates indices ii, ^2, is such that ii + i2 + h = The gap equation becomes 
Uh = ^4oo(m^-mgJ^ + (mg-mgJiii + (m2-m2j2^i2 + (m2-mgj3ii3 (1.42) 

+E E «?(p-Pc)^K-o^+E E a^A^ip-Pcnmi-miy. 

e>4 A:IA] = 1 e=0 A:[A]>2 
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From this equation one can read the correct scaUngs. We have mg — mg^ ~ u]/^^ and, as a 
consequence, Uj ~ Uh{ml — rnU'^ ~ u]^^^^ . Moreover, expressing pi —pic in terms of Ui, we 
obtain pi — pic ~ ~ fJ^J,^ ■ Thus, close to the multicritical point we expect 

^0 -"^oc = (2^1, 3^2, 3:^3), (1-43) 

where the scahng variables X j are defined by Xj = UjU^fl^ ^ . If we now substitute this Ansatz 
in eq. (1-42) we find however a difficulty. All terms in the equation should vanish close to 
the MCP as Uh or faster. On the other hand, the terms with [A] = 2 and I — Q vanish only 
as u^h^ . Thus, they are the dominant ones and make the Ansatz (1.43) inconsistent. There 
is a very simple way out of this problem. It is enough to include these terms in the definition 
of Uh and consider the nonlinear scaling field 

UH = ^-^- Yl ^a\p-Pc)''- (1-44) 

A:[A]=1,2 

With this definition the dangerous terms are no longer present in eq. (1-42) and in the scaling 
limit, M/i — > 0, — > 0, at fixed Xi, we obtain 

4oo/' + ^3/' + + xi/ - 1 = 0, (1.45) 

which is exactly the mean-field scaling equation of state for a generic 0® theory with Hamil- 
tonian 

(5) 

^ = -0 + y 0^ + f 0^ + f 0^ + (1.46) 

Note the absence of the 0^ term, which, if present, can always be eliminated by performing 
an appropriate shift of the fields. 

The generalization to generic values of J\f is straightforward. We define 

Uh = (3-/3,- ^AiP-Pcf (1-47) 

A:[A]<mo 

^ ('AiP-Pcf (1-48) 

A:[A]<m,i 

X, ^ u,uf^^'^-\ (1.49) 

i = 1, . . . ,Af, which imply {pi — Pic) ~ m^^^"^"*"^^ and (mg — mg^) ~ u]/^''^~^'^\ The integers 
must be determined by requiring that 

Y 4\p - p^nml - mlY (1.50) 

A:[A]>m( 

is of order u'^ with a > 1 for < £ < jV'. It follows 

Tr 7. > 1, (1-51) 
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so that 



me 



Af-e 



+ 1, 



:i.52) 



where [x\ indicates the largest integer that is smaller than or equal to x. In terms of these 
variables we have therefore in the scaling limit 



^0 - ^Oc = 



i/(Ar+2) 



/(a;i,...,XAr), 



;i.53) 



where the scaling function satisfies the mean-field equation of state of a 0^+^ theory: 



(1.54) 



This discussion also clarifies why only the case J\f odd is relevant. If jV is even we obtain a 
0" theory with n odd, which is unstable, providing another indication for the instability of 
the corresponding solution. 



1.4 The 1/N calculation: propagator and effective ver- 
tices 

In order to develop the standard 1/N expansion we insert the fields expansion (1.6) into the 
action for the auxiliary fields (1.5). Taking a short-hand notation we define a five-component 
field 

* = (A, Ai,A2,pi,p2), (1.55) 



then H can be written as 



(1.56) 



where the indices Ai run from 1 to 5 and we have neglected the constant part -F(mg) (1.12). 
Notice that in (1.56) the linear term disappears due to the gap equations (1.7-1.8-1.9). The 
propagator can be explicitly written as^ 



P-'(P) 



/ 1^0,0 












\ 


-2^1,0 


2^2,0 




-\I3W'^ 







-2^0,1 




2^0,2 
























V 












/ 



(1.57) 



^It is useful to write the Fourier transform of as A(p) = e ^Pi^/^ ^ ^^'^^x- This makes all vertices 
and propagators real. 
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where W should always be intended as a function of T{ml), 

. / 9x f cos' Oj. cos-' o„ 

AAp,ml) ^ / 2 . (1-58) 

and p"^ = 4(sin^p2:/2 + sin^pj^/2). 

For p — > 0, by using the algebraic algorithm described in App. A of Ref. [74], it is easy to 
express the integrals Aij{0,ml) in terms of the integrals Bn{ml) defined in eq. (1.11) with 
n — 1,2. Exphcitly we have 

Aoo{0,ml) = B2, 

f m?\ 1 
Ao(0,m2) = Aoi(0,m2) = h + _o j ^2 - -B^, 

An(0, ml) = -^(4 + ml)B^ + ^(8 + Sm^ + m^)B2, 
o o 

^20(0, ml) = Ao2(0, ml) = \ + B2-\{A + ml)B^. (1.59) 

o o 

Vertices are analogously computed. It is easy to check that the only nonvanishing contribu- 
tions for which some Ai is equal to 4 or 5 are 

Vtlip,, . . . , Pn) = VtliPl, . . . , Pn) = -PW^^\f). (1.60) 

If all indices satisfy < 3, then 

<U„(Pl,---,Pn)5(EP')= (1-61) 

(-1)"+1 



n / ^(^i+i - - Pi)^ — 2^^i(Pi' 



+ permutations > , 



where 

-Ri(p,q) = l, -R2(p,q) = -cos(ga:+Pa:/2), i?3(p, q) = - cos(gy + Py/2). (1.62) 

The permutations should made the quantity in braces symmetric under any exchange of 
(pj, Ai) [the total number of needed terms is (n — l)!/2]. As already discussed in Ref. [55], 
at the critical point the inverse propagator at zero momentum has a vanishing eigenvalue. 
Indeed, a straightforward computation gives 

detP-i(O) = Ko,detSi, (1.63) 

where ^cdet is given by 

^o,det ^ - J^l^ [^B.W - (1 - ml{8 + ml)B,)W"] , (1.64) 

and 

d(3 1 



Si = 



dm, 
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Eq. (1.63) shows that the determinant vanishes at the critical point — hence there is at least 
one vanishing eigenvalue — since there dp/dml = 0. The corresponding eigenvector can be 
written as 

Aoi(O.m^) 1 1 \ , , 

'Aoo{0,miy ' '2W"'2W"J ^ ' 

computed at the critical point. Indeed, 

J](P-^)ab(0)zb = J^^ii^^ 1. 1. 0, 0)^. (1.67) 

B=l ^ 

Note that there is always only one zero mode. Indeed, at the critical point, we have from 
eq. (1.65) 

W" = -4^W', (1.68) 



Bt - B2 



so that we can write at criticality 



Bt[B ! - {8 + ml)mlBl] 1 
32(^2 - B2) W 



We have verified numerically that the prefactor of 1 /W'^ (that is independent of W so that 
is the same for all the interactions that one consider) is always finite and negative, so that 
-K^o.dct is always nonvanishing. Thus, if wc consider (multi) critical interaction, for pj = the 
determinant dct P^-'^(O) vanishes as (mg — mg)-^+^ that is the scaling of eigenvalue associated 
with the zero mode, cf. Eq. (1.67). Thus, there can only be a single eigenvector with zero 
eigenvalue. This means that the effective action for the zero mode will be (see chapter 3) a 
scalar interaction. However choosing different will change drastically the nature of this. 



1.5 The failure of expansion 

Now we want to show as the appearance of a zero mode (see sec. 1.4) invalidate the standard 
large-A^ expansion. Because of the zero mode, it is natural to express the fields in terms 
of a new basis. For each uIq and p, given the inverse propagator P^^(p), there exists an 
orthogonal matrix U (p; ml,p) such that U^P~^U is diagonal. If va(p; ^o,p) = Uai{p', Tnl,p) 
is the eigenvector that correspond to the zero eigenvalue for p = at the critical point and 
QAa(p;"^0'P) = ^A,a+i(p; "T.o,p), a — l,...A are the other eigenvectors, we define new fields 
^a(p) by writing 

B a 

where $ = (0, v^a)- Eq- (1-70) defines the fields up to a sign. For definiteness we shall 
assume va{p) to be such that, at the critical point. 



VAiO)^ZA/{Y,4y^'- 
B 



(1.71) 



1.5 The failure ofl/N expansion 



23 



We do not specify the sign of (fa since it will not play any role in the following. 

The new field corresponds to the zero mode, while the four fields (pa are the noncritical 
(massive) modes. The effective Hamiltonian for the fields $ has an expansion analogous to 
that presented in eq. (1.56) for ^. The propagator Pab{p) of $ is 

Pab{p) =Y,Pcd{p)Uca{p)Udb{p): (1-72) 

CD 

while the effective vertices are related to the previous ones by 

^iri.,A„(Pl> ■■■:Pn)^ ^BiLBn(Pl' " " " ' PuWe^aAPi] mo,p) • • • C/B„A„(Pn; mo,p). 

Bl,...,Bn 

(1-73) 

Note that Pab{p) is diagonal by definition, i.e., Pab{p) = SabPaa{p)- Relation (1.63) 
implies 

Ai(0) ~ ~ (p, - p,e), {ml - mlf+^ (1.74) 

close to the critical point. Now we want to show that due to the singular behavior of P (1.74) 
the standard 1/A^ expansion fails near the critical point. Suppose we want to compute the 
expectation value of the n-th point function of the critical mode Xn 

Xn = J dX2 ■ ■ ■ dx„ < Lp{0)lp{x2) ■ ■ ■ ^p{Xn) >connected (1-75) 

In order to do that one can use the effective action for the zero mode Ticfr, in which the 
massive mode (fa (1-70) have been integrated out. The study of Ties is the main subject of 
the chapter 3. There we will show that at the critical point Pi = Pd the effective vertices 
yin) behave as (in the infra-red limit, for j = 1 • • • A/" + 3) 

y(-)(p,, ml;p, = pe^) ^ Iniml - m^J^"" + «f^Pi + ' " ' (1-76) 

i 

where k = Af + 3 if one considers the zero mode defined in (1.70) or /c = A/" + 4 if one 
translate the zero mode (1.70) by a constant factor so that lj\f+2 = 0.^ If we consider then a 
general £-loops graph Dn/ entering into the computation of Xn [eq. (1-75)] with Nj j— legs 
vertices (j = 3, • ■ • , A/'+ 1) and A^int internal line, neglecting ultraviolet divergences, by using 
eq. (1.74) and eq. (1.76) the following scahng relation holds^° 

Dn,e ^ (mg-m^J-2"-2^-+'^+S/-"^'('=-^^) 

= (m^-m2j-2-+'=^-i-3M.t+(d+i)^^ (1.77) 

d being the dimension of the Euchdean space. In (1.77) we have used the topological relation 

Mnt = ^ + (1.78) 



^However in the considerations reported in this section we can take k undefined. 

^°In eq. (1.77)we neglect contributions coming from TV + 2-legs vertex, however the general claim remains 
unchanged. 
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and we have neglected possible logarithmic corrections. The expansion is a loop expan- 
sion around the saddle point solution. Consider for instance the d = 2 case, from (1.77-1.78) 

D^,e - {ml - m^J-2"+^^-i-3(E. iv.-D. (1.79) 

Looking at eq. (1.79) it is clear how if one considers higher order in the 1/N expansion 
new more severe algebraic infrared divergences appear while in principle one expect only 
logarithmic contributions. This is the clear sign of the breakdown of the standard 1/A^ 
expansion. In order to solve the problem one has to consider a non-perturbative approach 
based on the study of weakly coupled scalar theories which are studied in chapter 2 and 
will be applied to the effective zero mode interaction TleS in chapter 3. One can guess that 
divergences appearing in (1.79) cancel when explicitly Feynman diagrams are evaluated. One 
can convince that this is not the case considering for instance X2 at one loop. There are three 
different diagrams entering into the computation, one with a four-legs vertex (~ m~"^logm^) 
and two built with a couple of three-legs vertices (~ m~^, ~ m~^) so that there is no 
chance for any cancellation (we have defined m = ml — ml^). On the other hand this 1/N 
contributions cannot match with 1/N corrections of the N — 00 result X2 ~ ^/{ml — rnl^y 



X2 



(ml - ml. + f ) 



2+b/N 



, ; 9 NO \ - ^ T^TT log("^0 ~ "^Or)^ ~ ^7 ^ ^ I • (1.80) 

(m^-m2j2^ 2iV ° N{ml-ml)) ^ ^ 

In the previous equation (1.80) we generate (with respect to the leading order) only Am~^ 
and log Am divergences, while in principle, using the scaling arguments reported above, one 
would match m~^ and divergences. 

The problem of Infra Red Divergences is well known in large N expansion [49] and it is 
usually addressed as a failure oil/N expansion. They are not peculiar of the Heisenberg 
models introduced above (1.2) but are also present in other models like Gross-Neveu and 
Yukawa theory at finite temperature [12] [49] [96]. It is important to stress that for the 
previous models at zero temperature 1/A^ expansion works well [49] [24] [13]. In this case, 
concerning universal coefficients for instance, one finds the N — oo value and then standard 
perturbation theory is able to get 1/A^ corrections (e.g. [13]). This works exactly in the 
same way of 1/iV expansion of the vector model: one finds for instance^^ for = oo 
u = l/{d — 2) (e.g. [17] [18]) and fluctuations around the saddle point solution can be 
resummed giving a 1/N correction [19] [20] (e = 4 — d) 

1 2(3 -e) ^ 1 



4sin(7re/2)r(2-e) 
7rr(l - e/2)r(2 - e/2) 



However when the temperature is flnite the Infra Red Divergences become more severe so 
that the perturbative series cannot be useful [14] [15]. This can be understood also looking 



^^With 1/ commonly one refers to the critical exponent related to the critical behaviour of the order 
parameter M near the critical point M ~ {—t)'^, t being the reduced temperature t = {T — Tc)/Tc- 
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at the standard perturbation theory in which nonanalytic terms appear if the temperature 
is not zero [14]. For a discussion of the finite temperature see also sec. 10 of [16]. For the 
case we have investigated we suggest that one cannot correct the N — oo critical behavior 
because critical exponents are discontinuous for A?^ = oo. However one observes a crossover 
between the two different critical region that will be illustrated in the next chapter. 



Chapter 2 

Crossover in critical phenomena 



In the past decades a huge amount of work (experimental, theoretical and numerical) has 
been devoted to the study of the nature of critical points in physical system characterised 
by a large number of elementary components (since now "spins"). In particular -from the 
theoretical point of view- the possibility to apply field theoretic methods (e.g. [24], [25]) to 
implement basic physical assumptions like scaling invariance of the universal quantities at 
the critical point (e.g. [23]), has allowed to obtain accurate predictions for these observables 
in the scaling region r ^ 1 [r = (T — Tc)/Tc, T being the temperature].^ Outside the scaling 
region^ (physically dominated by correlations between the degrees of freedom of the system) 
one expects to recover the physical picture given by mean-field approximations which is 
supposed to be dominated by the short range interaction of the spins. Several works [27] [28] 
[29] [30] [31] [32] [33] [34] [35] [36] [37] have investigated the nature of the crossover between 
the mean field region and the non trivial critical region (for instance Ising in simple fluid). 
Using the language of Wilson renormalization group theory, the main effort of the theoretical 
works is to resum the infinity set of irrelevant operators that becomes important as soon as 
one leaves the scaling region (for instance 0^, (f)^, ... in the Ising theory), tuning a minimal 
set of adjustable non universal parameters [27] [28]. This description is valid because the 
introduction of the irrelevant operators leads simply to a multiplicative renormalization of 
this set of free parameters. Accordance with experiments is under investigations [38] [39] 
[40] [41] [42]. For simple liquid it seems to work, however for complex fluid the field theoretic 
methods is questionable.^ 

In this chapter we will investigate the universal crossover behavior that is observed for a 
family of scalar models in the weakly coupled limit. We do that because similar interactions 
will be obtained considering the effective interactions for the zero mode in large- expansions 
affected by infrared singularities that are the main topic of this work. In order to understand 



^For a review on the argument with a complete hst of theoretical and experimental references see [26]. 
^Instead of "scaling region" we could use also "critical region" . 

^In condensed matter system (contrarily to field theory) the cut-ofF A"-"^ has a physical origin related to 
some size of the basic constituents. In simple fluid we have A « 1 and one can reach the condition ^ 1 

that is necessary in order to apply RG. In complex fluid A^^ is not simply related to inter-molecular distance 
(think for instance to a polymer solution) and the previous condition could not be satisfied [42]. Maybe 
the understanding of these problems (for which techniques presented in this work could be useful) deserves 
future investigations. 
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the nature of the problem, let us consider for instance the following interaction 



(2.1) 



Tuning both m — > or i = r — Tc — > in (2.1), the two universal behaviour (Ising for r — rc, 
Mean Field for u — 0) enter into competition; keeping a proper scaling variable x = tu""' 
fixed one expect to compute the observables (in particular critical exponents) as functions 
of X 



and to recover Ising (Mean Field) behavior by taking x ^ {x ^ oo). The previous 
statements have been demonstrated in a series of theoretical works [29] [30]. More interesting 
it has been shown that the crossover functions fo arc independent of the regularization used 
(i.e. they are universal). This means that in principle one can compute this functions using 
e.g. field theory [29], spin system with medium range interaction [32] [35] [36] [37], or polymer 
systems [35]. 

In this chapter we consider more general interactions than (2.1) or that considered in 
[29] [30]. The reason is that in our large- computation we will obtain an effective action 
for the zero mode without Z2 symmetry [cp — >■ —0 in eq. (2.1)], so that we need to include 
also odd-legs vertices and higher than four-legs vertices. Our goals will be to show that 
including proper counterterms, in the limit in which our crossover is defined, the crossover 
functions remain unchanged with respect to the symmetric case. In some sense we recover 
algebraically the Z2 symmetry that is then dynamically broken."^ This will permit us to 
investigate the nature of the Ising Mean-Field crossover in 0{N) models with a critical zero 
mode (chapter 3). The presentation given follows [56]. In section 2.2 we will also consider 
the case in which the four- and five-legs vertices are tuned to zero at the critical point so that 
one needs to generalize (2.1) including also a 0^ operator. This will permit us to describe 
(sec chapter 4) the crossover in 0{N) models with multicritical zero mode which effective 
theory will be described by an interaction similar to that presented in sec. 2.2. We stress that 
the basic point that permits us to apply the results presented in this chapter to our large- iV 
models, is that the appearance of a zero mode [i.e. r — > in eq. (2.2)] is accompanied by 
the relations V^*^^^(0, 0, 0) = (and similar relations, involving higher order vertices, in the 
case of multicritical zero mode). 

2.1 Critical crossover limit 

We wish now to discuss the critical behavior of the following generic Hamiltonian 



{0)n ~ foitu-'^) 




p 




V^^^ (p, q, -p - ci)ip{p)(p{q)(p{-p - q) 



3! 




u 



[ [ [ V^'^ (p, q, s, -p - q - s)(^(p)<^(q)<^(s)<^(-p - q - s) (2.2) 



"^This is similar to what happen when one wants to investigate chiral symmetry using Wilson fermion. 
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in some particular range of parameters that we will specify in the following. We use a square 
lattice to regularize (2.2), however most of the considerations reported in the following 
easily apply also to sharp-cut off regularization (used in chiral models) simply replacing the 
integration over the first Brillouin zone with an integration over p < A. On the other hand 
the magnetic and thermal regularization hc{u) rc{u) depend of the regularization chosen. 
One is free to impose both V"*^^'*(0, 0, 0) = and K(p) ^ normalizing and translating 
9? by a constant factor ^^(p) — >■ ^^^(p) + kd{p) while V*^'')(0, 0, 0, 0) = 1 can be obtained 
normalizing u. In the models to which we have applied our formalism, H and r are functions 
of the two tunable parameter (for instance p and T in Heisenberg models introduced in 
chapter 1, or M and T in Yukawa models chap. 5), while the coupling parameter u goes to 
zero like in the spherical limit. Simple considerations show that for u ^ Hamiltonian 
(2.2) exhibits a nontrivial limit only if both H and r go to zero. 

Let us first review the case ^(^^(O, 0, 0) = 0. There is an interesting critical limit, the so- 
called Critical Crossover Limit (CCL) [27] [28] [29]. Indeed, one can show that, by properly 
defining a function rc{u), in the limit m — > 0, t = r — rc{u) — > 0, and if ^ at fixed H/u 
and t/u, the n-point susceptibility x„, i.e. the zero-momentum n-point connected correlation 
function 

Xn — ('/-' (0)) connected 

has the following scaling form 

Xn^ut-yr'^{t/u,H/u). (2.3) 

The scaling function f^^"^{x,y) is universal, i.e. it does not depend on the explicit form of 
K(p) and of V'^'^^p, q, r, s) as long as the normalization conditions fixed above are satisfied. 
The definition ofrc{u) has been discussed in Ref. [29] for the continuum model and in Ref. [32] 
for the more complex case of medium-range models. Developing the standard expansion in 
u one can compute rdu) which is obtained requiring that X2 scales according to eq. (2.3) in 
the critical crossover limit. If rdu) — riu + 0{u^), then at one loop we have 



_ 1 _ M 



, l^(^)(0,0,p,-p) 
ri+ ' 



K{p)+t 



+ 0{u^). (2.4) 



For t ^ 



/p K{p)+t JpP^ + t ■ 47r 

where p^ = 4sin^(pi/2) 4sin^(p2/2). Thus, we obtain 



yW(0,0.p.-p) f 1 1 , ,^ , 

~ ' constant mi constant, (2.5) 



_ 1 _ M 



1 , 1 , ^ 

ri m M — — m — h constant 



+ 0{u^). (2.6) 



Therefore, if we define 



u 

r,{u) = —\nu, (2.7) 
an 

the susceptibility X2 scales according to eq. (2.3). Actually one can show that (2.7) is 
sufficient to regularize the theory to all order in the perturbation theory. Indeed in the 



30 



Crossover in critical phenomena 



continuum two-dimensional theory there is only one primitively divergent graph, the one- 
loop tadpole, and therefore only a one-loop mass counterterm is needed to make the theory 
finite. However the exact knowledge (to all order) of rc{u) is a non-perturbative problem 
and leaves some ambiguity in the definition. This will be discussed in sec. 2.1.4. 

In this work we will consider the more general case in which also the three-leg vertex is 
present in (2.2). We will show that, if one properly defines functions rc{u) and hc{u), then 
in the limit m — >• 0, t = r — rc{u) — > 0, and h = H — hdu) — at fixed h/u and t/u, the 
n-point susceptibility Xn has the scaling form (2.3) with h replacing H and with the same 
scafing functions /^^™°^ of the symmetric case. The previous statement easily follows from 
the condition ^'^^^(O, 0, 0) = 0. Indeed considering the continuum fimit of the action (2.2) 

Kff = y"d^(^(x)(n + r)(^(x) + ^(^(x)2n(^(x) + ^(^(x)^ (2.8) 

and making the change of variable x = yu^^^'^, defining 4>{x) = ip{x/ \/u)^ we obtain 

In the previous equation (2.9) it is clear how in the critical crossover limit {u — > 0) the 
three legs vertex can be neglected, so that the critical crossover coincide with that of the 
symmetric theory (2.3) . However in the previous naive discussion we do not have taken into 
account divergences that come from this contribution. This divergences must be carefully 
taken into account by proper counter-terms that define hc{u) introduced above. 

We will first prove this result at two loops and then we will give a general argument 
that applies to all perturbative orders. Notice that it is enough to consider the case h — 
Indeed, if eq. (2.3) is valid for /i = and any n, then 

m=0 m=0 

- ut--Y-^fn;z{t/umM-"'ih/ur. (2.10) 

ml 

m=0 

which proves eq. (2.3) for all values of h. 



2.1.1 Explicit two-loop perturbative calculation 

At tree level, all contributions to Xn that contain three-leg vertices vanish because V^^^ (0, 0, 0) = 
0. Therefore, Xn is identical to Xn in the Z2-symmetric ip'^ theory; in particular X2n+i = 0. 
At one-loop order, graphs contributing to Xn ^-re formed by a single loop made of a three-leg 
vertices and of b four-leg vertices, with a + 2b — n. Each of them contributes a term of the 
form 

where Vs{p) = 'K*^^^(0, p, — p) and V4{p) = V(^)(0, 0, p, — p). The leading contribution for 
i — > is obtained by replacing each quantity with its small-p behavior, i.e. V^lp) ~ 1, 
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Figure 2.1: The four topologies appearing at two loops. Dots indicate parts of the graphs 
with additional external legs. 



^3(p) ~ P^) and K{p) p"^, with — 4sin^(pi/2) + 4sin^(p2/2). Then, we obtain a 
contribution proportional to 



Now, for t — > we have 



2.a (I for 6 = 



(2.12) 



Int for 6=1 (2.13) 



Therefore, for t — > 0, m — > at fixed t/u, we have 

^n/2-l for & = 

t'^Xn/u ~ { int for 6 = 1 (2.14) 

^n/2-b ^ ya/2 for 6 > 2. 

Thus, all contributions vanish except those with: (i) n = 1, a = 1, 6 = 0; (ii) n = 2, 
a = 2, b = 0; (iii) n even, b = n/2, a = 0. Contributions (iii) are those that appear in the 
standard theory without (p^ interaction. Let us now show that contributions (i) and (ii) can 
be eliminated by redefining rc{u) and hdu). Consider first Xi- At one loop we have 

i^,^J}^.J-l ^YM_ + 0i^u). (2.16) 
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For t — > we have 



Thus, if we define 



/_1«EL= /^ + 0(tln(). (2.16) 



then txi/u ~ 0(tlnt/y^, y^) ^ in the critical crossover hmit. 
Now, let us consider the two-point function. At one loop wc have 



t rJu) 1 



-X2 ~ t; 

u u u 2 



K{p) + t {K{p) + ty 



(2.18) 



The first one-loop term is the contribution of the tadpole that has to be considered in the 
pure theory and which requires an appropriate subtraction to scale correctly, cf. eq. (2.7). 
The second one is due to the three-leg vertex and is finite for t — > 0. Therefore, if we define 

rc(^^) = r7-M + ^|M!, (2.19) 

where r^y™™(-u) is given by eq. (2.7), we cancel all contributions of the ip^ vertex. 

Now, let us repeat the same discussion at two loops, in order to understand the general 
mechanism. There are four different topologies that must be considered, see Fig. 2.1. 

Topology (a). 

A contribution to Xn is proportional to 

~ ip X ""^ (g(p) + t)°^-+'(j^(q) + ■ 

where 2a -\- 2h -\- c -\- d — n. The leading contribution for t — > is obtained by setting 
1/W(p, -p,q, -q) 1/(^)(0, 0,0,0) = 1. Then, the integral factorizes and wc can use 
eq. (2.13). Ignoring logarithmic terms, we see that the corresponding contribution to f^Xn/u 
scales as 

'«A"/%n/2-a-.^/«y/%(e+.)/2_ ^2.21) 



tJ \t. 

Thus, a non- vanishing contribution is obtained only for c-\- d — Q. Three-leg vertices can be 
neglected in the critical crossover limit. 
Topology (h). 

A contribution to Xn has the form 

yn/2+l 



^ VA{vYVA{cifVA{p + q)-l^3(p)'^^3(q)-^3(p + q)^ ..^ 
(ir(p)+i)«+<^+i(X(q)+i)''+-+i(ir(p + q)+t)-+/+i' ^ ■ ' 

where 2{a -\- h -\- c) -\- d -\- e -\- f — n — 2. The leading infrared contribution is obtained by 
approximating all expressions with their small-p behavior. Therefore, we can write for i — > 

I I {py{<?n^qy ^223) 
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Integrals of this type can be easily evaluated. By writing = {p"^ + t) —t in the numerator, 
we obtain integrals 

Imnr — I I 2 ' (2.24) 

with m > 0, n > 0, r > 0. Then, we can rescale p t^^'^P, q t^^^q and extend 
the integration over all R^. This is possible since the corresponding continuum integral is 
finite. As a consequence the integral scales as ^^-m-n-r- ^j^g corresponding contribution 
to t^Xn/u scales as 

(^|y/%n/2-a-6-c-l ^ (^^y/%(.+e+/)/2_ (2.25) 

Thus, a non vanishing contribution is obtained only for d + e + / = 0. Three-leg vertices can 
be neglected in the critical crossover limit. 

Topology (c). 

A contribution to Xn has the form 



u 



n/2+l 



[ I vW(0,p,q,-p-q)\/(3)(p,q,-p-q) 

^ VA{pYVM)'VAiv + q)^V^3(p)'^V^3(q)^V^3(p + q)^ , . 

where 2(a-|-6-|-c)-|-(i-|-e-|-/ = n—1 and we assume without loss of generality a > b > c. We can 
now repeat the analysis performed for topology (b). We replace each quantity with its small- 

momentum behavior. In particular, we can replace V^^\p. q, — p — q) with (p^ + g^+p + g ). 
Then, we rewrite each contribution in terms of the integrals Imnr, cf. eq. (2.24). However, 
in this case it is possible that one (and only one) of the indices vanishes. If this the case, 
Imnr factorizes and we can use the one- loop result (2.13). A careful analysis shows that in 
all cases the integral scales as t~°'~^~'^ for t — > except when b = c — Q. In this case, if a 7^ 
the integral scales as t~"lnt, while for a = it scales as log^t. Thus, ignoring logarithms all 
integrals scale as t~"'~''~'^. Therefore, the corresponding contribution to t^Xn/u scales as 

These contributions always vanish. 
Topology (d). 

A contribution to Xn has the form 

! /fyWfn o n n\? ^-i(P)"^4(q)^V^4(p + q)-V:3(p)^V^3(q)^V^3(p + q)^ 
ypiq (i^(p)+t)«+'='+i(i^(q)+t)''+-+i(i^(p + q) + t)<=+^+^' 

(2.28) 

where 2(a-|-6-|-c)-|-d-|-e-|-/ = n, a>b>c. We repeat the analysis done for topology (b) 
and (c). We find that the integral scales as 

for a = 6 = c = 0; 
^-a-h-c+i I 1/t for 6 = c = 0, a > 1; (2 29) 

^ for c = 0, 6 = 1, a > 1; ^ ' ^ 

otherwise. 
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Thus, if 6 7^ and c ^ 0, ignoring logarithms, the contribution to f^Xn/u scales as 

(^«y/%n/2-a-.-c+l ^ ^]!_Y'\(a^e^m^. (2.30) 

These contributions therefore always vanish in the critical crossover limit. If, however, 
6 = c = 0, then 



tJ \t. 

and thus one may have a finite (or a logarithmically divergent if a = 0) contribution for 
d = e = f = 0. Let us now focus on this last case in which a = n/2 > with n even. Let us 
now show that these contributions are cancelled by the one-loop counterterm due to rc{u). 
Indeed, at two loops we obtain a contribution to Xn the form 

where 2a + b — n. Thus, contributions to f^Xn/u scale as — t^f^ and thus vanish 

unless 6 = 0. Thus, for each even n there are two contributions that should be considered: 
one with topology (d) and one associated with the counterterm rc{u). Taking properly into 
account the combinatorial factors, their sum is given by (of course, we only consider here 
the contribution to rc{u) due to the three-leg vertices) 



u' 



n/2+l r 2 



\/(''^)(p,q,-p-q)2 V'-^\0,q,-q) 



{K{ci) + t){K{p + ci) + t) {K{ci) + ty 



(2.33) 



For a > 1, the subtracted term improves the infrared behavior, and indeed the integral scales 
as Int and is therefore irrelevant in the critical crossover limit. 



2.1.2 Higher powers of the fields 

It is interesting to consider also Hamiltonians with higher powers of the field (p. The standard 
scaling argument indicates that these additional terms are irrelevant for the critical behavior, 
but in principle they can contribute to the renormalization constants like he and Vc- For this 
purpose let us suppose that the Hamiltonian contains also terms of the form 

J,k-2)/2 r r 

= j:^/ ■■■/ <^(Pi + ---Pfe)^^'HPi,---,Pfe)<^(Pi)---<^(Pfe)- (2.34) 

Note the particular dependence on the coupling constant which is motivated by the large- 
calculation and is crucial in the argument reported below. 

Let us consider the contributions of these additional terms. At tree level we have an 
additional contribution to Xn given by 



y(n-2)/2 
^Xn^—^V^-\0,...,0) 



(2.35) 



2. 1 Critical crossover limit 



35 



Thus, t"Ax„/-u ~ u'^'' "^^Z^, that vanishes as ^ for k > 4. At one-loop order there are 
additional contributions of the form 

j (X(p) +t)-^«^y3(p)"^ • • •K+2(p)""+% (2.36) 

where Vjk(p) = V^^\—-p, p, 0, . . . , 0) and ~ 2)ajk = n. Proceeding as in the previous 

section, using eq. (2.13), and noting that 'Ylik>4^k > by hypothesis, we obtain 

r AXn/w ~ M"/2-^t^-^^>3 ~ ^^"^+5 T,k>^{k-4)ak _ (2.37) 

Here, possible logarithmic terms have been neglected. Then, since some with /c > 5 is 
non-vanishing by hypothesis, we find that this correction vanishes in the crossover limit. 
Therefore, no contribution survives at one loop. The same is expected at any perturbative 
order. 



2.1.3 The general argument 

The discussion reported above shows that at two loops one can define a critical crossover 
limit with crossover functions that are identical to those of the symmetric theory. This is 
expected to be a general result since formally the added interaction is irrelevant. This result 
can be understood diagrammatically. 

Consider the continuum theory with Hamiltonian 

^ = / ^ + l^^^^f + ^^i^f^^i^) + Im'- (2.38) 

Given an /-loop diagram contributing to the zero-momentum n-point irreducible correlation 
function, we can compute the superficial degree of divergence of Feynman integrals ~ A"^ 
as 

D ~ KM+^N^-2N^ (2.39) 

obtained rescaling each momenta with A. In (2.39) £ is the number of loop of the diagram. 
As (A4) is the number of the three (four) legs vertices and Aj is the number of the internal 
lines. Using the topological relations 

n + 2Ni = 3A"3 + 4A"4 
Ar3 + 2A"4 = n + 2£-2 

we obtain for d = 2 

D ~ a2(i-^4) (2.40) 

Thus, there are primitively divergent diagrams for any n: those with A4 = (and corre- 
spondingly A3 = n -|- 2/ — 2) are quadratically divergent, while those with A4 = 1 (and 
A3 = n -|- 2/ — 4) are logarithmically divergent. The previous considerations suggest that 
in the evaluation of the generic Feynman diagram with n external legs A^ and log A^ diver- 
gences will be generated. To regularize the theory a counterterm 5„ • 0(p)"' is introduced. 
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However in principle this could not be sufficient. Indeed we observe that other divergences 
comes from second derivatives of the n-point irreducible correlation functions with respect 
to the external momenta 4>{p)'^~^Op(f>[p) . Following the same step done before we find 

□p(0(p = 0)<P{Or-') - (2.41) 

which diverge logarithmically for every n when A^4 = 0, so that also these operators must 
be included in the renormalized Hamiltonian which will contain an infinite number of coun- 
terterms 

+ J2[Zn{A)r + Cn(A)0"-'n0]- (2.42) 

n 

Now, let us show that all counterterms except those computed in the previous section can 
be neglected in the critical crossover limit. Suppose that we wish to compute Zn{A) at / 
loops. Keeping into account that the divergence may be quadratic or logarithmic, we expect 
the Z-loop divergent contribution to Xn to be of the form 

n/2+l-l 

^^^(aiA^ + iPi[ln(AVi)] + P2[H^Vm. (2.43) 

where Pi{x) and P2ix) are polynomials and ai a constant. Therefore, the contribution to 
f^Xn/u vanishes unless n/2 + I — 2 < 0. The only two cases satisfying this condition (of 
course n > 1 and / > 1) are / = n = 1, / = 1 and n = 2, which are the cases considered 
before. Let us now consider the contributions to Xn,i, which is the ffist derivative of the 
n-point connected correlation function with respect to the square of an external momentum 
p computed at zero momentum. Since momenta scale as t^^^, in the critical crossover limit 
we should have Xn,i ~ '>J't~"'~^fn(t/u,H/u). The divergent contributions are logarithmic 
(diagrams with A^4 = 0) and therefore we expect an Z-loop contribution of the form 

— (P[ln(A7t)] + finite terms). (2.44) 

t 

Considering t"-~^^Xn,i/u, we see that this contribution always vanishes in the critical crossover 
limit. Therefore, the renormalization constants Ci(A) can be neglected. Thus, the only renor- 
malizations needed are those that we have considered. Finally, let us show that correlation 
functions computed in the renormalized theory have the correct scaling behavior. Indeed, 
in the renormalized theory diagrams scale canonically with possible logarithmic corrections. 
Therefore, D introduced at the beginning of this section scales as u^*'^^^^^t^~^*x\ogs, so 
that the contribution of D to f^Xn/ u scales as 

X i X „^4+iV3/2^i-Ar4 ^ ^iV3/2_ (2.45) 

Therefore, the only non- vanishing diagrams have A^3 = 0, confirming the claim that three-leg 
vertices do not play any role. 

2.1.4 A unique definition for the renormalization functions hdu) 

and rc{u) 

In this section we wish to discuss again the definition of rc{u) and hc{u). It is obvious 
that these functions are not uniquely defined, since one can add a term proportional to u 
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without modifying the scahng behavior. We wish now to fix this ambiguity by requiring that 
t = h = corresponds to the critical point. 

It is easy to see that no modifications are needed for hc{u). Indeed, with the choice (2.17) 
one obtains the correlation functions of the symmetric theory and in this case the critical 
point is uniquely defined by = by symmetry. The proper definition of rc(-u) requires 
more care, since we must perform a non-perturbative calculation in order to identify the 
critical point. For this purpose we will use the fact that in the critical crossover limit the 
perturbative expansion in powers of u is equivalent to the perturbative expansion in the 
continuum 0^ theory once a proper mass renormalization is performed. 

In the continuum theory, if t = tcont/'^cont is the a-dimensional reduced temperature 
defined so that i — corresponds to the critical point, we have at one loop, cf. eq. (2.10) of 
Ref. [32], 

iicontX2,cont ^ l + TTl^ ( + ^ + 87ri?2 ) + ^(r^ In?), (2.46) 



t Snt^ 

where D2 is a non-perturbative constant that can be expressed in terms of renormalization- 
group functions, cf. cq. (2.11) of Ref. [32]. By using the four-loop perturbative results of 
Ref. [86], Ref. [32] obtained the estimate D2 = —0.0524(2). It is not clear whether the error 
can really be trusted, since in two dimensions the resummation of the perturbative expansion 
is not well behaved due to non-analyticities of the renormalization-group functions at the 
fixed point [87, 88]; still, the estimate should provide the correct order of magnitude. 

The expansion (2.46) should be compared with the perturbative expansion of X2 in the 
lattice model. We write rdu) as 



K{p) p 



1 , 2567r 3 
A = -D2 - — In 



87r 3 87r 2 



n2 



(2.47) 



where A has been determined comparing X2 with a lattice regularisation eq. (2.18) 



u (. t „ . . /■ ri^4(p) 1 



"^^ = i + 8;^V""32^"^"'^"'" 



0{u^). (2.48) 



with X2 in the continuum limit given by eq. (2.46). If we use definition (2.47), the critical 
point corresponds to t = 0. 



2.1.5 Higher than two dimension 

In this section we want to generalize the previous results (obtained for d = 2) to higher 
dimension 2 < d < 4, for which Ising criticality is still present. Cosidering the continuum 
action (2.49) one is able to predict the general scaling relations (2.52) and (2.53). On the 
other hand the critical parameters hc{u) and rc{u) will include more contributions with 
respect to the two dimensional case [56]. This is simple related to the fact that in three 
dimension for instance, the mass counterterm must be computed up to two loops. However 
in the case in which the theory is not symmetric, anomalous terms appear also for ^3; these 
bad contributions cannot be deleted by using hc{u) or rc{u). In this section we show that the 
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expected scaling relation (2.53) still survive if the critical field is translated by a constant 
factor^ (i.e. ip ^ (p + u^^'^kji), so that the anomalous three legs diagrams can be deleted 
fixing kR. 

Let we start considering the continuum interaction in d < 4 dimension 



«SeontM = / d'r \h^{v) + \ [d^{v))\ '-^{vf + ^^{vfn^{v) + ^<^(r)^ 



, (2.49) 



generalizing the considerations done in (2.8), i.e. making a change of variable r = yu ^^^^ 
and defining a new bosonic field 

V'(x) = ii(<^-2)/PM(^(xx,-V(4-ci)). (2.50) 

we can rewrite (2.49) as 



-"cont [ 



where we have defined. 



, (2.51) 



H = //t,-(<i+2)/P(4-<i)]^ f = ru-^m-^l (2.52) 

Thus, formally, once the action is expressed in terms of ip, the bare parameters appear only 
in the combinations H and f while the three fields coupling can be neglected. Then, consider 
the zero- momentum connected correlation function Xn, neglecting the three- legs vertex, we 
have 

Xn = j c^S . . . d'r^ (^(0)^(r2) . . . (^(r„))-- 

^ ^[2d-n(2+d)]/[2(4-<i)] J ^dg^ _ _ _ ^dg^ (^(0)^(S2) . . . ^(s^))™"-^ 

= «Pd-n(2+d)]/[2(4-d)]_^^(^^~)^ (2.53) 

i.e. xi~P''-"(2+'^)]/[2(4-d)]^^ jg ^ scaling function that is universal given the normalization 
conditions for H and f. The above-reported discussion is valid only at the formal level 
because we have not regularized the theory, or in other words a cut-off regularization breaks 
scale invariance. 

Let us first discuss the symmetric case because it will be used in this work. In a similar 
way of what happens for d — 2, the previous relation (2.53) remains true if proper countert- 
erms are included in the theory. For d < 4 we have to include only a mass counterterm rc{u) 
[hc{u) = due to the symmetry]. Having done this the correlation functions Xn satisfy the 
scaling relations (2.53) with i = tu~'^^^^~'^\ t = r—rc{u) replacing f, and h = /i,ii~('^+2)/[2(4-rf)]^ 
h = H — hc{u) replacing H 

Xn = iiPd-n(2+d)]/[2(4-d)]_^^(^^ (2.54) 



^The normalization of kR by using u^/^, has been taken for convenience. 
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However in more than two dimensions more than one loop terms enter into the definition of 
rc{u) and hc{u). This is related to the fact that in the continuum theory, divergent diagrams 
are present beyond the one loop order. For instance in three dimension {hc{u) = 0) we have^ 

ri^){u) ^-^u+-^\nu + Ku' (2.55) 

where K is chosen so that t = corresponds to the critical point (see sec. 2.1.4). 

The non-symmetric case demands a more accurate analysis. Following the discussion of 
sec. 2.1.3 we have that a generic contribution to Xn scales as^ Xn ~ G{A, i)i~"ii"/^+'~^. Using 
the scahng relation (2.54) t ~ xi2/(4-d) observe that 

n(2+d)-2d „_4 , . 

u 2(4-d) ^ u^+^G{A,u^) (2.56) 

so that the anomalous diagram can be obtained by studying + i < 0. If = 3 there are 
the following dangerous contributions: a) n = 1, £ = 1, 2, 3 (so that he can be obtained by a 
three loops computation), h) n — 2 i — 1,2 [cancelled in the definition of rc{u)], but also c) 
n = 3 ^ = 1 that in principle give us some problems. Now we show that the one loop n — 3 
anomalous contributions can be deleted by introducing a new bosonic field 0(p)^ 

<^(p) = (/.(p) - «^/2M(p) (2.57) 



and taking kji so that the effective three-legs vertex [computed up to 0{u) in c? = 3] for 

the new bosoni 
fields, we have 



the new bosonic fields is zero. Defining V (pi, • • -p^) the n— legs vertex for the translated 



^^'\pi,P2,P3) = V^^\pi,p2,P3) -uknV^'\p,,p2,P3,0) (2.58) 

V'^'\P1,P2) = V^'\p^,P2)-ukRV(-'\p,,P2,0) 

+!^y(4)(p^,p^,0,0) (2.59) 

-^^^V^(^)(0, 0,0,0) (2.60) 
Imposing that the first radiative correction for the three-legs vertex is null we get the result^ 



kR = 



VJO) 



1 f ,3Vsip)V,ip) 1 (-.Vsip 



2 J ^ k{py 6 J ^ k{p) 



\3 



(2.61) 



^The form of Tc and he depends on the explicit regularization used. Here we use a cut-off because it will 

be used for a fermionic model investigated in this work. 

''G„(A) for d = 3 computed up to £ loops diverges as A^+^. Indeed generalizing (2.40) we have that for a 
generic £ loops diagram D ~ A^+2(i-Ar4) = f^3-N,+Ns/2-n/2 

^This operation is exactly what we do when we obtain the effective action for the zero mode in the large 
N limit where in principle the condition V^^\0, 0, 0) = is not satisfied. For d = 3 one needs to take into 

account also radiative corrections. 

^Notice that both the integrals are infra red finite in three dimensions. 
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that delete the the anomalous c) contributions. We have regularized the anomalous a) and b) 
contributions introducing respectively hc{u) and rc{u) similarly to what done for o? = 2 (but 
with more diagram). But now new anomalous contributions appear due to the new terms 
(proportional to ku) in (2.59) (2.60). It is straightforward to understand that this terms can 

be deleted recomputing hc{u) and rc{u) as a function of the new vertices V . Indeed this 
is possible because the leading behavior of the vertices remain the same Vn{0) ~ K.(0). 



2.2 Multicritical crossover limit 

Now we want to generalize the previous considerations, considering more general interactions 
than (2.2). We want to take into account interactions in which we tune to zero at the same 
time all Vj{0) for j — 0,1- ■ ■N'+2 (with jV odd), while VV+3(0) remain finite and positive.^" 
In the continuum limit we want to study the following Hamiltonian 



Af+3 



n=2 



(2.62) 



with aA/'+2 = 0. Here we have disregarded the momentum dependence of the vertices — it 
is formally irrelevant in the continuum limit — and higher-order terms in the kinetic term. 
We are interested in the limit in which — > for i = 2, ■ ■ -M + I, so that the leading 
non- vanishing term is aj^+sX^^^- In analogy to what have been done in the previous section 
2.1 we wish now to compute the crossover from mean-field to multicritical behavior. The 
analysis is quite complex and thus we only consider explicitly the case J\f — 3. 
Using the lattice regularization we start from the following Hamiltonian 



n 



Uoi(, 



(0) + ^(i^(p) 



%2)0(P)0(-P) 




U03 + ul/'V('\{p,}) 

P. 1 v('\{p,}) 



+1 [til ^^p^))^(tp^) [l+^^'Hip.})] 



where we have defined 



y(')(o,...,o) = 



(2.63) 



(2.64) 



for i = 1, • • • , 6. We have written each vertex as a constant plus a function of the momenta 
multiplied by a power of Uq. This is indeed the expression that one would obtain from the 



10 



The positive condition guarantees us that the stability of the theory. 
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effective Hamiltonian for tlie zero mode by setting uq ~ l/N"^. Moreover, we liave dropped 
the contributions of vertices witli more tlian six fields. We will discuss them in sec. 2.2.4. 
The normalization of the field is fixed by requiring: 

K{p) « p' + 0{p'), (2.65) 

while we introduce ccj to name the leading contribution in p for higher order than two vertices 

y»(p,-p,0,---0) « a,p2 + 0(p^). (2.66) 

We wish now to show that four critical parameters {ud, i = 1, . . .4) can be defined, such 
that for Ui = UQi — Ud —>■ 0, {i = 1, . . .4), uq —>■ 0, the n-point susceptibility has the scaling 
form: 

Xn^t^~''fn{xi,X3,X4,X(i), (2.67) 

where we write for notational convenience t = U2 and we set Xi = Ui/t. Moreover, we 
will show that the function /„ is independent of the vertices y^'^ and it is thus a universal 
scaling function of the theory. The non-universal features of the model appear only in 
the functions Ud-, which can be computed in perturbation theory up to an additive constant. 
In our case a two-loop calculation provides the full answer. Notice that, as pointed out in 
section 2.1 for the N = 1 case, in order to prove eq. (2.67) it is enough to consider the case 
in which ui — Q. This is the case we shall consider below. 



2.2.1 One-loop computation 

We begin by considering a generic one-loop diagram with three-leg vertices and so on up 
to Uq six-leg vertices contributing to Xn- For t ^ its contribution is 

Xn - + ^^V^(p)r (2-68) 

where Vi(p) = V^^\p, — p, 0, . . .), Wj = Uq~'^^^^ {j = 3,4,5,6), and — 0. Expanding we 
obtain contributions of the form 

where < ki < Hi and = n^- If X](?T'j — kj) > 2 the integral behaves as with 
A = —1 + X^(nj — /cj), so that 

r-^Xn~n<'"''^^- (2.70) 

i 

Thus, in the scaling limit, the contributions of the formally-irrelevant vertices can be ne- 
glected and thus we can set ki = 0. The remaining contributions (with ks = k4 = k^ = = 
0) scale according to eq. (2.67). Note that in this case the integral can be computed in the 
continuum, by replacing K{p) with and by integrating over the whole plane. Therefore, 
the result is independent of the explicit form of K{p). 
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W^{ni — ki) = 1, we have similarly 

i"-'Xn~^r"'^«^flogt. (2.71) 

The formally-irrelevant vertices can still be neglected. However, in this case there is an 
additional logt that breaks the scaling law for ki = 0, i.e. for 11,3 + + Uq = 1. Anomalous 
contributions appear therefore in xi, X2, and X4- 

Finally, if '^{rii — ki) — (this implies ris — k^, — k4^, rie — ke) the integral is constant 
as i — > 0. Therefore, we have 

t-'Xn-jU^-'^. (2.72) 

i 

In the scaling limit this contribution is irrelevant except when n < 4. 

In conclusion, at one loop, the only contributions that scale anomalously in Xn are that 
whith n < 4. We shall now show that we can absorb these anomalous terms in the definition 
of the functions (in this section contributions to Ud given by one-loop diagrams will be 
indicated as wj'). 

Let us begin by considering the one- loop expression of xi' 



Xi 



11 



2tL K{p)+t 



t 2t 

with 



and = 4 sin^(pi/2). Thus, if we define 

= C/1 + C/2 



the anomalous contributions in xi cancel out and the result is independent of V'^^'> and of 
K(p) 

x}' = -|^log(32x6). (2.76) 

In order to render easier the analysis of the theory at two loop level, in (2.75) we have 
introduced Ui and t/2, defined as 

C/i = ^(log«6-47rXi) (2.77) 

OTT 
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Let us now consider the susceptibility with n = 2. We have 



"c2 



U 



c2 



1 

2t 



Thus, if we define 



Uj + »'iV'i(p) 



^c2 — 



1 

+ 



b'3+»'3l;i(p)]- 



(2.79) 



MS 



+ 



u, 



1/2 



2i 



^3(P)' V^ip) 



K{py K{p) 



+ 0(«;/Mogt) 



D1+D2 + L>3 

M4 

Stt 

1/2 



M, 



(log Its - 47rXi) 

V^3(P)' 



,1/2 



^4(P) 



we obtain 



11 



X 



-J-|ilog(32x6), 

on on 



(2.80) 



(2.81) 



that is independent of the vertices and of K{p). 
Now let us consider xs- At one loop we have 



u 



[m3 + m;3V3(p)]^ 



+ 



2t 



[K{p) + 1]3 

[Ms + W3\/3(p)] [Ui + W4V4(p)] 



[^(p)+tP 



W5t^5(p) 

2t J^K{p)+t 



u 



11 

c3 



3X3X4 



3/4 



8n 



2t 



2V,{pf 31^3(P)V^4(P) , V^5(P 



K{pf 



K{py 



+ 



K{p) 



+ 0{uf\ogt). 

(2.82) 



Thus, the correct scaling is obtained by setting 



u 



11 — 



c3 



= T1+T2 



To = 



3«fi3/4 



1*6 



3/4 



^3 

21^3(P)^ 
ir(p)3 

V^3(P)V^4(P) 

K{pr 

pK{p) 



(2.83) 



Finally, let us consider the four-point function. There are five graphs contributing to it. It is 
easy to see that an anomalous contribution arises only from the insertion of the 6-lcg vertex. 
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(I) 



q 




P 



(II) Vi 




p + q 



Figure 2.2: The general two loop topologies for Xn- We report only the vertex {vi and V2) 
for which more than two internal legs are contracted. 



However, if we wish that X4 does not depend on any lattice detail we should also subtract 
other contributions that arise from the other graphs. We define therefore 



2.2.2 Two loop computation 

At two loops there many graphs contributing to Xn- In general they can be grouped in 
two different topologies reported in Fig. 2.2. Graphs are obtained from these topologies by 
adding external legs to vertices vi and V2 and to the internal lines. We will then distinguish 
the contributions in different classes, according to the nature of the vertex appearing in Vi 
and V2- We will thus consider contributions (lA) and (IB): in the first case we consider at 
the vertex v the zero-momentum part of the vertex, while in the opposite case we consider 
the remaining part. Analogously we distinguish contributions (HA), (HB), (HC). 




Q1 + Q2 + Q3 + Q4 + Q5 + Qa 




(2.84) 
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Case (lA) 

In this case the contribution to Xn can be written as 



Xn ~ Uj 



,r" f [K{p) + t]-^'^^-'l[[ui + wM{p)r 

f [K{q) + Yliu, + w,Vi{ci)r^ (2.85) 

Expanding we obtain a sum of terms of the form 

Xn ~ ii,-r"J]«,''+'^^<+"^^-*^'-'^' /" [[K{p)+t]-^''^-'Y[[v,{p)p-''' 

i JpJq j 

[K{ci)+t]-^"''-'ll[Vi{ci)r'-''\ (2.86) 

i 

The integral scales as t"^ with A — '^^ih + hi) apart from logarithmic factors. Therefore, 

i""'Xn ~ Xj II x\^+^^wT^'^^~'''-''\ (2.87) 

i 

Thus, the contributions of the additional vertices can be neglected and we should only 
consider the case rii = ki, rrii = hi. In this case, if there are no additional terms proportional 
to logi, these contributions scale according to eq. (2.67). Logarithmic terms occur when 
^rii — or '^rrii = 0, i.e. when there are tadpoles. Thus, an anomalous contribution 
occurs for the terms of the form 

with n = j + ma + 2m4 + Auiq — 4. However, at the same time we must consider the 
contributions of the one-loop counter-terms that give rise to terms of the form 

Xn ~ / [K{q) + Him + WiViiqT^ (2.89) 

^q i 

The analysis presented above can be repeated for these terms here, showing that one can 
neglect the contributions of the vertices V^. Thus, we have a contribution 

Xn ~ ul^t-^l[u, [ [i^(q) (2.90) 

Now, ulj is the sum of several contributions: we only consider here the logarithmic terms 
[Ui, Di and Qi in (2.75,2.80,2.84)]. The additional ones will be relevant for case (IIC) below. 
Taking into account the combinatorial factors we obtain 

Xn~r%n«^ /^ ^J)^^ + ^log^6-i^i j[K{c^)+t]-^^^-' (2.91) 
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The subtracted term replaces the logt term that arises from the integration over p with 
log{t/uQ). Thus, if no other logarithmic terms are generated by the integration over q this 
term scales correctly. However, when rrii — 0, another logt is generated. This occurs for 
n — 2 and j — 6. Thus taking into consideration this new anomalous contribution to X2 
(and neglecting other contributions): 



21;IA 



At 
'At 



47r 



log 2,2xq 



log 32x6 - log uq + Ki 
47r 47r 



(2.92) 



In order to cancel the previous anomalous contribution to X2, ^ correction to Uc2 will be 
introduced in sec. 2.2.3. 



(IB) 

The general contribution scales according to 



Xn ~ Wjt "J]^ 



ki+hi^,ni+mi-ki-hi 



W, 



[ [ V^^\p, 

«/ p «/ q 



-p,q, -q,o, ...) 



mi— hi 



(2.93) 



Keeping into account the fact that ^^'■-'■'(p, — p, q, — q, 0, . . .) ~ (p^ + g^) for p.q ^ 0. wc find 
that, if ki and hi 7^ 0, the integral scales as with A = X^j(A;i + hi) — 1 apart 
from logarithmic factors. Thus 



n 



X, 



hi+hi rii+rai-ki-hi . 



(2.94) 



such a term always vanishes in the critical crossover limit. 

Let us now consider the case in which /cj = but hi ^ Q There are several 
anomalous diagrams that come from this contributions. We are going to show now that all 
these non-scaling terms are eliminated by one loop diagrams with insertion of counter terms 
defined in the one loop computation. We define: 



^^'■^(P, -P, q, -q, 0, • • •) = ^,-(p) + ^i(q) + ^S-(p, q) 



(2.95) 



then inserting (2.95) into (2.93) and following the discussion of the /ij 7^ and /cj 7^ 
case, we find that the anomalous terms only come from Vj{-) in (2.95), while Vj terms 
are regular. At this point the integrals factorize and the discussion reduces to one loop 
computation of the previous section with the properly counter terms introduced. Explicitly 
we have 

^ JJa;i'*'Wi"'+"*'-'^' (2.96) 



f'-'Xr 



so that we have to take care all the diagram for which: 

6 



5:ii^(„.+d.)+i^<i 



(2.97) 
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having defined di = rrii — hi and with j = 4, 5, 6. 

If J = 6 we find a scahng diagram with di = and = 0. This non-universal contribution 
is cancelled by Qq term of ul[ (2.84), indeed for this terms (including combinatorial factors): 



Wj 



Y[ xtwi"^^"^^-^^ log t (2.98) 



that can be neglected in the CCL limit. In a very similar manner using T3 (2.83), we can 
delete the terms with j = 5, = 0, = and j = 5, = 0, = 5j,3; while using D^, 
(2.80) the diagrams with j = 4 and rij = are cancelled (anomalous one are di = 0, di = (5j^3, 
di — 25i^s s-nd di — 5j_4). The last diagrams with j = 5 we have to discuss is the ^3 = 1, 
di — 0. This diagrams are regularized by one loop-insertion of Q5 (2.84). Indeed we have: 

.. TJ^k. f ^3(P)V^5(P) ^3(P)V^5(P) 



{K{p) + ty K{py 

llxi'^nogt (2.99) 



In a very similar way the remaining anomalous diagrams for j = 4: are regularized by T2 
(n3 = 1 and di = 0, = I and di = Si^s), Qi = 1 and di = 0) and Q3 {n^ = 2 and 
Si = 0). In the discussion above we do not have used Q2, Ti and D2, we expect that these 
terms enter into topologies of type (II) . 

In the case in which J2i h — hi = other counterterms must be included in order to 
regularize the theory. In this case (neglecting log): 



n+2 

f'-'Xn - ^Uw-^^"^^ = (2.100) 



so that we have to take care Xi X2 in (IB). From now we take the following short-hand 
notation for the propagator 

A{p)=K{p)+t, (2.101) 

then 

21;IB_ 'tJ'T f ^5(P,q) , uT f ^3(p)^4(p,q) 



/ ^5(p,q) , ur f ^3(p)^4(p,q) . . 

4, A(p)A(q) ^ At ip,, A(p)2A(q) ^ • ' 



8t 

Therefore in the CCL (i.e. for t ^ 0): 

3/4 / , \ 

X.''''^ = ^(^^;iB-%flog|J (2.103) 

,(1) ^ _1 /■ V^5(p,q) If V;(p,q)V:3(p) f YM 

8 ip,, X(p)X(q) ^ 4 KipfKiq) 16n K{q) 

_l f YMf^_l_\,lf V^3(P)V^4(P) f I i\ 

4 ip.d ^(p) U(q) qV K{pr U(q) qV 
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4 4^ i^(p) U(q) A'(q) J ' 

,(1) _ 1 f V^3(P)^4(P) 1 f VM , as /• V^4(P) 



where ctj was defined in (2.66). Now we compute the IB diagrams that enter into the 
computation of X2' 

. 21;IB ^ f VeCP, q) ^5(P, q)V3(p) f ^4(P, q)V'3(p)^ 

8t ip,, A(p)A(q) ^ 2t ip,, A(p)2A(q) 2t 7^,^ A(p)3A(q) 

_«£ /- ^3(p)^3(q)^4(p,q) , u, f V,ip)V,{p, q) 
4t 4^ A(p)2A(q)2 ^ 4t 4, A(p)2A(q) ^ " > 

Therefore in the CCL: 

21;IB _ ^6 ( .{2) ^21;IB . t\ (9 ^(^7^ 

tX2 - y (^•^21;IB-^l0g-J (2.107) 

1 V-6(p,q) l V5(p,q)V^3(p) l V^3(p)^V^4(p,q) 1 V^4(p)V;(p, q) 
8ir(p)X(q) 2 X(p)2K(q) 2 X(p)3X(q) 4 X(p)2K(q) 



•^21:IB 



1 ^3(p)^3(q)V^4(p,q) l 1 f YMf^-l\+lf 

4 X(p)2K(q) J 4 ip,^ K(p) U(q) qV 2 4^ X(p) 

U(p) ^(q) ^(q) ^(q) ^(p) q' ""'qV stt 7^ x(q) 

_1 / ^4(p)/ l"3(p)V 1 1\ V"3(q)^ 1 «3^\ 3^3^ /■y4(q) 
2 4^ X(p) lvX(p)2 U(q) q^; + X(q)2 X(q) q2 J + le^r K{q) 

^3(p)^4(p)/^3(q) aA fVs{q)V,{q) a, f V,{q) 



2yp,q ^(P)' V^(q) qV 87r,/q ;^(q)2 167r^i^(q) 

4 ip,, /^(p)2 U(q) q2y'^4 4qir(p)U(q)^ qV ^' ^ 



;.(2) _ 1 / v^6(p) ^ 1 /- v^5(p) / ^ ^3(P)\ 1 f v,ip) ( , v,{pf\ 

- -4yp^+2yp^i,"^+^j-2yp^i,"^ 

+ 04) (2.109) 



c^3 7 ^3(q)^4(q)' ^ 1 [Va{v)(V,{v) 



2 7q K{qf Aj^K{p)\K{p) 
Case (IIA) 

In this case the contributions of vertices Vi and V2 are simply Uj^Uj^. The contribution to Xn 
has the general form 

(t + ir(p))-i-^"'(t + ir(q))-i-^'"»(t + ir(p + q))-^-^"' 
\[[Vi{p)r-''mqT^-^^[V,{p + q)]^'^-^ (2.110) 




p -'q 
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with the topological relation n = ji + j2 — 6 + J2ei^ ~ 2)(?^£ + me + pi). The analysis of this 
contribution is completely analogous to that presented in sec. 2.1.1 for the two-loop topology 
(b). The integral scales always as t~'^ with A = 1 + J^iih + hi + U) in the infrared hmit. 
Thus, we have 

i"~'Xn ~ Xj.Xj^ ^]^^k,+h,+k^n,+m.+P^-k.-h.-k^ (2.111) 

i 

In the scaling limit, the only non- vanishing contributions are those with rii = ki, rrii = hi, 
Pi = li- the formally irrelevant vertices can be disregarded. The remaining terms scale 
correctly. 

Case (IIB) 

In this case the vertex vi enters into the diagrams with the zero momentum part of a vertex 
Uj^, while V2 gives the contribution V^'^\y>i ~P ~ 0; ■ ■ ■)■ Then Xn has the following 
form 

Xn - n[«^^''^ V'^"^^^''""^^^ /0^-^)(p,q,-p-q,O,...) 

i -'p -'q 

n[T/,(p)]"-'=^[T/,(q)]-^-'^^[T/,(p + q)p-'^ (2.112) 

i 

Here n = ji + ^2 — 6 + Yliii^ ~ + + pe) as before. The analysis follows from the 
analogous one presented in sec. 2.1.1 for topology (c). Disregarding logarithmic factors the 
integral scales as t'"^ with A = '^^{ki + hi + li) in the infrared limit. Thus, we have 

^""'Xn ~ Xj.Wj, ^[x\i+hi+h^ni+m,+pi-ki-hi-U^^ (2.113) 
i 

which shows that these contributions can always be neglected. 
Case (IIC) 

Using a similar notation of case IIB, we have that Xn has the following form 

i 

[ [ V(^^Hp,q,-p-q,0,...)y(^^)(p,q,-p-q,0,...) 

./ p «/ q 

{t + ir(p))-i-S"^(t + X(q))-i-^"^^(i + K{p + ci))-^-^P' 

n{['^^(P)]"'"''['^^(q)r'~''1'^i(P + q)P"'0- (2-114) 

i 

Here n — ji + j2 — 6 + ^^(^ — 2)(n£ + + pi) as before. The analysis follows from the 
analogous one presented in sec. 2.1.1 for topology (d). Assuming without loss of generality 
that — J2hi < should consider three cases: (a) ki = hi = k = 0; (b) 

ki = hi = 0, > 0; (c) Ehi>0 and ^Z, > 0. 
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In case (c) the integral behaves as t ^, with with A = '^^{ki + hi + k) — 1 in the infrared 
hmit. Thus, 

^""'Xn ~ wj.wj, (2.115) 

i 

which can always be neglected. In case (a) the integral behaves as logi, while in case (b) it 
behaves as t~^^\ Therefore, 

f'-^Xn - -Wj,wj,Y[w^'^'^'^P'-^'x'^ - w-^'\ (2.116) 

t t 

SO that we have to discuss some cases. First let us discuss the anomalous diagrams that can 
be cancelled by one-loop insertions. If di — Pi — k, they are the following (with li > 1): 
(A) di = 5j,3, m ^ rrii ^ 0, ji = = 3 {f'^Xn ~ uq'^/'^); (B) di = 2(5^,3, m ^ rrii ^ 0, 
ji = J2 = 3 (r-^Xn ~ (C) di = Si,4, ni^m, = 0, ji = j2 = 3 (r-^Xn ~ (D) 

di = Si^3, ni = mi = 0, ji = 3, j2 = 4 (r~V„ ~ ^6°); (E) di = 6i^3, Ui = Si^s, rui = 0, ji = 3, 
J2 = 3 {t"'~^Xn ^ uq^)- Diagrams (A), (B) and (C) are regularized by insertion on respectively 
one-loop diagrams of D2 counterterm of Uc2 (2.80). Indeed, including combinatorial factors 
(let us consider A): 



^n-i !fL^TT„'i / ^3ip;v --ip,q, -p-q; ysypjysy^) ^9117^ 



t 



rl^3(p)\/(3)(p, q, -p - q)2 V^{p)VM) 



A(p)E«.+2A(p + q)A(q) A(p)^^^+2ii'(q)2 



If we write T^^^^(p, q, — p — q) =: Kslp) + ^^^(q) + V3(p, q) in the first term in the r.h.s. of 
the previous expression (2.117), we find that only V3(q) terms are relevant [indeed for the 
other contributions the integral scale as with A = /i — 1 as in the case (c)] but these 
are regularized by the second term on the r.h.s. of (2.117). The same happens for (B) and 
(C). Similar considerations follow for (D) and (E) using respectively T2 and Ti of Uc3 (2.83). 

Now we have to compute exactly the previous two loop contributions to xi and X2- We 
define: 

v.o)(p, -p - 0. . . .) ^ v;-(p) + v-(qHv;-(p + q) ^ ^^^^^ ^, ^^^^^ 

We notice in particular that = for p = 0, q = or p -|- q = 0. One-point function is 
given by the following two diagrams: 



Xi 



I 



At A(p)2A(q)A(p + q) 

, / V^'Kp, q, -P - q, 0)V^'Kp, q, -p - q) ...... 

^ 6t A(p)A(q)A(p + q) ^ " 



then in the CCL: 



3/4 / , , R^^^ + 



^ 21;IIC _ ^6 .(1) "21;IIC , ^ , ^9 190^ 

XI - ^(^•^2i;iic-^log-j (2.120) 

1 f (.^M' , ^ V3{p?\ , 1 /• vMVsip) ...... 
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The explicit expression for ^li jic quite long and we not report here. 



Two points function computation give us: 



tX2 



t 



I V^^'Hp,q,-p-q 
4A(p)2A(q)A(p + q). 



1 



V:3(p)^3(q) 



+ 



^3(P) 



2A(p)2A(q)2A(p + q) 2A(p)3A(q)2A(p + q) 
V^'Kv, q, -P - ^)V^'\V, q, -P - q, ^ziv) 



A(p)2A(q)A(p + q) 



+ 



l\/W(p,q,-p-q,0)2 , lV^(3)(p^q^_p_q)y(5)(p^q^_p_q^0^0^ 



6 A(q) A(p + q) A(p) 
and in the CCL: 



A(p)A(q)A(p + q) 



(2.122) 



«(2) 
*-*21;IIC 



t 



•^2I;IIC 



»(2) 

^21;IIC 1 t 

■ — log — 

47r ^32 



^3(P)V^2^. ^3(P) 



2 K{^Y 

V4(P 



'^(P) 



1 1 



Kip? 



(2.123) 

(2.124) 
V^3(p)/5 \/4(p)V^3(p) 

Vsipmp) (2.125) 



also in this case we have reported explicitly only the log term of the previous expression. 



2.2.3 Correction to uic and U2c at two loop 

In order to complete the two loop analysis we have to compute the contributions to xi 
to X2 of one loop diagrams with the insertion of counterterms Ud- Then corrections to Ud 
and Uc2 are needed in order to delete non-scaling or non-universal terms. In the analysis of 
lA we have just used the logarithmic part of Ud {Ui and Di) that will be neglected in this 
section. Using (2.75) (2.80) (2.83) (2.84) we have: 

^ _^ [ I + [ y^(P) .2 126) 



and in the MCCL: 




^ r t \ 

^log— (2.127) 
47r ^ 32 ' ^ ' 



3 



^3(P)^4(P) VM 

T ^ / \ o ~r* 1 / \ 



. K(p)3 K{py ^ k{pj 

>{1) 1 / 0^3(P)^ ^ ^5(P) ^ ^ ^3(P)^ ^ V,{p) 
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For X2 we have: 



-i^3)^^6^' / ^4(P) 



+ (^^±|>^/^ (2.130) 



and in the MCCL hmit: 

" ' ^ Me^/^ ip K(p)2 p2 StTMs^/^ M6V2 ypi^(p)3 p2 

3_{D2±Ds)a^ D2 + D3 f Va{p) _ 04 _ a^D^ + D^ 
^6 ^ „, ^, / n I n \ ^, 2 



K>(2) _ Ei=2<3i , ^c3a3 (£'2 + Aj)«3 , (£'2 + ^3)a4 ^ „ 

^^1'^ " 2^6 «6=^/^ u,'/' 2 ^^-^"^ 

where we have used (ref. to Ud) At the end the only log correction comes from lA, while 



ul = -^6 (4?;iB +4?;nc + 4?;c) " 1^ 32x6 i^K, - —\ogUej (2.134) 

ui = -^^6^/' (4;?iB + 4;?iic + 4;?c) (2-135) 

where we have used the fact that: 

= ^gk + ^Sic + ^S?c (2.136) 

= sS;1b + ^2;lic + M;?c (2-137) 

The previous equalities are expected, indeed log terms come from the factorization of two 
loops integrals, so that the one-loop result is recovered. 



2.2.4 Higher powers of the fields 

Wc want to show that higher power of the fields in the starting interaction (2.63) docs not 
affect the MCCL: indeed we want to give evidence of the fact that including higher power 
vertices the scaling equations (2.67) and the critical parameters Ud do not change. In order 
to do that we have to show that every diagram, in which a more than six leg vertex appear, 
is suppressed in the MCCL. Let us consider an extra term in the Hamiltonian (2.63): 

m = Y,u,^^-'y^ j j _^5(p, + ---p,)b, + l^(^)(Pi,---P,)] (2.138) 

At tree-level we have the following additional contributions Ax„ to Xn- 

e-'^Xn ~ ~ u,^^-'^" (2.139) 



2.2 Multicritical crossover limit 
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so that the contributions are null in the MCCL. 

The generalization to one and two loops diagrams is direct. Suppose to have a diagram 
Qn with n-external legs and at least a vertex V^^ with j > 6. ^*_j_,_g is the diagram obtained 
replacing V^^^ in with V^^^ in such a way that the topology remains unchanged. In the 
previous section we have shown that the general behaviour for a diagram in which a six- legs 
vertex is present Q^_j_^.Q is: 

^:_,+6 - t'-^''-'^'\A + B\ogUe) (2.140) 

where A and B are constants. Otherwise the loop-structure of Qn is exactly the same as the 
loop structure of Q*_j^Q; it follows that 

f'-^Qn ~ Ue^-'^y^A + SlogMg) (2.141) 

that shows that graphs with higher order vertex are suppressed in the MCCL. 

In this section we have obtained similar results of sec. 2.1 for Multicritical interactions eq. 
(2.63). In particular several critical parameters can be defined (sec. 2.2.3) so that taking the 
multicritical crossover limit (2.53) the connected n-point correlation functions Xn behave in a 
universal way (2.52). The universality property is not affected by lattice details [^'^■'^(p, ■ ■ ■) 
in eq. (2.63)] or by higher order vertex sec. 2.2.4, but can be defined using only field 
theory. 



^^This can be done cutting j — 6 external lines of V^^^ in Q. We note that this operation is always possible 
up to two loops. 



Chapter 3 



Large- A/^ expansion of 0{N) models: 
the critical zero mode 

In sec. 1.5 we have shown as the standard 1/N expansion fails for certain 0{N) models 
at the finite temperature phase transition. Indeed in the 1/N expansion severe infra-red 
divergences appear that cannot be resummed including corrections to the leading {N = oo) 
order. This means that one -in principle- needs to consider all the orders. However the 
failure of the expansion is suggested because, using a symmetry argument [82], one expects 
an Ising behaviour (if Af — 1) for every N finite. However, we have just pointed out, as in 
principle there is also the problem of the expansion: in [11] the study of 0{N) models 
in one dimension points out that the N = oo limit is affected by unphysical phase transitions 
that strictly disappear when N is taken finite.^ 

This chapter is organised in the following way. In sec. (3.1) we discuss the effective 
interaction of the models with a critical mode (i.e. J\f = 1). We will find an interaction 
similar to what introduced in the past chapter 2, so that introducing proper scaling fields, 
we will be able to describe the crossover between the Ising criticahty with the Mean Field 
one. More interesting, we will present also some numerical results (for the correction to the 
iV = cxo critical temperature) for a pair of model considered in literature [82] [81]. For the 
mixed 0(A^)-RP^~^ model [81], there are no numerical evidences of the existence of a phase 
transition for N — 3; however we predict the possibility that for N < Nc {Nc ~ 100, see sec. 
3.5), the critical point disappears. On the other hand, for the model simulated in [82], there 
are strong evidence that for = 3 the system undergoes a phase transition. Our numerical 
prediction for the non universal critical constant Pc seems to have the same magnitude order 
of that measured in [82] , however the accordance is not very good due to the big difference 
between Pc{oo) (fsi 5, see [55]) and]9c(3) (~ 20, see [82]). More stringent check will come from 
other J\f — 1 models (see chapter 5), for which simulations with several N are available [50]. 
We stress that this could be an important test for the scheme presented in this work, and 
due to its universal setting, it could be applied to all the model with Mean Field crossover. 



phase transition (i.e. a singularity in the partition function) can appear only if infinite degrees of 
freedom appear in the system (Yang-Lee Theorem [21] [22]). This can be realised both with a (physical) 
thermodynamic limit 1^ ^ oo or taking the (unphysical) N ^ oo limit. 
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Large-N expansion of 0{N) models: the critical zero mode 



3.1 Effective Hamiltonian for the zero mode 



We have discussed in the first Chapter that, for the class of phase transitions we are interested 
in, the propagator P has a zero mode at the critical point. In order to understand the role 
of fluctuations, a careful treatment of the zero mode is required. For this purpose, we 
are now going to integrate out the massive modes, obtaining an effective Hamiltonian for 
the critical field 0. More precisely, we define 



-nm 



(3.1) 



The effective Hamiltonian 7ieff[0] has the following expansion 



N 



(3.2) 



n=3 ■ ''Pi "'P" i 

where vertices and propagator also depend on and have an expansion of the form 



m=0 



m=0 



1 ~ 



m=0 



(3.3) 



We report here the explicit expressions that we shall need in the following: 



^^^IJ E^i»f(0'P'-P)^«f(P)' 



ab 

Po\p)^Pu\p), 

1 



(3.4) 
(3.5) 



Y.vSUp,-p,Ci,-q)Pab{ci) 



ab 



E ^naiP: -P, 0) P„(.(0)V;S(0, q, -q) Pcd(q) - 

abed 

E (p, -p - ^WSl (p, -p - q)i^ac(q)A<i(p + q) 



abed 

T>(3)/ A _ T>(3) 



CHP,q>r) = V;\^/(p,q,r), 
^^0^^ (P, q, r, s) = Vlfl-^ (p, q, r, s) - 

E ^i?a(P' q' ~P ~ q)^n6V> s, -r - s)P„6(p + q) + two permutations. 



(3.6) 
(3.7) 

(3.8) 



ab 



Vq^^ (p, q, r, s, t) = Vl°l^ (p, q, r, s, t) 



(5) 
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Yl ^i2(P' ~P ~ q)^m6(r, s, t, -r - s - t)Pabip + q) + 9 permutations 



ab 



+ 



J2 Vl?a{p, q, -P - q)V;2 (r, p + q, s + t)V^/2(s, t, -s - t) x 



.abed 



X-Pa&(p + q)-Pcd(s + t) + 14 permutationsj , 



(3.9) 



where a, b, c, d run from 1 to 4 over the massive modes. 

The vertices introduced above are not independent, but near the critical point can be 
related to the scaling fields introduced to parametrize the gap equation eq. (1.10). The 
identities presented in sec. 3.1.1 allow us to derive several relations among the effective 
vertices at the critical point. We have for p — Pc and ml — 



V^J'^ (0,0,0) 



drriQ 
= 0, 



d 



^0^'^ (0,0,0) 



1 2 



^V-r(0. 0.0.0). 



d{ml) 



(3.10) 
(3.11) 
(3.12) 



Relations (3.10), (3.11) and (3.12) are the first step to show the effective interaction for the 
zero mode is a weakly coupled interaction studied in the previous chapter. On the other 
hand (3.10) clearly confirm that the standard 1/N expansion fails close to the critical point. 
Indeed, outside the critical point, Pq^{p) is nonsingular and for large N it is enough to 
expand the interaction Hamiltonian in powers of 1/N. On the other hand, this not possible 
at the critical point. Since also the three-leg vertex vanishes at zero momentum in this 
case, cf. eq. (3.11), the zero- momentum leading term is the quartic one. Since the coupling 
constant is proportional to 1/A^, the model effectively corresponds to a weakly coupled (f)'^ 
theory. In order to have a stable (j)'^ theory, we must also have V"q*^''''(0, 0, 0, 0) > 0. For a 
generic solution of the gap equations satisfying eq. (1.14), this is not a priori guaranteed. 
Note, however, that if V'o^^^(0, 0, 0, 0) < 0, then, for p — p^, we have -Po~^(0) ~ a(ml — ml^y 



■^OeJ 

with a < 0, as a consequence of eq. (3.12): the propagator has a negative mass for N = oo. 
We believe — but we have not been able to prove — that such a phenomenon signals the fact 
that the solution we arc considering is not the relevant one. We expect the existence of 
another solution of the gap equation (1.10) with a lower free energy. 

In the previous chapter we have seen that the weakly coupled 0^ theory shows an in- 
teresting crossover limit. If one neglects fluctuations it corresponds to tune p and ml so 
that H and P~^{0) go to zero as oo, in such a way that HN and P~^{0)N remain 

constant. In this limit, Ising behavior is observed when the two scaling variables go to zero, 
while mean-fleld behavior is observed in the opposite case. Fluctuations change the simple 
scaling forms reported above and one must consider two additive renormalization constants 
{hciu) and rdu), in chap. 2). 

However in order to strictly apply the results of the previous chapter (including also the 
computation of the renormalization constants with the lattice regularization) we need to 
impose the condition V^^\0,0,0) = for every values of the parameters, while in the case 
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under investigation this happens only near the critical point (3.11).^ In order to do that we 
change again the definition of the field so that the effective zero-momentum three-leg vertex 
vanishes for all p and ml in the limit N ^ oo. For this purpose we now define a new field 

axip)=Hp) + kS{p), (3.13) 

where a and k arc functions of p and m^ to be fixed. The function k is fixed by requiring that 
the large- zero-momentum three-leg vertex vanishes. Apparently, all V^"^ contribute in this 
calculation. However, because of eq. (3.11), t^Q*^^'*(0, 0, 0) vanishes at the critical point. As 
we already mentioned the interesting limit corresponds to considering = — m^^ — > 
and Ap = p — ^ together with N — > oo. We will show in sec. 3.3 that this limit should 
be taken keeping fixed A^iV^'''^ and ApN, so that yQ''^'*(0, 0, 0) is effectively of order N~^/^. 
Therefore, the equation defining k can be written in a compact form as 

^ + V =0 (3 14) 

n>4 

where a„ ~ a„o + o,ni/N is the contribution of the n-point vertex at zero momentum. 
Eq. (3.14) can be rewritten as 

n=l ^ ^ 

which shows that k has an expansion of the form 

k = koN^'^il + kiN-^'^ + C>(7V-2/3)]. (3.16) 

The leading constant /cq depends only on the three- and four-leg vertex, the constant ki 
depends also on the five-leg vertex, and so on. The explicit calculation gives 

k = VivjfteL + 



T/o^'^(0, 0,0,0) 



N vi'\Q, 0, Qfvi'\Q, 0, 0, 0, 0) ^ o(Ar-V2). 



2 [Vo(^)(0,0,0,0)]3 

By performing this rescaling, the /c-leg 0-vertex that scales with N as N^'^/"^ (except for 
k <?>) gives a contribution to the m-leg x vertex (of course m <k) ol order ]Sf'^-'>n/'^H-m-k)lz ^ 
which is therefore always subleading. Taking this result into account we obtain 



T^eff = i^X(0) + ^^x(p)x(-p)P-^(p) + 

/ ^^^Hp>q>-p-q)x(p)x(q)x(-p-q) 

J D,a 



3!VAr 



■T^ / V^^\p, q, r, -p - q - r)x(p)x(q)x(r)x(-P - q - r) + . . .(3.18) 

Jp,q,r 



^However the fact that y^'^^ (0,0,0) goes to zero at the critical point is a fundamental topic in order to 
apply the following considerations. In particular this implies that k, defined in (3.13), can be computed 
perturbatively (3.17) near the critical point. 
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where 



H — , — — ak 



+ ^v^J3)^o,o,o)-^v^r(o,o,o,o) 



P-\P) = a%-\p) + ^Pr^(p) - ^V<'\0, p, -p) + ^V,^'\0, 0, p, -p) 



2 1.3 



a k 



vi'\0,0,0,p,-p)+O{N-'/'), 



6Ar3/2 



l-(^)(p,q,r,s) = «Vo(')(p,q,r,s) + 0(7V-V3), (3.19) 
where a — aQ-\- aiN~'^/^ + 0{N~^l^) is fixed by requiring that 

P-^(P) - P-^(O) = K{p) ^ p^ (3.20) 

for p — > 0. 

3.1.1 Identities among effective vertices 

In this section wc wish to derive a general set of identities for the zero-momentum vertices 
in order to justify the relations introduced above (3.10,3.11,3.12). For this purpose we first 
rewrite the propagator and the vertices by considering the constants a, 7, and r appearing 
in eq. (1.6) as independent variables without assuming the saddle-point relations. We define 
the integrals 

^'-•■"'"'"'^A (7-«E„cos,,.)- - 

It is easy to verify that, if a and 7 are replaced by their saddle-point values, a = 2W'{f) 
and 7 = (4 + ml)W'{f), then 

In terms of a, 7, and r the propagator is simply obtained by using eq. (1.57) and replacing 
.,(0, mQ)/[W^'(r)]^ with Aij^2{(^il)- Let us now consider a generic n-leg vertex Vj^'^ at 
zero momentum. It is easy to verify that the only nonvanishing terms with = 4 or = 5 
for some i are (in this section we do not explicitly write the momentum dependence since in 
all cases we are referring to zero-momentum quantities) 

vt}A = vt, = -PW^''\T). (3.23) 



If aU Ai < 3, eq. (1.61) gives 

~ 2 



vit.Au = k-ir^''^''^\n - l)!A,,3,„(a,7), (3.24) 
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where ^2 (resp. i^) is the number of indices equal to 2 (resp. to 3). These expressions allow 
us to obtain simple recursion relations for the vertices. 

T/("+l) _ At/W 

T/("+l) I T/(n+l) _ _^T/H 
^2,Ai,...,A„ ^3,Ai,...,An — Ai,...,An^ 

..(n+1) _ d 

r{n+l) _ d (n) 
'5,5,.. .,5 — >^5,...,5> 



nt',i = 5> (3-25) 



where a, 7, r, and /9 are considered as independent variables. These relations also apply to 

(2,) 1 

the case n = 2, once we identify V^^ = P^b- 

Let us consider the projector on the zero mode va, cf. eq. (1.70). Keeping into account 
that the symmetry of the matrix at zero momentum implies V2 — and V4 — V5, we obtain 

Close to the critical point we can write, cf. eq. (1.66), 

VA = ZA + 0{si) (3.27) 

with 

= K (2'^,1,1, , y (3.28) 

where K is such to have YIib = 1- Now, let us note that at the saddle point we have (we 
now think of a, r, and 7 as functions of mg and p) 

=2iy"#I, = ^#I,i2 + O(.0, (3.29) 



9mo dml K dml 

where we have used eqs. (1.65) and (1.59). Thus, if we define 



we can rewrite 



T/(n+i) _ n f d da d dr d\ („) 

2^vbVb^a,,...,a^ - + + + 



C^ytLA. + 0{s,l (3.31) 



''0 
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where in the last term a, 7, r, and j3 take their saddle-point values in terms of tUq and p. 
To go further we compute the derivative of za with respect to mg. Since 

Y^PXk^B^Ois,), (3.32) 

B 

we have 

E^I^H-Eff^^ + OW. (3.33) 

where S2 = d'^f3/d{ml)'^ also vanishes at the critical point. 

Since ^^5^ = 1, dzB/dml belongs to the subspace orthogonal to za- Thus, if P^^ is 
the inverse of in the massive subspace, we obtain 

= - E P^B^M^c + 0{s,) = -i 5: PLviin^cZu + 0(s,). (3.34) 



Finally, let us compute 

7 > 7'. A 7. n 

' dm; 



= C Y: ^M...ZA.^Vtl,A. + 0{s,) 







ab 



where the indices a and b run from 1 to 4 over the massive subspace. Since S2 is of order 
p — Pc and itLq —ttLq^ this relation implies an identity only at the critical point. However, for 
n = 2 we can obtain another relation. Since V^^ = P^^ and Pab^b = 0(si) we do not 
need eq. (3.34) and obtain directly 

Vlii = C'^Al' + 0(.i). (3.36) 

The presence of corrections of order si gives rise to two critical-point identities. Since 
P{^ ~ si we obtain 

Viii = 0, 

- — 2^111 — ^"^7 — 2^2-^11 ' (3.37) 

arriQ o{mQy 

where all quantities are computed at the critical point. 



3.2 Critical crossover limit in the large-N case 

The main result of the previous section 3.1 is the identification of the zero mode Hamiltonian 
with the weakly coupled theory introduced in the previous chapter (2.2). One can work in 
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order to apply the same results of chapter 2 by using the following identifications 

chapter 2 chapter 3 

u = ^F^'^(0, 0,0,0) (3.38) 

^^"Hpi,---Pn) = 4) ^'"''"''i. (3-39) 

^(0,0, 0,0) V 

K{p) = P(p)-P(0) (3.40) 
r = P(0) (3.41) 
H = H (3.42) 

so that the normalization conditions assumed in sec. (2.1) are satisfied. In particular one can 
obtain the mass- and the magnetic-term regularisation (r,. and he) in the large limit, using 
the function hc{u) and rc{u) obtained in the previous chapter'^ for the lattice regularisation 
using u defined in eq. (3.38). 

However in our large- A?" case we have to take into the mind that the propagator and the 
vertices of the theory are function of = ttlq — m^^ and of Ap = p — p^., and moreover all 
quantities, beside r and H, depend on these two variables. We have seen in the previous 
chapter that the Critical Crossover Limit (CCL) is defined having fixed two scaling variables 
Xt and Xh 

P(0) - r,{u) 
1^^(0,0,0,0) 

,r H — hJu) /„ , ,S 

= (3.44) 

0,0,0) 

Expanding the prcvioTis couple of equations (3.43,3.44) near the critical point one can obtain 
Am and Ap as function of Xf x^ and A^. As we will show in the next section A^ and Ap 
scale respectively as 1/A^^/^ and 1/A^, so that we can assume an additional dependence on 
^^1/3 for the vertices, the propagator and the counterterms of the theory: i.e. we consider 
X(p; u^/^) and ^(^^^(pi, . . . , Pit; u^/^). Note that, by definition, X(0; v}'^) = for all values 
of u. Moreover, we assume, as in the case of interest, that K^^^(0, 0, 0; w^^^) = for all 
values of u. Following the calculation presented in sec. 2.1.1, it is easy to realize that the 
dependence of the vertices on u is irrelevant except in hdu). In this case, eq. (2.15) becomes 

t K{u) 1 f Vs{p;uy^) , 

— Xl = ;= / TT7 TTTx hC(VM)- (3.45) 

Because of the prcfactor 1 / ^/u we must take here into account the first correction pro- 
portional to u^/'\ Thus, if V3{p;u^/^) ^ T^3,o(p) + u^^^V3,i{p) and K{p;u^/^) ^ Ko{p) + 
u^^^Ki{p), we obtain 



t hJu) 

-Xi = 

u u 



Ko{p)+t {Ko{p)+t)^ 



+ 0(7/1/^). (3.46) 



Vc(u) and hc{u) are given by eq. (2.47) and by eq. (2.17). 



3.3 Crossover between mean-Geld and Ising behavior 
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Now, t can be set to zero without generating infrared divergences, neglecting corrections of 
order u~^^'^tlnt ~ u^^'^lnu. It follows 



hc{u) 



VsM . l/3 ^3,l(p)i^0(p)-^3,0(p)i^l(p) 



(3.47) 



Eq. (3.47) represents the only equation in which the explicit dependence of the vertices on 
Am should be considered (Ap is proportional to u and thus it can always be set to zero). In 
all other cases, we can simply set = A^ = 0. 



3.3 Crossover between mean- field and Ising behavior 

In this section we wish to apply the above-reported results for the critical crossover limit to 
the Hamiltonian (3.18) of the zero mode. With the identification (3.38-3.42) the effective 
Hamiltonian (3.2) corresponds to the Hamiltonian discussed in sec. 2.1. 

We begin by reporting the additive renormalization constants. The mass renormalization 
constant is given by, cf. eqs. (2.47), 

u 



r(N) = — In f - — (3 + 87rL'2) 



+- 



f [l^(3)(0,p,-p)]^ yW(0,0,p,-p) 1 

2 I Vi^) (0, 0, 0, 0)ir(p)2 (0, 0, 0, 0)K{p) p2 



(3.48) 



where D2 is a nonperturbative constant defined in Ref. [32] (numerically D2 ~ —0.052). As 
discussed in sec. 3.2, we can compute all quantities appearing in eq. (3.48) at the critical 
point and keep only the leading terms for N ^ 00. We thus obtain 



TciN) 



u 

8^ 



In 



3m 



2567r 



u 



— (3 + 87rL'2) 

OTT 



(3.49) 



a 



2N 



[Vi'\0, p, -p)Po(p)]' - C^(0, 0, p, -p)Po(p) + ^K''(0, 0, 0, 0) 



a ~ 



(4), 



p. 



where all quantities are computed at the critical point. 

Analogously, we should introduce a counterterm for the magnetic field: 



h(m--l r f ^^^Ho,P, -P) 

nc[iV 2^''yp[t/(4)(0,0,0,0)]V2x( 



p)' 



(3.50) 



As discussed in sec. 3.2, such a quantity should not be simply computed at the critical 
point, but one should also take into account the additional corrections of order N~^^^. Since 
k ~ 0(ArV6) ^j^g critical crossover limit, we have 



hc{N) = 



2VN 

a 



V^,(')(0, 0,0,0) 



v;,(')(o,o,p, -p) 



^0^'^ (0,0,0) 



2\/iVVo(^)(0, 0,0,0) 



[vi'\o,p,-p)Po{p)r. 



(3.51) 
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The second integral should be computed at the critical point, while the first one and the 
prefactor of the second one should be expanded around the critical point. Since, as we shall 
show below, = ml — ttIq^ ~ N~^^^ and p — Pc ^ N~^, it is enough to compute the first 
correction in A^. In practice, we find that the renormalization terms have the form 



rc{N) 
hc{N) 



N 



1 



[rcolnN + rci), 

(hco + ^mhcl)- 



(3.52) 



Once rc{N) and hc{N) are computed, we can define the scaling variables ~ t/u and 
~ h/u. Choosing the normalizations appropriately for later convenience, we define 



N 

Xh = -[H-h,{N)\. 

a 



(3.53) 



Now, the critical crossover limit is obtained by tuning p and around the critical point in 
such a way that Xt and Xh are kept constant, i.e. p — >• Pcrit, ""^o ~^ '^o.crit' N ^ oo aX. fixed Xt 
and Xh- Note that here Pcrit and tUq^^-^^^ correspond to the position of the critical point as a 
function of N (thus Pcrit ~^ Pc ^-nd tt^q ^-rit "f^oc as A?" — > oo) and are obtained by requiring 
Xt = Xh = Q (cf. sec. 2.1.4). 

In order to compute the relation between p, and Xt-, Xh we set A^ = to^ — m^^ and 
p — Pc and expand eq. (3.53) in powers of A^ and Ap. In the following we shall show 

N~'^, so that the relevant terms are 



that A^ ~ 



N-^l^ and Ap ~ 



Xh 



N{btoAp + btiApArr, + 6t2A^) + dilnN + d2 



•hoA^ + 6wApA, 



+ 6,,2ApA^ + h, 



7i3 



Al + hh.At) + 



(3.54) 



(3.55) 



The constants are obtained by expanding P~^(0) and H around the critical point and by 
using the expansions (3.52). All quantities are analytic in A^ and Ap and several terms are 
absent because of identities (3.10), (3.11), and (3.12). In particular, a term proportional to 
A^ is absent in the equation for xt- This is a consequence of eq. (3.12). Indeed, we have 



1(0) 



(0) 



2\/o^^^(0, 0,0,0) 



+ 0(A^,iV-i) 



d{mlf ■(/^('^(0, 0,0,0) 



d 



dml 



Al + 0{Al„N- 



= + O(Ai,N-'), 



(3.56) 



where in the last step we have used eq. (3.12). 

We wish now to determine the behavior of A^ and Ap that is fixed by eqs. (3.54) and 
(3.55). We assume 



NO 



(3.57) 



3.3 Crossover between mean-Geld and Ising behavior 
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where 6mo and 6pQ are nonvanishing constants and a and P exponents to be determined. If 
/3 < 1, eq. (3.54) implies 

btoSpoN'-^ + bt25moN'-^'' = o(7V'-^) , (3.58) 

which requires (3 = 3a. Now, consider eq. (3.55). The term N^^'^Af^ is of order iV^/^~" and 
cannot be made to vanish. Therefore, we must have (3 > 1. Considering again eq. (3.55), 
it is easy to see that all terms containing Ap cannot increase as fast as A^^/2_ Therefore, 
cancellation of the term d^oN^^'^ requires a = 1/3. Consideration of eq. (3.54) implies finally 
[3=1. This analysis can be extended to the subleading corrections, obtaining an expansion 
of the form 

A ^mO , ^ml , ^m2 



jYl/3 jY2/3 jY5/6- 



= % + ll^ + 0(iV-/^). (3.59) 



The coefficients are given by 



5, 



mO 



"30 



^5 

+ {bhibto - bhibt2), 



5po = --^(Ot2 + c^ilniV + (i2-Xt). (3.60) 

We are not able to compute 5m2 and 5pi but we can however compute a relation between 
these two quantities. We obtain 

35TO05m2&/i3 ,1 fn 

dpi = h T-^^ — Xh. (3.61) 

Correspondingly, by using eq. (1.21), we can compute the expansion of Uh'- 

UhO ^ Uhl Uh2 „„s 

where w/ji depends on Xt and on ^^/j. We thus define a new scaling field by requiring that 
no term proportional to XtN~'^/^ is present. For this purpose we set 

Uh = Uh + 1^ (p - Pc) (3.63) 

where the coefficient Xmix is given by 
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The new scaling field has an expansion of the form 



UhO Uhl 



+ 



Uh2 

Ar3/2' 



where 



Uhl = 



«03()m0' 



Uh2 



0.03 



-003*^31 + (O04&/i3 



a03&M)^mo] 



-Xh- 



(3.65) 

(3.66) 
(3.67) 

(3.68) 



Interestingly enough, in 11^2 all terms proportional to the unknown quantity 5^2 cancel. At 
this point we can easily compute the expansion of the critical point Pcnti /^crit- It is 
enough to set Xh — Xt — Q va. the previous expansions, obtaining 



Pcrit ^ P<^^ 



+ 0{N-'/^), 



xt=0 



It follows that 



P - Pcrit 



Uh ~ Uh,cnt 



2 btibhi6lo 



+ 1 



3 rfao / btoN 



(3.69) 
(3.70) 

(3.71) 
(3.72) 



where Uh,crit is the value of Uh at the critical point. Therefore, the 1/A^ corrections modify 
the position of the critical point and change the magnetic scaling field which should now be 
identified with Uh = uh ^ w/i.crit (as we already discussed the thermal magnetic field is not 
uniquely defined and we take again p — pcnt)- Eqs. (3.71) and (3.72) also indicate which are 
the correct scaling variables. Ising behavior is observed only if N(j) — p^nt 

) and N^''^{uh - 

Uh,ci:it) are both small; in the opposite case mean-field behavior is observed. Therefore, as 
N increases the width of the critical region decreases, and no Ising behavior is observed at 
N — 00 exactly. 



3.4 Critical behavior of (cr^ • (Tx+pt) 

In this section we wish to compute the large- A?^ behavior of 

E^{(T^-(T^+^). (3.73) 

Such a quantity does not coincide with the energy. However, as far as the critical behav- 
ior is concerned, there should not be any significant difference. In order to perform the 
computation, note that the equations of motion for the field A-^^ give 



{Px„) = l + E. 



(3.74) 



3.4 Critical behavior of {a^ ■ cr^+n) 
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Thus, we have 

E = T-1 + -^ip.,) = T - 1 + ^ ^ U,Bm^B.). (3.75) 
ViV Viv ^ 

Therefore, we need to compute the correlations (0) and {(fa)- For the zero mode we have 
immediately 

{(t^x) = Oi{Xx) - k = —fr^'^ixt, Xh) - k, (3.76) 

Xt 

where fl^^^{xt,Xh) is the crossover function for the magnetization in the Ising model and 
the constant k diverges as N^^^, modulo some obvious normalizations. 

Let us now consider (lPo) and show that such a correlation vanishes as iV — >■ oo in the 
critical crossover limit. Indeed, we can write 

{cpa) = ^^6(0) /[fUO,P,-p)i^cd(p) + (Hii(0,p,-p)-nii(0,0,0))Ai(p)] 

V W Jp 

+-^A6(O)Hii(O,O,O)(0^). (3.77) 



Note that in the last term we have replaced 

Ai(p) = (0^), (3.78) 



a necessary step, since the integral diverges at the critical point and therefore should be 
computed in the effective theory for the zero mode. Now, we have 

{<f>l)^e-2ak{xx)+a'{xl). (3.79) 

In the critical crossover limit, the two expectation values are replaced by the crossover 
functions for the magnetization and the energy and by a regular term. Therefore, for — oo, 
the leading behavior is {(pl) = A;^ ~ N^^^. It follows that (ipa) ~ N"^/^. In conclusion we 
can write 



— [xt,Xh)-k 
Ereg + ^U,,{0)-fr'^{xuXu) + 0(iV-2/3), (3.80) 

V -/V Xt 



where E^^„ is the regular part of E: 



Ereg = r-l-4^[/4i(0)A; 

V W 



= ro-l + [n-koU^im5moN-'/^ + 0{N-^/^), (3.81) 
where we have written r ~ tq + tiA^ and k ~ koAmVN- 
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Equation (3.80) shows that the singular part of E behaves as the magnetization in the 
Ising modeL For Xt ^ I and Xh ^ 1 one observes Ising behavior and thus 

E — Ej-eg ~ \xtf' for Xh — 0, low t phase 

E - E,,g ~ \xh\^/^' for Xh ^ 0, (3.82) 

where (3i = 1/8 and 5/ = 15. On the other hand, in the opposite limit we have 

E — Ereg ~ \xtf'^'' for Xh — 0, \xt\ — > oo, low t phase 

E -E,^gr^ \xh\^/^^^ for \xh\^oo (3.83) 

with (3mf = 1/2 and Smf = 3. Note that the limit \xt\ — > cxo and \xh\ oo should always 
be taken close to the critical limit. Therefore, \xh\ oo means that we should consider 
N ^ oo, Uh ^ Uh,crit, P Pcrit in such a way that N^^^{uh - Uh,crit) oo, i.e. N should 
increase much faster that the rate of approach to the critical point. 



3.5 Numerical results for selected Hamiltonians 

In this section we present some numerical results for some selected Hamiltonians. First, as 
in Ref. [82], we consider 

W{x)^^{^y. (3.84) 
Second, we consider the mixed 0(A^)-RP^~^ model with Hamiltonian [81] 

n = -NPv J](<T, • <7,+^) - ^ J](<7, • o-,+^)2. (3.85) 

This Hamiltonian corresponds to the function 

W{x)=px+^{l-p)x'^, (3.86) 

where we set /3y = (1 +p)/3/2 and /3r = (1 — p)/3/2. This Hamiltonian is ferromagnetic for 
p> -1. Note that for p = -1 we obtain the RP^'^ Hamiltonian [89], [81], [90], [91], [92], 
[93], [94], [84], [95], [85] 

^^-^EK-'^-^-^)'' (3.87) 

that has the additional gauge invariance cr^ ^x^^x, = =tl- Under standard assumptions, 
the large-// analysis should apply also to this last model: the local gauge invariance should 
not play any role.'' 



''The irrelevance of the Z2 symmetry for the large-/? behavior of RP^~^ models have been discussed in 
detail in Ref. [94, 95, 84]. Thus, in spite of the additional local invariance, RP^^^ models are expected to 
be asymptotically free and to be described by the perturbative renormalization group. 
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Table 3.1: Numerical estimates for interaction (3.84) (Wi) and (3.86) {W2). 

In the large- iV limit [55], the first Hamiltonian has a critical point at f3c ~ 1.335 and 
Pc ~ 4.538. By using the numerical results reported in Table 3.1 we can compute the first 
corrections to the critical parameters. We obtain 

/3c ^ 1.335 + ^(36.127 + 2.093 In AT), (3.88) 
Pc ~ 4.538+ (87.92 + 4.80 In AT) . (3.89) 

Note the presence of a In N/N correction due to the nonanalytic nature of the renormalization 
countcrtcrms. The correction terms are quite large, indicating that the large- A?^ results are 
quantitatively predictive only for large values of A^. This is not totally unexpected since [82] 
Pc ~ 20 for A^ = 3, that is quite far from the large- A^ estimate Pc ~ 4.538. If we substitute 
A^ = 3 in eq. (3.89), we obtain pc ~ 35, which shows that the corrections have the correct 
sign and give at least the correct order of magnitude. 

For A^ = 00 the second Hamiltonian has a critical point at [81] /3c ~ 0.918 and Pc ~ 
—0.971. Including the first correction we obtain 

f3c = 0.9182 --^(11.062 - 1.437 In A^), (3.90) 

Pc = -0.9707 +-^(2.905 - 1.265 In A^). (3.91) 

In this case the corrections are smaller. However, for A^ < 99.4 they predict Pc < — 1, 
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although only slightly. This is of course not possible, since for p < —1 the system is no 
longer ferromagnetic. Thus, we expect the transition to disappear for some N = N^, with 
Nc ~ 100 (this is of course a very rough estimate). Since Pc ~ —1, we also predict RP^~^ 
models to show a very weak first-order transition for large N . 

3.6 Conclusions 

In this chapter we have explicitly presented our technique to solve the problem of Infra Red 
Divergences in large N expansion, giving a detailed description of a class of one-parameter 
0{N) models that present this kind of problem. This models show a hne of first-order 
finite-/? transitions. We focus on the endpoint of the first-order transition line where energy- 
energy correlations become critical, while the spin-spin correlation length remains finite, in 
agreement with Mermin- Wagner theorem [62]. In sec. 1.5 wc showed that, at the critical 
point, the standard 1/N expansion breaks down, since the inverse propagator of the auxiliary 
fields has a zero eigenvalue. A careful treatment shows that the zero mode, i.e. the field 
associated with the vanishing eigenvalue, has an effective Hamiltonian that corresponds to 
a weakly coupled one-component 0^ theory. Thus, the phase transition belongs to the Ising 
universality class for any A^, in agreement with the argument of Rcf. [82]. In Rcf. [55] it was 
shown that for N = oo the transition has mean-field exponents. We reconcile here these two 
results. If Ut and Uh are the linear thermal and magnetic scaling fields, in the critical limit 
a generic long-distance quantity O has a behavior of the form 

(e'>sing ~ uUoiutN, UhN^/^), (3.92) 

where f{x, y) is a crossover function. Only if UfN <^ 1 and UhN^^"^ ^ 1 does one observe 
Ising behavior. In the opposite limit one observes mean-field criticality. Therefore, the width 
of the Ising critical region goes to zero as — > oo and, even if the transition is an Ising 
one for any A^, only mean- field behavior is observed for A^ = oo. The behavior observed 
at the critical point for A^ — > oo resembles very closely what is observed in medium-range 
models [36], [37], [34], [35], with A^ playing the role of the interaction range. Our analysis 
fully confirms the conclusions of Ref. [82]. 

From a more quantitative point of view, we give explicit expressions for the critical values 
Pc and Pc and for the nonlinear scaling fields to order 1/A". Numerical results are given for 
the Hamiltonian introduced in Ref. [82] and for mixed 0{N)-RP^^^ models [81]. We plan 
to compute corrections to the critical parameters for other model for which more available 
data are present (for instance Yukawa models, chap. 5 and ref. [50]). This would provide a 
more evident check for the scheme proposed in this work. 

The main result of this paper, i.e. the fact that the width of the Ising critical region 
goes to zero as A^ — > oo, holds in any d < 4. However, in generic dimension d the natural 
scaling variables are (p — Pcrit)A^^^^^~'^^ and {uh — ttfe,crit)A^^^'-^~'^\^ while the scaling equation 
eq. (2.3) is now replaced by 

^ ^d/(4-d)^-n(d+2)/4jsymm^^^-(4-d)/2^^^-2(4-d)/(2+(i) 

'""Thus, for d > 2 the Ising critical region shrinks faster as N increases: in three dimensions as iV~^ and 
N^'^ in the thermal and magnetic directions respectively. 



(3.93) 



3.6 Conclusions 
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The claim simply follows from the fact that for any d the zero mode has a Hamiltonian that 
corresponds to a weakly coupled theory. Thus, for d < 4, the phase transition is always 
Ising-like. One only needs to take into account the different definition of the scahng fields 
(3.93). Also the expressions for the renormahzation constants Tc and he should be changed: 
for instance, in three dimensions we also expect contributions from two-loop graphs, as it 
happens in medium-range models [32]. More important, following the discussion of sec. 
(2.1.5) in order to get the previous scaling equation (3.93) one need to correct the definition 
of the zero mode (2.57) to impose that the three-leg vertex renormalized must be equal to 
zero unless contributions that are irrelevant in the Critical Crossover Limit. 



Chapter 4 



Large- expansion of 0{N) models: 
the multicritical zero mode 

In this chapter we want to generahze the considerations of the previous chapter to the case 
in which a multicritical zero mode appear J\f > 1} In this case the study of the effective 
interaction for the zero mode would require a more involved algebra which gives several non- 
trivial cancellations in the expansion of the scaling fields near the multicritical point. This 
cancellations generalize what found in chapter 3 [where for instance the term disappears 
in the expansion of Xt eq. (3.54)] and appear as soon as one consider the translated theory 
(3.13). The main effects of this cancellations will be the identification of the leading order 
of the effective vertices for the translated theory with the scahng fields for = oo -i.e. 
without 1/N corrections- (1.48) introduced in order to parametrize the gap equation near 
the critical point (1.47).^ Following the notations of the previous chapter, if we define 
V^*^-')(p) = V^^\0, ■ ■ ■ 0, p, — p) the vertices of the zero mode, the first result will be to show 
that near the multicritical point 

~ uj_i, A^uj, Ajuj+i---Arn-^-^+^ (4.1) 

for j = 0, ■ ■ - A/". In the second line of (4.1) we have written all the contributions that scale 
in the same way in the N = oo critical limit defined by (1.48) 

Ar+2-j 1 

In the next section we will show that introducing a proper translation for the zero mode 
[similar to what introduced in (3.13)], one is able to cancel in eq. (4.1) all the contribution 

that factorize A^ so that V (0) ~ Uj-i-^ This result is very important in order to include 

^Referring to the notation introduced in the first chapter, we remember that only the J\f odd case is 
investigated, the other case being unstable for A'' = oo. 

^This identification will be exact apart some contribution that are irrelevant in the scaling limit taken in 
the past chapters (also called Critical Crossover Limit ~CCL-). 

^As in the previous chapter V^^^ refers to the vertices for the translated theory. 



73 



74 



Large-N expansion of 0{N) models: the multicritical zero mode 



1/N fluctuations to the N = oo theory of chapter 1. Indeed in chapter 2 we have obtained 
that the scaling fields describing the crossover to the multicritical point are given by 

u,{N) = V^'^'\0)-Uci (4.2) 

where regularise the theory and cancel the nonuniversal details. If the Am contributions 
of the second line of eq. (4.1) do not cancel, inserting (4.1) in (4.2) one would obtain an 
expression similar to the following 

u,{N) = «,(oo)+ J] c(%(oo)+o(-) (4.3) 
j=i+i 

If cf^ ^ for every i and j > i, then the scaling fields would be discontinuous for = oo, 
so that our mechanism to explain crossover between the two different critical regime would 
fail. We remember that there is some freedom in the definition of the scaling fields Ui (1.48). 
Indeed also 

i-l 

u[ — Ui + UjUj (4-4) 

j=0 

satisfy the scaling relations of the N — oo theory; however eq. (4.4) can never match with 
eq. (4.3). We stress how in principle the require to cancel J\f{J\f — l)/2 parameters [cj''^ 
in eq. (4.3)] tuning a single parameter (the translation constant of the zero mode) looks 
an impossible work! The answer to this paradox relies on the fact that the vertices of the 
effective theory arc not independent but arc strictly related each other and at the end to the 
gap equation. In the first part of this chapter we want to elucidate this relations by using 
similar techniques of sec. 3.1.1, while in the last part of this chapter we will describe the 
crossover behaviour to mean field in a similar way of what done for the J\f — 1 case, chapter 
2. In particular we will describe (in a general way), how to obtain the scaling field in the 
large limit and how to include 1/N corrections to their N — oo definitions computing 
explicitly the J\f = 3 case. 



4.1 Effective Hamiltonian for the zero mode 
4.1.1 Basic definitions 

As in the previous chapter the starting point is the definition of the effective interaction for 
the zero mode through the integration of the massive mode of the theory. However in this 
case, being interested in general relations, we need to take a more general point of view. If 
va{p] Pi) is the normalized eigenvector that corresponds to the zero eigenvalue for p = 
at the MCP we have that we define the critical field (rescaled) 

<^(P) = -J=X]^^A(p;mg,Pi)*A, (4.5) 
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while the effective Hamiltonian becomes 



where 



Zlf] = / n H-i-"- vVN)e-^, (4.7) 



where the functional 5 function guarantees the identification (4.5) and \E' ■ f = "^aVa- 
In this expression we assume that VA{p',iTLQ,pi) is normalized, i.e., X^a"^! — 1- Note the 
additional factor \/N in the definition oi cp. It guarantees that, in the large- limit, the 
effective Hamiltonian becomes N independent. 

The effective Hamiltonian can be expanded as before: 

HM = H<p+^j^p-\pMp)<p{-p) (4.8) 

n=3 "^Pi -^P" V i / 

The magnetic field, the propagator, and the vertices acquire now an explicit N dependence 
and we can write 

H = J^iJnA^-", (4.9) 

n=l 
n=0 

V^(^)(pi,...,p,) = 5^V;(^)(pi,...,p,)7V- (4.11) 



n=0 



4.1.2 Identities among the vertices of the effective Hamiltonian 

Let us now derive some important relations concerning the zero-momentum vertices Vq''\q, . . . , 0) 
and propagator Pq~ (0) [that below will be denoted by Vq (0)]. Since we will discuss only 
zero-momentum quantities, in this section we will not write the explicit momentum depen- 
dence. 

The leading behavior of Ticfriv'] for — ^ oo is obtained by considering the tree-level 
diagrams that arise from the perturbative expansion of (4.6) in powers of N . If we are only 
interested in zero-momentum correlations of the field </? at leading order in 1/A^ expansion^, 
in Ti, we can simply consider the zero-momentum vertices and in Z we need only to integrate 
over the five zero-momentum fields ^ . In the following we take the point of view of sec. 
3.1.1, in particular taking the short-hand notation 

/a = (7, a, a, T, r) 



''in order to consider 1/N corrections to the vertices, one have to include radiative contributions in which 
vertices appear at every momenta. 
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we have that, using the gap-equations (1.7,1.8,1.9,1.10), the vertices of the theory that in 
principle depend on /3 and are parametrized only by ttIq. In particular we can write^ 



rin) 



(4.12) 



In (4.12) /^(mo) and (3{ml) are given in (1.7,1.8,1.9,1.10). We are not interested in the 
explicit form of V^'' that is given in [56], otherwise in the following we need only the 
following relations: 



(n+l) 



dfA 



Al,---An 



= 



(4.13) 



in particular the (4.13) easily follows from the fact that the action of the auxiliary fields is 
linear in /3 [55]. In section 3.1.1, it was shown that the derivative of with respect to ml 
near the multicritical point is related to the zero mode va of the theory 



0(ml,pi) 



C{ml,pi) 
dp{ml,pi) 



{vA{ml,pi) + 0{ml,pi)U{'ml,Pi)) 



Si[mQ,Pi 



(4.14) 



is a no-morc specified vector while va is the zero mode at zero external momentum. O is 
a function of Pi and mg that goes to zero at the MCP (4.14). Since now we could neglect the 
explicit dependence on ttIq and Pi of O, ^a, va and C. Using (4.14) and (4.13) we recover 
the following equality 



d yin) 



dm?, ^1'-'^" 



(4.15) 



Identity (4.15) can also be written in the compact form 



dml 



N 
C 



J2 (^^ + ^^b) 



dn 



B 



-Y.'^Sl^A 



m 

dp 



(4.16) 



We begin by computing dTies/ dm^. We have 
dHe& 1 



dmQ 



J ^^*5'(*•^-^^/iV)*•^. 



-n 



NZ[(p] 

-ll[d^5i^.v-^VN)^^e 



■H 



((^^0). 



(4.17) 



The second term can be related to dHes/dif, by using identity (4.16). Indeed defining (•) as 



J 



-n 



(4.18) 



^Here and in the following we are neglecting the dependence on W oi pi. 
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we find 

A simple calculation allows also to write 

(4.20) 

Thus, we obtain the following representation 

— = 211 -J /xp . \ -L/xj/. \ 

dml C dip iV3/2rf^\ dml' iVV2 dip ^ dml' 

Now, let us show that the second term in the first line of (4.21) is of order 1/N and can 
therefore be dropped. Since v ■ dv/dm^ = 0, dv/drn^ projects over the massive subspace, i.e. 
the field ^ appearing in the mean value is a massive field. Thus, the leading contribution 
to the coefficient of 99" is due to tree-level diagrams with n lines corresponding to the field 
if and one line to the massive field ^ ■ dv/drn^. The A'"-dependence of a tree- level graph 
contributing to the coefficient of </?"■ is given by N'^/^~°- with 

00 

a^-J2ik-2)nk, (4.22) 

k=3 

where Uk is the number of k-leg vertices appearing in the diagram. Now we use the fact that, 
for a tree-level diagram. Next = '^ + Ylk=3i^~'^)''^k, where iVcxt is the number of external legs. 
Since N^xt = n -|- 1, we obtain a = A^ext/2 — 1 = (n — l)/2, independently of the diagram. 
Thus, 

(4.23) 

independently of the number of (p legs. This proves that at leading order in 1/N the second 
term can be dropped, while the third one must be kept. Consider now the first term of the 
second line of (4.21): 

§ E E y^T-AjB^A. ■ ■ ■ (4-24) 

n>2 ■ B{Ai} 

In (4.24) wc decompose ^ on the eigenvector of P (one of which is the critical mode va, 
while the other four are massive modes v^'^^ with m — 1, 2 • • - 4) 

4 

= (t;-*)t;^ + ^(^;M .^)^M (4.25) 

m=l 
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Considering only mass- less component in (4.24) we find the following (p"' contribution: 

^E^^ E yl^T-AjBVA.----A. (4-26) 

n>2 B{Ai} 

Otherwise inserting / (1 < / < n) massive vectors into (4.24) we have a contribution that 
looks like: 

n>2 B{Ai} 

A leading order graph that enters into the evaluation of (4.27) is a tree graph with / 
external massive legs and an arbitrary number of massless one. It behaves independently of 
the massless legs like iV~^/^+^ so that we recover an 0{N^) terms in (4.27) taking / = 1. As 
a consequence of (4.26-4.27) (up to l/N) we have that the first term in the second line of 
(4.21) goes as </7^ for (/? — > 0. The same considerations follows for dH/Ndp in (4.21), indeed 
following the same previous steps one can show that the leading order is 0{N'^) and that for 
</? — > this behaves like ip'^. In the following we will write: 

with Ok ~ Si near the MCP. Collecting all the previous results, at leading order in 1/iV we 
have 



dTies 1 dHes 1 dHes , J dv 
dml ^ C dip ~ ArV2 dip ^ "d^l 



n>2 

Since • dv/dml) ~ (/?^ for (/? — > we obtain a set of differential equation for the vertices 

7-j~r(n) T X (0) 

"'^0 _ -I- - <Jn,l T7(n+1) , V-^ JVl-fc T>(fc) 



^o'^^ -^'E7^^^o'' + «- (4.30) 



dml C " ■^^(^■-1)' 



where coherently with (4.28) we have implicitly assumed 01 = 02 = 0. In (4.30) we have 
defined 

' '*-^) - E4V«.o(l). (,31) 



Eq. (4.30) allows us to compute the behavior close to the MCP of the zero-momentum 
vertices Vq'^'' for k < J\f + 3. Using eq. (1.36) and the definition of the scahng fields (1.49) 
we have 

^ = X;^«^A^'-' + «i^^'^(-A^ + 2)A^^+^ 

k=l 

rrii+l 

+0 [A^^+^ \aJ-' T\uk,] , {A^ V}fe=i,..A/l , (4.32) 

L I, J i=l,2---M J 
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where = rriQ — m^^. All the terms in the first line of (4.32) go as near the MCP; the 
next to leading contribution (~ Uh^+^) is obtained considering the three terms of the second 
line of (4.32) with i even and = uj\f concerning the second term, and with Uk = uj\f in 
the third term. We will now show that for 2 < /c < jV + 3 zero-momentum vertices can be 
expanded close to the MCP as 

Vi'^ = a,A^^+3-^ + X: bf^u, + J2 bfuAJ-'"-' + . . . (4.33) 

j=l j=k-l 

For k = 2 the second term is absent. 

The (4.33) can be obtained starting from the following equality 

n^'^=^o-^ = /K,P.)^ (4.34) 

(that was demonstrated in [56] and recovered in chapter 1) and using recursively (4.30). We 
notice that for k = 1 the r.h.s. of (4.30) cancels in agreement with the condition Vq^^ = 
given by the gap equation. For k — 2 expansion (4.33) follows from eq. (4.32) and it is 
enough to identify a2 — f{0,0)a^^(^'^\N' + 2), b^^^ ~ j/(0, 0). For A; = 3 the expansion 
follows from eq. (4.30) and indeed all coefficients can be related with those appearing in the 
expansion of Vq'^^ and O. We will need the explicit expression for Ofe and b^''^ for j > k — 1. 
Applying eq. (4.30) recursively we obtain 

a, = (AT +!)! ,_, 
" {Af-k + 3)1 ^ ^ 



(k) 

where the second relation applies only for j > k — 1. For j < k — 1 the expression for bj is 
not so simple, indeed it will appear also contributions coming from O and higher derivative 
of /. However we are not interested in the exphcit expression of bf^^ ior j > k — 1 because 
this coefficients are not involved in the cancellations of the translated theory. This will be 
clarified in the next section 4.1.3. 



4.1.3 Translated zero- mode Hamiltonian 

In the previous section we have shown that all zero-momentum vertices with k < J\f+3 vanish 
at the MCP. We shall now perform a field redefinition in order to have ^^'^■'^"'"^^(0, . . . , 0) = 
for all values of the parameters. Proceeding in a similar way of what done in chapter 3 eq. 
(3.13) we write 

^{p) = ax{p) + kS{p). (4.37) 
In terms of the effective Hamiltonian has the expansion 

HMx] = ^X+^^"'(P)X(P)X(-P) (4.38) 

„=3 -^Pl -^Pn \ i J 
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The parameters a and k are determined by requiring that 

F^^+'^O, . . . , 0) = 0, (4.39) 
P"'(p) - P~'(0) = p2 + O(p^) (4.40) 
for all values of the parameters. Condition (4.39) is satisfied if 

^^^(Ar+2+,)(o,_._,0) = 0. (4.41) 

i=o ^' 

Since F(-^+^)(0, . . . , 0) vanishes at the MCP for N ^ oo, close to the MCP we can expand 

k as 

oo 

k = Y,kn[v^^^'\0,...,0)Y. (4.42) 



n=l 



The first term is easily computed: 



by using eq. (4.30) and eq. (4.33) it is easy to see that 

k = -^Am + 0[Am', Ui, 7V-1]. (4.44) 
Since k vanishes at the MCP for A?" ^ oo, we have at the MCP 

F^^'^O, . . . , 0) = \/(^')(0, . . . , 0) + 0{1/N), (4.45) 

for any j > A/' + 3. In principle one can observe that due to the leading behaviour of k (4.44) 

V will have the same general expansion of V^(^) (4.33), for j < TV + 2: 

1=1 i=j—l 

On the other hand in this case an important cancellation occurs (for j < J\f + 2). Indeed 
one can show that 6^ = for i > j and that aj = so that at leading order in Uh one has 

V ^ Uj-i. This implies that the translated vertices vanishes faster as — > 0. 
We start from 

V^'^^a^ J2 ^V^,("+'■^ + 0(A;^+^-^A^-l). ^^^^^ 

n=0 

The leading contribution independent of the scaling fields is given by [using the notations 
given in (4.33)] 



n=0 



"-^^^^^^{Ck + AJ^-^'-^ + 0[(A^)^+^-^ = 0[(A^)^+^-^ AT-i], 



(4.48) 
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where we have used eqs. (4.36) and (4.44). As for the contributions proportional to the 
scahng fields Ui, no cancellation occurs for i < j — 2. If i > j — 1 the coefficient of is given 

by 



E ^^^1— (A„.)'-^-"+' + 0[(A^)^-^+2,Ar-i] = 

77/. 

n=0 



a- 



(^-J + 1) 



— -(Cfc + A„)^-^+^ + 0[(A^)-^+2, AT-i] = 0[(A^)-^+^ 7V-i](4.49) 



The previous equations (4.48-4.49) show that = and 6^ = if i > j in (4.46) confirming 
the starting claim. Otherwise other cancellations appear at next to leading order, so that 
the general structure of the translated vertices is not given by the corrections of (4.46). 
This is investigated in the next sub-section 4.1.4 and in 4.1.5 where the structure of k is 
investigated in detail, in particular it is shown that the A^^ correction of k in (4.44) is 
due only to contribution coming from C. The strategy of the next sub-section 4.1.4 is 
to generalize the set of differential equations for V^^^ (4.30) to to the translated one v'^^\ 
Near the MCP this set of differential equations decouple so that we can easily recover the 
cancellations discovered in this section. Otherwise using the new set of equations one is 
able to discover other cancellations. One can think to have in mind the final results for the 
translated vertices (4.63,4.72, 4.73) and in particular the translated vertices for the A/" = 3 
case (4.74,4.75,4.76) that will be used in the final section of this chapter 4.2, where 1/N 
corrections for the scaling fields and the scaling variables will be explicitly built for the 
A/" = 3 case. 



4.1.4 Structure of the translated vertices: another derivation 



For convenience in this section we will take a = 1 in eq. (4.37) (at the end V must be 
rescaled by a^). One can try to solve in a more direct way the structure of the translated 

vertices V . In a similar way of what done in previous sections we define 



Z[X\ = 



J &^S{ip ■ V - 
N ^ Z[0] 



(4.50) 
(4.51) 



and we follow the same steps of sec 4.1.2. The first change is in equation (4.17). Using an 
obvious definition for (■) and (4.20) we have: 



dU 



eff 



dmg 



1 dn 



1 dU 



cff 



N dx 



dv , dHcs d/c 
) + 



drriQ 



dx dm5 



+ 



X 







(4.52) 



On the other hand (4.19) does not change apart the substitution ip —>■ x + k (notice however 
that the explicit k dependence disappear in (4.52) because the subtraction in the second 
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line). So we generalize (4.29) in: 

dHeS 



d^eff _ fl dk \d7Yeff Vo^'\, V^^fv I M 



1 dHcs , _dv_ 



N dx 



+0 



(4.53) 



Using (4.31) we are now able to generalize (4.30) for the translated vertices: 



dml 



+ n\ 



(0)y7("-1) (0)T7(n-2) 



(n-2)! (n-3)! 



{n-l)\ 



+ 3 



+ --. + (n + l)dy 



(0) t7(i) 



(0)T7(n+l) 



n! 



Oj with J > 2 + 



(n-2)! 

+ 6|^ + --- + (n + l)ncilF(^^ 



On + {n+ l)on+ik +{n + 2){n+ l)o„+2y + 0{k^) 



t7(2) 



In (4.54) we have defined 



ks{ml) 



1 dfc 



C dml 



k. is studied in detail in sec. 4.1.5. There we have shown that 



ks{ml) = 



AoC 



(4.54) 
(4.55) 

(4.56) 
(4.57) 



in (4.56) the term cancels in a non trivial way, while in (4.57) also the A^^ term 
disappears. 

Wc have just noticed that V^^^ in principle has the same structure of V^^^ (4.46) near the 
MCP. Otherwise we notice that at leading order in Uh due to the relations (4.56) and (4.57), 
eq. (4.54) becomes 

'^'^ -U) 

= ^ V^^'^Uj_i (4.58) 



dV' 



dml 



so that we recover the cancellations discovered in the previous section 4.1.3: 



6?^ = 



Gj = 0. 



(4.59) 
(4.60) 
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We can complete (4.59) with 

6^ = 6?) (4.61) 
that follows from the fact that k goes to zero at the MCP. If n = 1 (4.54) becomes: 

= -^V-o''^' - 20.* - -^Vr F'" = 6<'>«. (4.62) 

where we have used the gap equation (in the following uq = Uh). Eq. (4.62) is expected if 
compared with similar results for higher order vertices eq. (4.58). 

Now we consider corrections to (4.58). In particular we have to compute at what order 

(in Am) Uj enters into V^^"^ ior j > k — 1, i.e. we want to compute x^^^ ^-^d Xi^ ( j < ^ < J^) 
in the following expression 

V^'^ = J2 ^'u, + Vi^^^j'' + X: Vl^KArJ'. (4.63) 

1=1 i=j 

Writing (4.54) for n = jV + 2 (taking into account the constraint y^"^"*"^^ — 0) we have: 

(ksiml) - eW^"^"-'^ = 4°)(Ar + 2)(AA+l)(l + 0(A;))F^^+'^ 

+c^^\Af + 2){Af+l)Af{l + 0{k))V^^^ + ■■■ (4.64) 



this show the importance of cancellations that happen in (4.56) and (4.57), indeed the 

previous equations tell us that the A^^ term in v^'^^^^ disappears according with (4.60). 
However in principle solving (4.64) one is able to compute the following expression 

AT 



i=l 



that can be recovered simply using the leading behaviour of the translated field. Then one 
can try to find . If we consider n — 1 case of (4.54) [defining r]{m'^ = ks{rn^ — c^q^ k{m^'^] 

f(i) 



where a is an appropriate constant that can be obtained including all the proper terms 
in (4.54). Remembering that F^^^ = f{ml)si, we easily obtain that the correction to the 

leading behavior V^^^ ~ Uh is of order A^ and comes from f(ml) so that 

xf^=z + l x^i)=AA+3. (4.67) 

Considering the n — 2 case of (4.54): 

VV^'^ = ^ + oJl + Oik)] + 2crF^'^ + 2kc^^^V^'\ (4.68) 

QTTT-Q L J 
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The leading corrections comes from 02 terms so we find 



(4.69) 



Using (4.54) is then quite easy to compute the terms that appear in the expansion (4.63). In- 
deed we can rewrite the differential equation (4.54) including only the terms we are interested 
in as 



dV, 



JV+3 

^ _L + A 377(^+1) I ^ A j+l-nTrU) 

j=n+2 



(4.70) 



— U) 

Integrating the previous equation, using Vq ~ Uj^i, we get 



-'0 -I 4 ^ J + 3 - n 



j=n+l 

_^A(+3_t7(-^+3) A M+5-n 

Al +0 — n 



(4.71) 



Using the previous equation (4.71) we are able to give our final result (for j > 2) 

X? = Min[7-,4] 



J 
U) 

i 

iJ) 



Min[i, i + 3 — j] j ^ i 
Min[A/' + 2,A^ + 5-i]. 



(4.72) 
(4.73) 



To finish this section we give an explicit expression of the translated vertices for the 
A/" = 3 case that will be studied in detail in the rest of this chapter: 



V 
V 



(2) 



= hf\^ + Vf «2A„2 + vf^u^Aj + V(2) A„^ 



;(2). 



T7(3) z,(3)„. , ^(3)„_ ^ v(3).,.A 3 ^ y(3)^ 5 



V 



(4) 



= hf\^ + + ht\, + V(^) A^^ 



2 



(4.74) 

(4.75) 
(4.76) 
(4.77) 



4.1.5 k{ml) near the MCP 

In this section we want to investigate the structure of k{m^) (4.37) as a power of A^ in 
order to recover the claims of the previous section [see eqs. (4.56) and (4.57)]. This is the 
fundamental issue in order to investigate the algebraic structure of the translated vertices 
In particular we want to show that the A^ and A^^ terms that enter into the expression 
of ks{rn^) near the MCP are due only to the C{m^) expansion near the MCP. We want to 
investigate only the structure in A^ [notice that the scaling fields Ui trivially appear linearly 
in k{m^y\ so in the following we can think to put Ui = 0. This allows us to think to the 
expansion (4.78) as an expansion in Am- Otherwise the considerations are general in the 
MCCL, the expression (4.78) being an expansion in Uh- 
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An easy computation gives us: 

+0 (A^^) (4.78) 

In this subsection we limit ourself to A/" > 1. In this case we have that d/5/dmo^ is taken 
into account in the A^^ term in (4.78). This allows us to neglect in (4.30) the terms that 
are order si near the MCP 



dv: 



in) 



-Vr 



(n+1) 



(4.79) 



dmo^ C{ml) ° 

where the dots mean that in the following the corrections can be neglected. Using (4.79) 
and (4.78) we get a very easy expression 



dk{ml) 



1 T/^ 

^ + ci°)(AA2 + 3AA+2)-^ 



{^f+l) 



+o(aj,1 



Vr 



■{Af+2) 



{M+3) 



Vr 



(4.80) 



In the previous expression we have taken the convention used in this work 
As a straightforward consequence of (4.79), we have that 



Vr 



Vr 



iAf+2) 



1 



T/(A/+3) 

so that we recover the wanted result 

dk{ml) 



A^2 + o(Aj), 



JO) 



+ 



ks{ml) = 



C{ml) Ciml 
dk{ml) 



2\2 



(4.81) 
(4.82) 

(4.83) 



+ 



drriQ 



C(m^)2 



(4.84) 



and finally 



k,iml)-kimlf = OiAj). 



(4.85) 



In (4.83) the A^ is due only to C{mQ), this disappears with the A^^ term in the (4.85). 
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4.2 Scaling Fields and the large limit 

In the previous section we have obtained a detailed description of the zero mode Hamiltonian 
near the MCP. We are now in a position to apply the results of chapter 2 concerning the 
multicritical points. In particular we have verified that at the AT-th order multicritical point 
all the vertices Vq''\o) (with j = 1, - ■ -Af + 2) go to zero in a proper way that has been 
careful characterized above.^ This is a necessary condition in order to apply the results of 
chapter 2 concerning the crossover to multicritical points. 

Taking a pedagogical point of view we will give a detail description of the A/" = 3 case 
for which one can apply the results of sec. 2.2. Following the same step of sec. 3.2, we want 
to investigate i) the scaling variables that describe the crossover between classical behavior 
to Mean Field behavior eq. (2.67), ii) the correction to the critical point Pi = Pic and P = Pc 
given by (1.15) and in) the 1/N corrections to the classical scahng fields defined in (1.48) 
and to the classical scaling relations (1.49). We have just pointed out at the beginning of 
this chapter that the cancellations we have found in the previous sections are very important 
to solve the point in). 



4.2.1 J\f = 3 case in two dimension 

For the A/" = 3 the translated vertices near the MCP has been explicitly given in eqs. (4.74), 
(4.75) and (4.76). Using them we are now in a position to apply the result of sec. 2.2, in 
particular the results for the Hamiltonian (2.63) can be used with the following identification 
(for i = 1,2,3,4) 

Uoi^V^\o) u,^^^^^. (4.86) 

Notice in particular that the normalization conditions assumed in eq. (2.63) are satisfied due 

to the definition of uq (4.86) and to eq. (4.37) that ensures V (0) = and K{p) ^ p^. 
We will consider only the two dimensions. Otherwise, as pointed out for jV = 1, for d > 2 
one has to translate the theory to impose that the renormalized five legs vertex is zero. The 
scaling variable are then simply given by 

Xh = - — (H-K) (4.87) 

a 

xi = ^{V^'\0) - u,^) (4.88) 

X2 = ^{V^'\0) - u,2) (4.89) 

^3 = ^{V^'\0)-u,3) (4.90) 

where Ud have been computed using a lattice regularization in chapter 4. We are not 
interested in the exact form of this variable but we will develop them in the most general 



^We are neglecting corrections and the countcrtcrms that one needs to include in the effective action 
in order to regularize the theory. This contributions will be deleted including corrections to the N = oo 
multicritical point. 
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way around the MCP. Notice that in the computation of Ud and Ucs we are not able to give 
terms: indeed it would require an hard fine tuning not available until now. On the 
other hand he and Ud can be obtained imposing Z2 symmetry at the transition. 

Scaling fields equations. Using the results of the tree-level translated theory we have: 

+N^/\Ho + HiAm + •••) + N^^^ ■ ■ ■ (4.91) 
X, = N'{b?m + viX2Arn' + viX,Aj + V^'^Aj + ...) 

+N{Uo + U,A^ + ■■■) (4.92) 
X2 = N'/-' {bf^u, + bfK, + vjt^usAj + V(^) + • • • ) 

+N'/\Do + D,A,n + ---) (4.93) 

+(To + TiA^) (4.94) 

In the previous expression we do not have considered the possibility of cancellations at one 
loop order, this is not necessary for the considerations we are going to do. Otherwise if one 
is searching for particular symmetries some cancellation could be model dependent.^ 

Wc notice that at the classical level, the scaling fields Ui are not defined in a unique way. 
Indeed we are able to re-define: 

i-l 

u[ = kiUi + ^ fjUj (4.95) 
i=o 

so that we can eventually resum at classical level the 6^'^ contribution in (4.91-4.94). On the 

other hand in order to cancel the V^*^*^ terms in (4.91-4.94) we have to include 1/N correction 
in the definition of the scaling fields. 

Correction to the critical point. Solving (4.92-4.94) in Ui and putting the result into 
(4.91) we argue that ~ mo/N'^/^ in the MCCL. Using then (4.92,4.93,4.94) we find: 

-^-ir + ^(]V^j ^ = 1,2 (4.96) 



with: 



„(0) / , 



W = (4.98) 

^0 "000 



''For instance if = 1 and the action is Z2 symmetric, then in chapter 5 we will show that some 
cancellations are present also at one loop order. These cancellations give the obvious result /ic = (i.e. 
explicit symmetry is not broken). 
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= (4.99) 
= (^"^-^"^^7/ + ^°-^° (4.100) 



.(3) 
^2 



4°' = -= (4.101) 

O3 ' 

The previous equations (4.96-4.101) are the leading correction to the position of the critical 
point (notice they are independent of Xi). Otherwise the position of the critical point (ud) 
can be obtained studying (4.91-4.94) with Xi = 0. It is not hard to see that for Ud we have 
the following developments in powers of N~^^^: 

^mc - ^1/5 + ^2/5 + ^ K^.^y^^} 

(0) (1) (2) 

(0) (1) (2) 

«2c = — + ^ + ^ + --- (4.104) 

(0) (1) (2) 

= ^ + ^ + ||, + - (4.105) 

in particular moc = ''tio and ■u-^'' = ■Uj-^'' (4.98-4.101). 

Redefinitions of the scaling fields. Now we take 7^ in (4.91-4.94). Defining 
m = — Ui = Ui — Uic, taking into account (4.102-4.105) the (4.91-4.94) now become: 



■y^, = A^^i + A^^2+A^^3 + 5aooo^m+--- (4.106) 

= ^1 ^1+ ^2/5 ^2+ ^3/5 i^3 + 5 m + --- (4.107) 

= b[>u, + b'2>U2+ ^3^, % + 5 ^^^^ m + --- (4.108) 

— = b\ 'ui + b'2 'U2 + 6^ ^-»3 + ^3^g m + • • • (4.109) 

where dots indicate sub-leading contributions. In principle one can solve the previous equa- 
tions (4.106-4.109) obtaining the following expansion: 

m = ^ MijkiX3'x2^xi''xh^ (4.110) 

i,j,k,l 

J2 U^x^'x^'x.^xu' (4.111) 



Up 



with ^ when N 00. Following chapter 3 and [56], then one tries to solve the 

previous mixing defining new scaling fields that have to take into account 1/N corrections re- 
spect the mean-field definition. As a first step we compute (at leading order) the linear terms 
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(in which only one Xi appears) in (4.110-4.111). We easily obtain the following expressions: 

Ui 



U(^) ^ + ^(1) ^2 (1) jc^ (1) JCf^ 



U2 



77(1) = /Vi9/i07y(i) - ^"^'^^'^ \ oi ^ 

- ^^0^00- (5) (3) ,(2) ,(3) ,(2) +^1;VV10 

"000 "2 "1 ^'2 ^\ ^ 

y.iX) = ^8/5^(1) mo,3V(2) ^0^3^(2) / 1 \ 

- "1000 - (5) (4) ,(2) ,(4), (2) + ^ I ^1/10 J 

"000^3 "1 "3 "1 ^ ^ 



^ 77(2)^ ,7/(2) J^,,/(2)J^, 7,(2) J^, 
^2 ^ ^3/2 ^ '^a;3 ^8/5 ^ ^^/. ^^5/2 ^ 

r27/(2) _ ^1 



^xi - ^V^ooio- ^(2)^(3) +^ \^ ^1/10 



2 

1,(3)^ 3-i;(2) 7,(2)i;(3)^ 3 ™ 3 

7/(2) = /V8/5zy(2) _ ^1 "^0'^ "^3 - bi V3 _ mp, ^fe(3)v(2) _ ^,(2)v(3)^ 

"^3 - "1000- .(2) ,(3) ,(4) (5) ,(2) ,(3) ,(4) ^ J 



+0 



^1 ^^2 "3 "000 "1 "2 " 



3 



^xn - "0001 - ^(5) ,(1), (2) ,(3) ^(5),(1),(3) 7VV10 j ^^■^^'^^ 

"000 "0 ^1 "2 "000 "0 "2 ^ ^ 

- ""0010 5 (5).(2).(4)^j(2),(3),(4) ,(2),(4) + " I JVl/10 

"J"ooo"i "3 "1 "2 "3 "1 "3 ^ 



(4.115) 
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= TVWio^v^ooio = 7^ ^ + 



-XI 



OU2 lllQc Uqqq \ 



Kirn) = N^/^MwOO^ T^J-^ + O^ ^ 



-X3 



Z^M ^ 7V^^/^°>loooi = rTM^-TTTT + O (^) (4.116) 



We see that at leading order in 1/A^ (4.114) becomes: 

u,^u,r.uif^ (4.117) 

so that the new scaUng field (u^) does not need any corrections with respect the older one 
{U3). The same happen for U2: 

= (4.118) 

In the expression for ui (4.112) we see that the terms in X2 and X3 are sub-leading with 
respect to xi. So one has to include correction to this scaling fields in the following way 
[using (4.117) and (4.118)]: 

_ , -I- Ux2 ~ -I- , 7/(2)^1 /A iin\ 

Ui = Ui- -T7^^U2 - TT^^^S ^ ^L' TT^ (4.119) 



Notice that in the previous expression Ux^J /Ux^J and U^] /lA^} need to be computed up to 
corrections of order A^^-^/"^ and N~'^^^ with respect the value reported in (4.112,4.113,4.114). 
This implies that if one wants to compute exactly the mixing coefficients in (4.119) one needs 
to include in (4.106-4.109) almost N^"^^^ corrections to the leading behaviour, in particular 
also radiative corrections will enter in the effective computation. Now we have to discuss the 
effect of having discarded the non-linear terms in (4.110) and (4.111). But it is easy to see 
that [including contributions that correct (4.107)] for instance the x^"^ contribution scales as 
l/A^22/^°. This shows that this contributions can be neglected. 

Now we pass to discuss the Uh field. Using the gap-equation and (4.91) we include 
corrections to (4.106) in the following way 

+ T7^«i + T7^«2 + ^77777% + + —TjKmm (4. 120) 



Notice that in the previous expansion (4.120), the M3M3 term (~ xs^N~^) is not present; 
indeed it has been adsorbed in the definition of Uh- At leading order in N we have (4.91- 
4.94) 

^1,^2,^3 = V;'^moe^V^'^moe^V;i')moc' (4.121) 
nil^ = 6miy^^+H, (4.122) 
n^^^ = I5mty^ (4.123) 
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while in order to compute completely the mixing-coefficients many other contributions (in 
l/iV^°) are needed. Using then (4.120) and (4.116) we can define the new scaling field Uh 



-17(1)7 /{"O \ 1 /-ly 7y(^) 7y(^) 7y(^) n-/^^h/^"^^\ 

rim l4x2 \_ J- / riiLAx3 n2L4x3 ri^LAx^ rim iAxa \_ 



Ar3/5 J "2 ^g) y ]^ + ^6/5 + ^ Ar4/5 ) ^3 



-1,(2) 



(^^(3))2 7V2/5 

SO that we can write: 



ul (4.124) 



4.3 Conclusions 

In this chapter we have investigate the critical behavior of 0{N) models with a general- 
ized ferromagnetic nearest neighbour interaction depending on J\f tunable parameters that 
undergo a phase transition at finite temperature (see chapter 1). In sec. 4.1 we have stud- 
ied in detail the algebraic structure of the zero Mode Hamiltonian near the critical point, 
parametrizing the vertices using the scaling fields of the N = oo theory (chapter 1). Using 
the results of chapter 2 in two dimension, in sec. 4.2 we have applied the theory of Mul- 
ticritical Crossover Limit to the specific case M — 3. The approach we have given is only 
algebraic (for instance we do not have taken care of the explicit form of but we have 
taken the most general behavior near the MCP) but shows the steps that one need to follow 
in order to study the critical limit when iV — > oo. In particular it shows the definition of 
the scaling fields and corrections to the critical point that in principle could be compared 
with available numerical results (in a similar way of what reported in sec. 3.5). Using the 
scheme presented in this chapter we believe several physical system could be investigate. 
For instance finite temperature Yukawa model with Nf fermions (for Nf — > oo) has been 
investigate in chapter 5. There we have shown as the system, for every Nf finite, undergoes 
an Ising phase transition that crosses on a Mean Field phase transition for Nf — oo in the 
same way of the J\f — 1 models introduced in chapter 1 and studied in chapter 3. It is usually 
claimed that if one tune also a chemical potential the phase diagram exhibit also a tricritical 
point. We plan to apply the method developed in this chapter for the investigation of this 
multicritical point. 



®The explicit form of Uic has been computed in chapter 2 



Chapter 5 



Fermionic Models with chiral 
symmetry 



In the previous chapters we have investigated several problems related to the 1/N expansion 
in models with global symmetry. In order to investigate such kinds of problems we have used 
generalized Heisenberg models. However we believe that the presentation given is general 
because in the effective action for the zero mode (translated or not) noone symmetries was 
assumed, but the most general expansions was taken. ^ For this reason we have tried to 
investigate the large limit of other Infra Red divergent models [49]. In this chapter we 
present our analysis of fermionic models with global symmetry at finite temperature. This 
kind of studies could be significant in order to investigate the nature of the QCD phase 
diagram. The finite-temperature transition in QCD has been extensively studied in the last 
twenty years and is becoming increasingly important because of the recent experimental 
progress in the physics of ultrarelativistic heavy-ion collisions. Some general features of 
the transition, which is associated with the restoration of chiral symmetry, can be studied in 
dimensionally- reduced three-dimensional models [44, 45]. However, a detailed understanding 
requires a direct analysis in QCD. Being the phenomenon intrinsically nonperturbative, our 
present knowledge comes mainly from numerical simulations [46, 47]. Due to the many 
technical difficulties — finite-size effects, proper inclusion of fermions, etc. — results are not 
yet conclusive and thus it is worthwhile to study simplified models that show the same basic 
features but are significantly simpler. 



^This is not completely correct. We have seen that in 0{N) models the sign of the effective four leg 
vertex at the Critical Point has the same sign of the mass term of the zero mode near the CP [the claim 
follows from eq. (3.12)]. In principle this fact could not necessary holds for other models which share with 
Heisenberg models the existence of a zero mode. 
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5.1 The model 

In this chapter we shall consider a Yukawa model in which Nf fermions are coupled with a 
scalar field through a Yukawa interaction. 

Nf 



DNf J d'^+^x (i {dcpf + f 0' + ^0') + E / d'^'^^/ {^ + 9<l> + m) 

(5.1) 



where ^ = Yl'jj^l^iln^n Itii 1^ — 1 are the generator of a Clifford Algebra 

satisfying the anti-commutation relations {'jfj,,'ju} = 2(5^^,^1 where tr7^ = D,^ and 1 is the 
identity DxD matrix. The integration is over M*^ x [0, T"^], and A > to ensure the stability 
of the quartic potential. Along the thermal direction we take periodic boundary conditions 
for the bosonic field (j) and antiperiodic ones for the fermionic fields -0/ 

0(x, 0) = 0(x, T-') ^(x, 0) = -^(x, T-i) ^(x, 0) = -^(x, T-i). (5.2) 

The theory must be properly regularized. We shall consider a sharp-cutoff regularization, 
restricting the momentum integrations in the spatial directions to p < A. However, the 
discussion presented here can be extended without difficulty to any other regularization that 
maintains at least a remnant of chiral symmetry. For instance an explicit lattice calculation 
involving staggered fermions has been considered. However all the computations strictly 
follow the presentation of this chapter without any new consideration;^ in particular the 
relations that protect the symmetry of the theory arc satisfied also for chiral fermions on the 
lattice (in the following this fact will guarantee that = for every Nf). 

Given a representation of the Clifford Algebra in even dimension (see e.g. [24]), 7c is 
introduced 

7c := 7i72 • • • 7(i, (5-3) 

so that it anti-commutes with all the generators of the algebra {jc , 7/*} — 0, fj, — 1,2 •■■ d. 
On the other hand an odd dimensional representation can be built by using the standard 
d — 1-dimensional representation, plus 7c (5.3) multiplied by a proper C number 



^(oH^^^ /^-(2rf) if//<2d+l 

7.(2ci+l) = | ^^^^^(2d) if^ = 2d+l. 



The interaction (5.1) has a discrete symmetry for M — 0. If 0? is even this can be realized 
by the following transformations 

'0^7c'0 '0^-'07c (5.4) 



^One can chooses a representation so that D = 2'^/'^ if d is even, D = 2(''+i)/2 if d is odd 

■^Thc only significant differences are related to the doubling problem (partially solved by staggered 
fermions) that gives additional degrees of freedom. 
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In odd dimension 7^ is trivial 7^ = (— Actually the previous symmetry (5.4) can be 
realized using a ji coupled with a reflection in the Euclidean space along the i-th direction 
[24]: 

i/j^'jiip i/j^-i/j'yi (f) ^ -(f) 

{xi, ■ ■ ■ Xi, ■ ■ ■ Xd+i} {xi, Xi,---Xd+i} (5.5) 

We present an unified picture using (5.5) with i — d + 1. We are interested to study the 
reduced-dimensional action which is obtained integrating out the thermal component so 
that the minus sign appearing in (5.5) can be neglected. The previous symmetries (5.4-5.5) 
protect by fermions mass generation and can be investigate studying the expectation value 
of the bosonic field < (j) > that will be used as an order parameter to investigate the phase 
diagram of the model (5.1). 

In the limit Nf — * 00 this model can be solved analytically and one finds that there is 
a range of parameters in which it shows a transition analogous to that observed in QCD 
[48, 49]. It separates a low-temperature phase in which chiral symmetry is broken from a 
high-temperature phase in which chiral symmetry is restored. For Nf — 00 this transition 
shows mean-field behavior, in contrast with general arguments that predict the transition to 
belong to the Ising universality class. This apparent contradiction was explained in Ref. [50] 
where, by means of scaling arguments, it was shown that the width of the Ising critical 
region scales as a power of l/iV/, so that only mean- field behavior can be observed in the 
limit Nf — 00. 

This picture is strictly like what we have observed in a generalized 0{N) a model in 
the first part of this work: for finite values of N the transition was expected to be in 
the Ising universality class, while the N = 00 solution predicted mean-field behavior. In 
chap. 3 we performed a detailed calculation of the 1/N corrections, explaining the observed 
behavior in terms of a critical-region suppression. The analytic technique presented in this 
work can be applied to model (5.1). It allows us to obtain an analytic description of the 
crossover from mean-field to Ising behavior that occurs when Nj is large and to extend the 
discussion of Ref. [50] to the case M 7^ 0. More importantly, we are able to show that the 
phenomenon is universal. In field-theoretical terms, it can be characterized as a crossover 
between two fixed points: the Gaussian fixed point and the Ising fixed point. This implies 
that quantitative predictions for model (5.1) can be obtained in completely different settings. 
One can use field theory and compute the crossover functions by resumming the pcrturbative 
series [29, 30, 31, 32]. Alternatively, one can use the fact that the field-theoretical crossover 
is equivalent to the critical crossover that occurs in models with medium-range interactions 
[33, 34, 32, 35, 26]. This allows one to use the wealth of results available for these spin 
systems [33, 36, 37, 34, 32, 26]. Finally, we should note that the phenomenon is quite 
general and occurs in any situation in which there is a crossover from the Gaussian fixed 
point to a nonclassical stable fixed point. For instance, similar considerations have been 
recently presented for finite-temperature QGD in some very specific limit [58], while for 
other models the scenario presented above is questionable [43] . 

In sec. 5.2 we review the behavior in the limit Nf = 00. In sec. 5.3 we consider the 1/A^/ 
fiuctuations and determine the effective theory of the excitations that are responsible for 
the Ising behavior at the critical point. These modes are described by an effective weakly- 
coupled 0^ Hamiltonian. These considerations are applied to the Yukawa model in sec. 5.5.2 
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and in sec. 5.5.3. We determine the relevant scaling variables and show how to compute 
the crossover behavior of the correlation functions. In the appendix B we carefully discuss 
the relations among medium-ranged spin models, field theory, and the Yukawa model. The 
presentation given strictly follows Ref . [59] and Ref . [60] . 



5.2 Behavior for Nf = oo 

The solution of the model for Nf = oo is quite standard. We briefly summarize here the 
main steps, following the presentation of Ref. [49] . As a first step we integrate the fermionic 
fields obtaining an effective action 5efr[0] given by 



-DNfS 



Nf 



JJd^yd^/e-'^[^''^'^], 
/=i 



(5.6) 



where 



5eff[(/.] = J d'^+^x(i(a0)2 + + ^04^ _ 1 y d'^+ix tr log + ^0 + m) . (5.7) 

For Nf ^ oo one can expand around the saddle point (j) — (j), that is determined by the gap 
equation 5Ses[(f)\/S(l) — 



A , 
um H — 
6 



Jp<A {27ry P^ + ujl + {m + My 



where we define the frequencies lo^ = (2n + l)7rT', and 
The action corresponding to a saddle-point solution m is: 



nez 



log " ^ 



(27r) 



p2 + a;2 



(5.8) 



(5.9) 



where we have added a mass-independent counterterm to regularize the sum [49]. Such a 
quantity has been chosen so that the action for M — m — vanishes. Summations can be 
done analytically using the identity^ 



E 



X 



(n + 1/2)2 + a;2 



= 7rtanh(7ra;)=7r(l- ^^J^ J . 



(5.10) 



''The (5.10) can be proved considering the identity = tt dz (^2+i/2)^+x^ ' that holds considering for 
instance C a rectangular in the complex plain with the vertex on semi-integer number [i.e. Vi — (m^ + 
1/2, mi + 1/2), mi,ni G Z and i = 1, 2 • • - 4], and taking the limit mi,ni ±oo, with the proper sign so that 
C encloses all the plane. Using Cauchy theorem one can evaluate the series (5.10) computing the residues of 
the previous integrand at z = —l/2±ix. 
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The gap equation can then be written as 

p{m) = (171 + M)g{m + M,T), 
where the functions p(m) and Q [x, T) are defined by 

A 



(5.11) 



p(m) 



11171+ —m , 
6 

r d'^p 1 



+ X^ \ 2 e^/^+x^/T 



(5.12) 
(5.13) 



Analogously we can rewrite eq. (5.9) as 



A 



5eff(m,M,r) = ^m' + -m 



4! 



-T 



f ^ 

Jp<A c^^r 



log cosh 



^yp^ + {m + My 



2T 



P 

log cosh — - 
6 2T 



(5.14) 



Using eqs. (5.11) and (5.14) we can determine the phase diagram of the model. Given JI and 
A, for each value of T and M we determine the solutions m of the gap equations. When the 
solutions are more than one, the physical one is that with the lowest action Ses{m, M,T). 
Note that eqs. (5.11) and (5.14) are invariant under the transformations m — > — m and 
M — > —M. Thus, we can limit our study to the case M > 0. In general, we can find either 
one solution or three different solutions mo, m+, and m_ with m_ < itlq < m+ (for some 
specific values of the parameters two of them may coincide). The previous claim can be 
proved using the following identities (see app. A) 



—xg{x,T) > 
dx 



dx'^ 



xg{x,T)e{x) < 0, 



(5.15) 



[0{x) is the step function] and straightforward consideration on the gap equation (5.11). 
There are four different possibihties: 



(a) liji> g{0, 0) = CoA"^-^ with 



Co = 



2V/^(d-i)r(^)]"\ 



(5.16) 



then, for every M > 0, there is only one solution m+ > 0; for M = we have m+ — 0. 

(b) If < /I < CqA*^"^, there exists a critical temperature Tc^jl) (see fig. 5.1). For T > Tciji) 
and any M there is only one solution m+ > (for M = we have m_|_ = 0). For 
T < Tciji) and M < M there are three solutions mo, m+, and m_ with m_ < mo < m+ 
and m+ > and m_ < 0. For T < TcIJl) and M > M there is only one solution 
corresponding to m_|_. The physical solution is always m+ so that M has no physical 
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-0.4 




0.4 



Figure 5.1: Plot of p(m) and m^(m) for CqA'^ ^ > > taking different temperature. Tfie 
symmetry is broken at low temperature (m+ 7^ 0) 



meaning. Moreover, for T < Tc(/i) and M = 0, m_|_ > 0. The critical temperature can 
be computed from the following relation: 



where u!c,n = (2?^ + i)T^T(. and Tc{ii) — tc{iJ>)A. For // ^ 0, we have Tc(//) 00. 
(c) If -Ci < /7 < (see fig. 5.2 left), with^ 



there is a critical mass M such that there are three solutions for M < M, two of 
them coincide for M = M, while for M > M the only solution is m_|_. The physical 
solution — the one with the lowest action — is always m+ so that M has no physical 
meaning. Note that m+ > for M = 0. 

(d) For fl < —Ci (see fig. 5.1 right) there are three solutions for all values of T and M. 
The relevant solution is always m+ > 0. 

In order to compute the more stable solutions (if more than one are present) one can use 
the previous consideration, the figs. 5.2 5.1 and^ 



^Ci is the solution of the equation p{—x) = limM—>oo MQ{M,0) with x = {2Ci/ Xy^"^ . The value x 
corresponds to the position of a maximum of p(m) for fl — — Ci. 

^c(M, T) can easily be computed using eq. (5.9), in particular c(0, T) — 0. 




(5.17) 




(5.18) 




(5.19) 
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Figure 5.2: Plot of p(m) and {m+M)Q{m+M) for several M and with JI < (the intersection 
points are the solutions of the gap equation 5.11). (Left) Case (c): there could be one or 
three solution, however m+ > is always present and is always the more stable. (Right) 
Case (d): there are always three solutions. 




Figure 5.3: S^s for M = and M > 



to evaluate the more stable solution. Using the symmetry of p(a:) and xQ{x) we get that 
in case (a) chiral symmetry is never broken, while in cases (c) and (d) chiral symmetry is 
never restored. Thus, the only case of interest — and the only one we shall consider in the 
following — is case (b), in which there is a chiraiiy-symmetric high-temperature phase and a 
low-temperature phase in which chiral symmetry is broken (see fig. 5.3). 
The nature of the transition is easily determined. We expand 

g{x,T) = Y,g^^x^^{T-T,Y\ (5.20) 

m,n 

where 



goo = 



n62 



d'^p 1 



^X<A(27r")^(p2 + ^2j2 
1 /■ d'^p eP'^^ 



;<A (27r)<^p2 + ^2^' 
d'^p - 



c,n 



gio = 



T2 {27rr {ePl^^ + T 
d'^p 1 



2' 



"^^gX<^(27r)^(p^ + 02' ^^-^^^ 
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and cUc.n = (2n + IjvrTc. Since goo = ^(0, = n [eq. (5.17)], the gap equation (5.11) 
becomes: 

^m^ = ]lM + goi{T-n){m + M)+gio{m + Mf + --- (5.22) 
where we have neglected subleading terms in m + M and T — Tc- Defining 

Uu^-^M u.^-^(T-T^), (5.23) 

6gio - A 6g-io - A 

and taking the hmit Uh,Ut ^ at fixed x = Ut/u^f/^ we obtain the equation of state 

m = Uh'/^fix) (5.24) 
= f{xf + xf{x) + l. (5.25) 

Note that the prefactor of T — Tc in Ut is always positive to ensure the existence of only one 
solution for T > Tc. Such an equation is exactly the mean-field equation of state that relates 
magnetization ip, magnetic field h, and reduced temperature t. Indeed, if we consider the 
mean-field Hamiltonian 

t u 

n^hip+-ip' + -ip\ (5.26) 



the stationarity condition gives 



h + t^ + -^^ = 0, (5.27) 

which is solved by (/? = Ah^^^ f {Bt\h\~'^^^) , where f(x) satisfies eq. (5.25), and A and B are 
constants depending on u. This identification also shows that M plays the role of an external 
field, while m ~ is the magnetization. 



5.3 Vertex, Propagator and the failure of expan- 
sion 

In order to perform the l/Nf calculation, we expand the field (p around the saddle-point 
solution, 

0(x) = + ^0(x), (5.28) 

gvN 

where N = DNf, and 0(x) in Fourier modes: 

0(x,, x,+0 ^tJ2 e-"'^-^^^^^ I 0„(p)e^P-<^ . (5.29) 

In the following we will refer to the integers n — or more precisely to 27toT — as frequencies. 
In this way we obtain the following expansion for the effective action: 



eff 



5.3 Vertex, Propagator and the failure ofl/N expansion 
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T 



p<a(2^ 



T 



i-i 



l>3 {m} '^P^<'^ ^ ' ^ ^ j=i 

^ni (Pi ) • • • $ni (Pi ) V^^^ {pi,ni]---]pi,ni). (5.30) 



Note the Kronecker 6 on the frequencies that ensures that vertices are nonvanishing only for 
Yli'>^i — 0- III particular — this property will be important in order to obtain the effective 
action for the zero mode — if only one frequency is nonvanishing y(')(pi, n; p2, 0; • • • ; p;, 0) = 
0. 

The vertices appearing in expansion (5.30) are easily computed. The fermion contribution 
is obtained by considering the one- loop fermionic graphs in theory (5.1). If we define the 
free fermion propagator 



AF{p,n;m) 



with ujn = {2n -\- l)7rT, the fermion contribution is 
y/^(pi,ni;---;P/,n/) = 

+ permutations, 



(5.31) 



■tr 



,i=l j=l j=l 



(5.32) 



where the permutations make the vertex completely symmetric [there arc (/ — 1)! terms]. 
Note that the frequencies a;„ never vanish and thus the vertices have a regular expansion in 
powers of m + M. Vertices V^''^ satisfy an important symmetry relation. First, note that 



7d+iAi7.(p, n; m)7d+i = -Af{p, -n - 1; -m). 



It follows 



E 

aez 



tr 



JJ AF(qi;a + 6i;m) 



i=l 



E 



tr 



JJ 7d+iAF(qi; a + m)'yd+i 



i=l 



(5.33) 



(5.34) 



n ^Ficii] -a -hi- 1; -m) 



,j=i 



= (-!)' E 



tr 



U Ai.(qj;a - h^] -m) 



.1=1 



In the second step we used 7^^^ = 1, while in the last one we redefined a — > —a + 1. This 
relation implies (we write here explicitly the dependence of the vertices on m and M)^ 



y/^(pi, m; . . . ; Pi, n/; m + M) = (-1) V/^(pi, -m; . . . ; p^, -m; -m - M) 



(5.35) 



^If d is odd, one can repeat the same argument using 7c defined in (5.3). It shows that vertices with I 
legs are multiphed by (—1)' if one changes the sign of m + M at fixed momenta and frequencies. 
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For every / > 4, the vertex is due only to the fermion loops, so that V^'^ = Vj''\ For / < 4 
we must also take into account the contribution of the bosonic part of the action, so that 

V^^\p,ni;q,n2;r,n3) = Am + V}^^^(p, m; q, ^2; r, na), 

y("^^(p,ni;q, 712; r,n3; 8,^4) = A + V}^''^(p, ni; q, ^2; r, 713; s, 714). (5.36) 

Finally, for the inverse propagator we have 

P{p,n) = ^+^^^+lI+^m' + v}'\-p,-n;p,n). (5.37) 

Vertices V^*^'-* also satisfy the symmetry relation (5.35), while the inverse propagator P(p,n) 
satisfies P(p, n; m, M) = P(p, —n; —m, —M). In order to clarify the equation of the critical 
crossover limit we need to explicitly compute P(p, n), V3(p, n) = T^^^^(p, n; — p, — n; 0, 0) and 
V3{p,n) = y^^^(p, 7^; — p, — n; 0, 0; 0, 0) for M, ttt. — > 0. Using standard formula Tr7^7i, = 
DS^^,,, Tr 7^7^75 = and Tr -f^-fp-f^-fg = D {Sa,f3S^,s - Sa,jSf)^s + Sa,sSf3,j) and defining 

Th{p,n) = 2vh{p,n) +Vh(O,0) 

-T(p2 + (27rTn)2)^ J ^ r;^(q, a) V(q + P, « + ^) 



aez 



Vh{p,n) = r]h{(l,a)r]h{q, + p,a + n) (5.38) 

with 77^(p, n) := (p2 + 7r^T\2n + if + /i^) ^ , we get the final result 

9 9 2 

+ (y + ^^)^.(p, n) + 2/i2^;,(0, 0) (5.39) 

V3{p,n) = Am + 2(m + M)r^+M(p,n) (5.40) 

-8{m + MfTY, / WM(q,a)V+M(q + P,a + n) (5.41) 

y4(p,n) = A + 2T^+M(p,n) + 0(m + M)' (5.42) 

The similarity of (5.41) with (5.42) simplifies the equations of the critical crossover limit 
with respect to the general case. This is due to the symmetry of the model (5.35) that is not 
present in other models investigated in this work. Indeed in the general case one has to tune 
the model at the critical point recovering the symmetry which is than dynamically broken. 
However in this case the critical point remain fixed at M = without 1/Nf fluctuations, 
technically this follows from (5.41) and (5.42). 

Finally, note that V'^^^ (0, 0; 0, 0; 0, 0; 0, 0) is positive at the transition. Indeed, one obtains 
exphcitly (note that A > to ensure the stability of the quartic potential) 

V<^) (0, 0; 0, 0; 0. 0; 0. 0) = A + 6T. E £^ 1^ > «. (5.43) 
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It is easy to verify that P{0, 0) vanishes at the transition. Indeed, for m = M = we have 

= -goi{T-T,) + 0{T-T,y, (5.44) 

where we used eqs. (5.20) and (5.17). Thus, the mode with n = is singular at the critical 
point. It is exactly this singularity that forbids a standard l/Nf expansion at T = Tc and 
gives rise to the Ising behavior. This type of singular behavior is the main topic of this work 
and has been described in the first chapter of this works in a general way using generalized 
Heisenberg Models. The strategy proposed there consists in integrating all the nonsingular 
modes 0„, n ^ 0, and study the effective theory for the zero mode 0o ^-nd is the same we are 
going to apply for the model under investigation. 



5.4 The effective theory of the zero mode 

Following the conclusion of the previous section and the standard protocol developed in the 
first chapter, we compute now the cfi:ective action for the zero mode integrating the non 
critical degrees of freedom. Integrating all fields (pn with n 7^ we obtain the effective action 



g-5eff [^0] 



f J]d(^„e-^««[^"], (5.45) 



with 

4fr = ^fNHUO) + ^l 7^0o(p)P(p)0o(-p) (5.46) 

+ g^^i<,(2^---(2^^2^)HgN 

xy(')(pi,...p,)0o(pi)---0o(pO 

where H, P, and V^''^ have an expansion in powers of 1/N. The computation of these 
quantities is quite simple. The contribution of order 1/N'' to V^''^ is obtained by considering 
all A;-loop diagrams contributing to the Z-point connected correlation function of 00 and 
considering only the nonsingular fields (i.e. propagators with n 7^ 0) on the internal lines. 
Frequency conservation y(")(p, n; 0, 0; ■ ■ ■ 0, 0) ~ Snfl implies that all tree-level diagrams 
with more than one vertex vanish.^ Therefore, H = 0{N~^), P(p) = P(p,0) + 0{N^'^), 
and V^*^''' = V^^'Lo + 0{N~^) (V^^'Lo is the vertex 1/*^'^ with all frequencies set to zero). For the 



^For a trco-lcvcl graph, the usual topological arguments give the relation X]n(^ ~ 2)iV„ = —2, where 
is the number of vertices belonging to the graph such that n legs belong to internal lines. Since n > 1 if 
there is more than one vertex, the previous equality requires A'^i > 2. But frequency conservation implies 
that Vi vanishes if all frequencies but one vanish. Therefore, each nontrivial tree-level diagram vanishes. 
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inverse propagator P and for the magnetic field H we shall also need the 1/N corrections. 
We obtain 

Hi 



H 
PiP) 



^ +0{N-') 



with 



Hi 

p(^)(o) 



T 



-E 

2 ^ 



d'^P 



1-1 



d'^p 



a ^3(p,n)P(p,n)" 

V4(p, n)P(p, n)-' - Vsip, n)2p(p, n) 



(5.47) 
(5.48) 

(5.49) 
(5.50) 



where Vi{p, n) = ^'^'•'(p, n; — p, — n; 0, 0; . . .). Note that relation (5.35) implies an analogous 
symmetry relation for 



V^«(pi, ...,pi;m,M) = (-1)'V^« (pi, . . . , p,; -m, -M), 



(5.51) 



where we have written explicitly the dependence on m and M. Analogously P(p) and H 
are respectively symmetric and antisymmetric under m, M ^ — m, — M. 

In order to obtain the final effective theory we introduce a new field x(p) such that the 
corresponding zero-momentum three-leg vertex vanishes for any value of the parameters. For 
this purpose we write 

axip) = r0o(p) + VNk5{p), (5.52) 

where a and k are functions to be determined. If we write ai = y*^'^(0, 0; . . . ; 0, 0), k is 
determined by the equation 



^(-iyk(m,My 
Yl M «/+3(m, M) = 0, 



(5.53) 



where we have written explicitly the dependence on m and M. Now, symmetry (5.51) implies 
also 



^{-iy[-k{-m,-M)y . ^ 

Z2 — —ai+3{m, M) = 0, 



1=0 



l\ 



so that k{m, M) — —k{—m, —M). Therefore, k has an expansion of the form 

k^ Y Kbm''M\ 



(5.54) 



(5.55) 



a6,a+6odd 



where the coefficients kab have a regular expansion in powers of 1/N. The leading behavior 
close to the transition is easily computed: 



03 



A;= — + 0(m"M^a + 6 = 3). 
04 



(5.56) 



^For a diagram D{m, M), with nj j— legs vertices, entering in the computation of V'-'''' wc have D{m. M) = 
{-l)^j"^W{-m,-M) = (-l)'£>(-m, -M). We have used the relation Y.j'^jj = I + 2n/, m being the 
internal line 
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In terms of x the effective action can be written as 

id 

Arl/2T7_ , / >^ 

>efr 



iV'/2-ffx(0) + i / j|^x(p)/'(p)x(-p) (5.57) 

- g • ■ ■ E P.) ^'"(P.. ■ ■ ■ P.)x(P.) ■ ■ ■ X(P,). 

Tlie quantities H, P, and V^''^ have an expansion in terms of m, M, and Explicitly we 
have: 

H = aT-^[H - kP{0) + —VaiO) - —V^O) + 0{rrfM\ a + 6 = 5)], (5.58) 

Zi o 

P(p) = a^T-\P{v) - + y K4(p) + 0{nfM\ a + 6 = 4)], (5.59) 

V^''^'\^u . . . ,P2.+i) = a^'+^T-^[V^(^^+^)(p„ . . . ,p,,+0 - A;V^(^'+^)(p„ . . . ,p.,+i,0) 

+0(m"M^a + 6 = 3)], (5.60) 

F^'^^pi, . . . , P2O = a2'T-^\/(2')(pi, . . . , P2O + 0{rrfM\ a + 6 = 2). (5.61) 

Up to now we have not defined the parameter a. We will fix it by requiring 

dP(p 



= 1, (5.62) 

p=0 



for all values of the parameters. The parameter a is a function of m, M, and The 
symmetry properties of k and of the vertices imply that a is invariant under m,M ^ 
—m, —M. As a consequence, under m,M^ —m, —M, the quantities H, P, and V^^^ have 
the same symmetry properties as P, and V^^\ 

In the following we shall need the expansions of if, P(0), and V^*^^''(p, — p, 0) close to the 
critical point. Using (5.39), (5.41) and (5.42) we get 

P(0, 0) ^ ^m^ - (T - T,)(?oi - 3(M + m)'<7io, 
1^3(0, 0) a; Am - 6(M + m)^io, 

14(0, 0) «i A - 6^10, (5.63) 

Notice the relation 

Vf'\p, n; -p, -n; 0, 0; m) = my/^^(p, n; -p, -n; 0, 0; 0, 0; m) + O(m^), (5.64) 

(we have explicitly written the mass dependence of the vertices) that we have just commented 
in the previous section. In the relations (5.63) we have used gio = —Vo{0, 0) that can easily 
obtained using the development of the gap equation near the critical point. We expand H 
and P(0) is powers oi 1/N as 

^ = K + ^ + 0{N-% (5.65) 
a iV 

^=Po + ^ + 0(iV-2). (5_gg) 
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By using expansions (5.63) we obtain 

^ -m + -^^M{T - TJ - ^ioA(6^io + A) ^3 ^ ^^^a^"^ a + 6 = 5), 
e^fio - A (6^10 - A) 2 

^ / 7^ [6^10 + (P, -p, -n; 0, 0; 0, 0)]P(p, n)'' 

Jp<A K'^'^) 



2(6c/io - A) ^qJp<a 

+0{m''M\a + b = 3), 
^3(0, 0)2 



Po ~ ^(0,0)- 



214(0,0) 



^ -goiiT - Tc) + ^^'"^- M^ + 0{rrfM\ a + 6 = 4), 

65(10 - A 

j)i ^ P(^)(0,0) ^ e + 0(m"M\a + 6 = 2), (5.67) 

where e is the value of P^^\Q^ 0) for M = m = 0. Note that several terms that are allowed by 
the symmetry m,M —m, —M are missing in these expansions. In the case of we used 
the gap equation to eliminate the term proportional to m^. This substitution is responsible 
for the appearance of the term linear in M and cancels the terms proportional to m(T — Tc), 
m^M, and mM'^. In the case of hi and po note that the terms proportional to m, and m^, 
mM cancel out. Finally, we compute the three- leg vertex. At leading order in 1/N we obtain 

|f-(p,-p,o)..3(p,o)-|M^,p.o) 

[65lo + V}^^Hp,0;-p,0;0,0;0,0)] + 0(m"M^a + 6 = 3). (5.68) 



65^10 - A 

Note that the term proportional to m is missing as a consequence of relation (5.64). 

Having computed the expansion of translated vertices of the theory with one incoming 
momenta (5.67,5.68) and at zero external momenta (5.63) we are now in a position to define 
the critical crossover limit for the model under investigation. 



5.5 The critical crossover limit 



The manipulations presented in the previous section allowed us to compute the effective 
action for the zero mode x(p). Far from the critical point -P(p) 7^ for all momenta 
and thus one can perform a standard 1/Nf expansion. At the critical point instead this 
expansion fails because -P(O) = 0. At the critical point, for N ^ 00 the long-distance 
behavior is controlled by the action 



5eff 



/ 



d'^x 



+ 0{N-''), 



(5.69) 
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where 

^^=^^^'^0, 0,0,0). (5.70) 

Here we have used the fact that vertices with an odd number of fields vanish at the critical 
point and the normalisation condition (5.62). Moreover, since the critical mode corresponds 
to p = 0, we have performed an expansion in powers of the momenta, keeping only the 
leading term. Since for N ^ oo, F^^\o, 0, 0, 0) = a^T'^y ^(0, 0, 0, 0), inequality (5.43) 
imphes u > at the critical point, eq. (5.69) is the action of the critical 0^ theory which 
should be studied in the weak-coupling limit u ^ 0. In this regime the model shows an 
interesting scaling behavior — we named it critical crossover — that describes the crossover 
between Mean-Field and Ising behavior and has been studied in detail in the chapter 2 
and applied for a specific model in chapter 3. In this chapter we want to recover the results 
obtained for Heisenberg models in chapter 3 for Yukawa model (5.1) using the general M —1 
theory presented in chapter 2. 



5.5.1 The general theory 

Following Ref. [56] in chapter 2 and 3 we extended these considerations to the general two- 
dimensional Hamiltonian 



d^p 



[i^(p)+r](^(p)(^(-p) 



u 



l\ 



(5.71) 

(/9(pi)...(/9(p,), 



Z>3 

where K{j>) =p'^ + 0{p*), V^^^ (0,0,0) = 0, and V(^)(0, 0, 0, 0) = 1. The presence of vertices 
with an odd number of legs requires an additional counterterm for the magnetic field. Indeed, 
we showed that it was possible to find functions rc{u) and Hc{u) such that for i = r — rc{u) 
(infrared limit), h = H — Hc{u), u ^ (weak-coupling limit), at fixed t/u, h/u one has 



(5.72) 



where the scaling function fn{x,y) is the same as that computed in the continuum theory. 
In particular, Xn^"^ vanishes in the crossover limit if n is odd. The counterterms are 
regularization-dependent. In the continuum theory with a cutoff we have 



hr — 



2 (27r)2 i^(p) 



u ^ u u 
— In — + - 



d^P V3(p) 



Stt A2 2 J^^^ (27r)2 K{pY 
with V3(p) = V(^)(p, -p, 0) and 



+ AqU, 



(5.73) 
(5.74) 



3 1,31 
-Do \ log 



d^p 

(2^ 



VW(p,-p,0,0) 
K{p) 



(5.75) 
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The nonpcrturbative constant D2 was estimated in Ref. [32]: D2 = —0.0524(2). The relations 
(5.73) and (5.74) generahze the lattice expressions found in chapter 2 [see eqs. (2.47) (2.75)] 
to a sharp cut off regularisation. However as will be clear later in this case (in which chiral 
symmetry is not broken in the starting action) he is irrelevant in the scaling limit. 



5.5.2 Scaling behavior 



In this section we wish to use the previous results to compute the crossover behavior of 
model (5.1) in 2+1 dimensions. Since u 1/N the relevant scaling variables are 



Xt 



{N^/^H - He) 



a 
NT 

9^(P(0)-re), 



(5.76) 



where the factors a/T^ and a'^/T^ are introduced for convenience. The critical crossover limit 
is obtained by tuning T, M, and close to the critical point so that and Xt are kept 
constant. The expansions of H and P{0) are reported in eq. (5.67). The expansions of H^. 
and Tc are easily derived. For He we have 



H. 



1 



2VArjp<A(27r)2 P(p) ' 

where all quantities are computed for M = m = 0. Using eq. (5.68) we have 

ahcoM 



(5.77) 



Hr 



+ 0(m"M^a + 6 = 3), 



TeVN 

where hco is a constant. Using eq. (5.74) we obtain for Tc the expansion 



(5.78) 



a 



TeN 



{rolnN + ri) + 0{m''M\a + b = 2), 



where 



^0 = 



1 j^3aV4(0,0) 

8^ ^ 87rT,A2 



D2 

Stt 



d=^p 



TeV4{p,0) aV4(0,0) 

P(P,0) p2 



(5.79) 



(5.80) 



(5.81) 



Note that the three-leg vertex that appears in eq. (5.74) does not contribute to this order, 
since it vanishes for m = M = 0. Thus, eqs. (5.76) can be written as 



Xh 
Xt 



7V3/2 

N 



-JIM + aoM{T - Te) + a^M^ + ^ + 



-goi{T - Te) + agM^ + ^ + . . .] _ In TV - n. 



heoMVN, (5.82) 
(5.83) 
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where a^, ai, a2, 03, and e are coefficients that can be read from eq. (5.67). These expansions 
show that 

-pM = -XhN-^/^, (5.84) 

Ngoi ^ ' 

The critical point is specified by the condition Xt — Xh = 0. The symmetry under m, M 
—m, — M guarantees that the critical point corresponds to M = 0. On the other hand, 
fluctuations give rise to a shift of the critical temperature. If Tc{N) is the finite- critical 
temperature, we obtain 

Tc{N)^T,+ ^ ' - -'\ (5.86) 

901 

Note that, beside the expected 1/N correction there is also a,\aN/N term that is related to 
the nontrivial renormalization. It follows 

r-T,(iV) = --^. (5.87) 

Note that Qqi is negative [see eq. (5.21)] and thus we have > for T > Tc{N), as expected. 
Using the gap equation we can also derive the behavior of m in the critical crossover limit. 
We obtain 

where mo is a function of In N that satisfies the equation 

^(A - 6gio)ml + {rolnN + n - e + Xt)mo + - 0. (5.89) 

For — > 00, mo has an expansion in inverse powers of \nN, the leading term being 

mo - — + 0(ln-2 N). (5.90) 
ro m JM 

Note that mo as Xh ^ 0. 

These results confirm the scaling predictions of Rcf. [50]. For the massless theory with 
M = 0, there arc two regimes: for N{T — Tc{N)) -C 1 one observes Ising behavior, while 
for N{T — Tc{N)) » 1 mean-field behavior occurs. If M 7^ the same considerations 
apply, the relevant variable being MN^^^. It is important to note the role played in the 
derivation by the symmetry m,M —>■ — m, —M, that is present because the regularization 
preserves chiral invariance. Even though vertices with an odd number of legs are present, the 
symmetry makes them irrelevant in the crossover limit. Thus, the additional renormalizations 
computed in Ref. [56] do not play any role here. 

The results reported above can be extended to d dimensions for d < 4, the relevant 
scahng variables being eqs. (2.52) 



^^(d+2)/[2(4-d)] ^ Xt^[T- T,(iV)]iV2/(4-'^) . (5.91) 
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In c? = 3, on the basis of eq. (2.55), we also predict for Tc{N) an expansion of the form 

where a, b, and c are constants that can be computed as in the two-dimensional case. However 
notice that, as pointed out in sec. 2.1.5, in three dimension (and in general for d > 2) one 
have to translate the critical field by a constant (2.57) that have to be taken into account 
in order to cancel some diagram for X3 that do not go to zero in the Critical Crossover 
Limit eq. (2.61). However in the symmetric case this is not necessary because the explicit Z2 
symmetry of the action is protected in the scaling limit. Indeed the two one loop anomalous 
contributions (= £>f-*,I?2 ) Xs [see point c) in 2.1.5] are-'^^ 

D2 = ^N^J^P (5-94) 

In the critical crossover limit in d = 3 from eq. (2.54) follows that X3 ~ A^^/^ and that 
t - so that 

Ar-9/2Df) ^ 7V-9/2L»f ^ M (5.95) 

that are obviously irrelevant in the CCL. We stress that the reason Dr and D'^ ' are negli- 
gible is the equation (5.64) that holds for all models which arc chiral symmetric. In general 
the effect of ku must be taken into account for all Hamiltonian which broke chiral symmetry 
as for instance in the case of Wilson fermions. This system is currently under investigation. 

5.5.3 Correlation functions 

The results reported in sec. 5.5.2 allow us to compute the scaling behavior of the correlation 
functions. For instance, we have 

(0(x<i, Xd+i)) = - - + ^=(x(xd)) (5.96) 
9 gv N 

Using eq. (2.54) with n = 1 and d = 2,we have (x(xd)) = fi{xh,Xt) in the critical crossover 
limit. The background term can be neglected in the crossover limit since 



pa - 

9 9 

Thus, we obtain 



_mo) 

V^4(0,0) 



(5.97) 

9 X- 6gio 



(X 

{<p{y.d,Xd+i)) ^ — j=h{xh,xt) (5.9^ 

9yN 



10 



In three dimension the integrals are both Infra Red finite. 
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The factor 1/yN is related to the particular normalization of used in (5.1) and disappears if 
we redefine (p = \fN(j) in order to have a canonical kinetic term for (p. The function fi{xh, Xt) 
is the scaling function for the magnetization in the Ising model. For instance, for Xh — 
and Xt < (low-temperature phase), we have /i(0, Xt) ~ {—Xt)^' and fi{0,xt) ~ (— x^)^^^ 
respectively for \xt\ ^ 1 and \xt\ ^ 1, where f3j = 1/8 and f3MF = 1/2 ai'e the magnetization 
exponents in the Ising and in the Gaussian model. The universality of the crossover allows us 
to compute the scaling functions in any other model in which there exists a crossover between 
the Gaussian and the Ising fixed point. In particular, we can use the results for systems with 
medium-range interactions [33, 36] (see also the appendix B). In Ref. [36] (LBB) the authors 
report {\m\R) versus tR'^ (see their Fig. 9), where m is the magnetization, t the reduced 
temperature, and R the effective interaction range. These results give us {(f){'Xd,Xd+i)) for 
Xh = 0. One only needs to take into account the different normalizations of the fields, of the 
couphng constant, and of the scahng variable. In the crossover limit A?" — > oo, T ^ Tc{N) 
at fixed A^(r - Tc{N)) we have 

g^/N {(l){^d, Xd+l)) = Kl,LBB(|m|i?)LBB (5.99) 

{tR\BB = Ki^bbN[T - Te(7V)] . (5. 100) 

The nonuniversal constants i^^LBB and i^i,LBB are computed in the app. B. 
It is customary to define an effective exponent Pes{T) as 

/?eff(r) = [r-r,(iv)]Ain(0(x,,x,+i)), (5.101) 

for M = and T < Tc{N). In the crossover limit T Tc{N), N ^ oo at fixed N[T-Tc{N)], 
the exponent Pes{T) interpolates between the Ising value Pi — 1/8 and the mean-field 
Pmf — 1/2. Again, this effective exponent can be derived from the results of Ref. [36]. The 
curve reported in Fig. 15 of Ref. [36] gives Pes in the Yukawa model once tR^ is replaced by 

i^LBB[T-Te(7V)]iV. 

The same considerations apply to the connected zero-momentum n-point function Xn- 

Xn ^ J d'^+^X^ . . . d'^+^X^ (0(O)0(X2) . . . 0(XJ)-- 

- J ^'^^ " " " ^'-n (x(O)x(x.) . . . x(x„)) — 



^4(0,0) 'V -/^-Vn(x.,x,), (5.102) 

gn 

For n = 2 the crossover function for Xh = can be obtained from the results of Ref. [36], 
since g^X2 — -f^2,LBB(x-R^)LBB- The constant i^^2,LBB is given in the Appendix. 

One can also use field theory to compute the crossover curves and thus use the results of 
Ref. [32]. For instance, in the high-temperature phase, for M = we have in the crossover 
limit 

9^2 = K2,FTF^{i), i = KftN[T - niN)], (5.103) 

where F^{i) is reported in Ref. [32] and Kp-£, K2^ft are nonuniversal constants computed in 
the appendix B. 
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In the discussion presented above we have focused on the case d = 2, but it is immediate 
to generahze all these considerations to the three-dimensional case. For d = 3 the universal 
crossover curves have been computed in Refs. [29, 30, 32, 31] (field theory) and in Ref. [37] 
(medium-range models). These results apply directly to the Yukawa model. 

5.6 Conclusions 

In this chapter we have considered the Yukawa model in the limit Nf oo, focusing on 
the crossover between mean-field and Ising behavior. For this purpose we have determined 
the action of the mode that becomes critical at the transition. In the long-distance limit, it 
becomes equivalent to that of a weakly coupled 0^ theory. This identification allows us to use 
the results available for this model summarized in chapter 2 ([29] [56]). In particular, we have 
identified a universal critical crossover occurring for Nf oo, M — 0, and T — Tc{N) 
at fixed Xt and Xh, see eqs. (5.84), (5.87), (5.91). In field-theoretical terms, this behavior 
represents the crossover induced by the flow from the unstable Gaussian fixed point to the 
stable Ising fixed point. Quantitative results for the Yukawa model can be obtained by using 
the field-theoretical results of Refs. [29], [32], [31], or the Monte Carlo results available for 
medium-range models [36, 37]. The necessary nonuniversal renormalization constants can 
be computed in perturbation theory. Results for = 2 arc reported in the appendix. 

We should stress that our results are not specific of the chosen regularization, but can 
be extended to other regularization as well. In particular, the extension to Kogut-Susskind 
fermions [61], the model considered in Ref. [50], is essentially straightforward and have not 
been presented here. The Wilson case is more involved. Indeed, the absence of chiral 
symmetry implies that the symmetry relations satisfied here by the effective vertices [see 
eq. (5.35)] are no longer valid. In turn, this may imply additional mixing as it happens in the 
generalized Heisenberg model [56]. Also the higher than two dimension is more complicate 
for Wilson fermions, due to the necessity of introducing and tune another parameter kR. 

Let us note that all calculations presented here refer to the model in infinite spatial 
volume. However, the crossover behavior can also be observed in the finite-size scaling 
limit. The discussion in sec. 5.5.1 can be easily extended to this case too. It is trivial to 
verify that the correct scaling variable is L = Lu^/^^~^\ i.e. L — LN~^I^^~'^^ in the Yukawa 
model. Again, one can use universality and obtain predictions for the Yukawa model from 
the results obtained in other contexts. In particular, one can use the finite-size scaling results 
of Refs. [33, 36, 37] that refer to medium-range models at the critical point. 

Finally, we should mention that one could also generalize the model and consider fermion 
fields il^af transforming according to a representation of a group G and a coupling of the form 
fT°'(f)°"^f, where are the generators of the algebra of G. The discussion is essentially 
unchanged, though in this case one would obtain the vector 0^ theory; several models have 
been just studied pointing out the reduction of the critical zone [51] [52] [53]. Field-theory 
results relevant for this case are given in Refs. [29], [32], [31]. More in general we believe that 
the considerations and the results obtained in this chapter are not peculiar of fermion models 
but could be recovered for instance also for vectorial models at finite temperature. -'^^ This is 
not an unexpected result in the light of chapter 6 where the large N scheme studied in this 



^^See for instance [54]. 
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work is applied to a vectorial 0^ interaction. Indeed using standard arguments one can see 
as the introduction of a cut off (the temperature in the (j)^ interaction) can be resummed in 
higher order vertices. 

We have just stressed that a careful numerical investigation of fermionic models (of the 

type presented in this chapter) could be a very precise test of the idea presented in this work, 
due to the fact that a lot of available simulations (also for several Nf) are present [50] [51] 
[52] [53]. The goals of such a studies would be the determination of the correction to 
the critical parameter Tc{N) [and Mc{N) in the case of Wilson fermions] to compare with 
the simulation data for several N. This step could be require few time and is actually under 
study. On the other hand is not clear to us if the data available in literature could confirm 
the claim that all the crossover considered in this work are universal, i.e. if for instance 
medium range models can be compared with Yukawa models. It is our opinion that such a 
question merit further investigations. 



Chapter 6 

Crossover in tricritical models 



In this chapter we want to give another example in which our technique can be used to 
elucidate some aspects of a tricritical phase diagram in the large limit. The Hamiltonian 
(6.1) is usually used in order to investigate systems that exhibit a tricritical phase transition 
like for instance metamagnet [97] or fluid ^He-'^He mixtures. The large N limit has been 
investigated long time ago [96]. The standard scenario predicts (for H ^ Qi) two hnes of 
continuous phase transitions (that usually are called wing lines) that end joining together in 
a tricritical point. However some questions remain open. Indeed the ref. [96] shows that the 
continuous transitions observed along the wing lines are Mean Field like while because of 
symmetry argumentation, the authors expect a Ising behavior for every N . In this chapter 
we want to show as the two different critical behavior can be reconciled defining the Critical 
Crossover Limit for Hamiltonian (6.1). The presentation strictly follows chapter 1 and 3 and 
some steps are missed. In particular we take care to define (sec. 6.2) the N = oo scaling 
fields of the theory (with which using results of sec. 1.3 one can recover the N = oo phase 
diagram) and to study the stability of the weakly coupled theory (sec. 6.4) that is the basic 
point in order to define the Critical Crossover Limit which implementation for interaction 
(6.1) is identical to what done in chapter 3 for jV = 1 Heisenberg models. 

6.1 The model 

In this chapter we will consider a vectorial model with a (j)^ interaction. Let us discuss the 
following lattice regularization 

(6.1) 

— * — * 

where H and are two vectors with components equal respectively to H and H^. 
Hamiltonian (6.1) has been used to investigate tricritical point. In particular the common 
accepted phase diagram exhibits two lines of continuous phase transition (that are usually 
called wing lines) that end in a tricritical point. Detecting the tricritical point in general is 
not an easy work and it requires an accurate fine tuning of r, m and v (while for symmetric 
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reasons at the tricritical point one has H = = Oy (for more details we refer to [98]). 

The large limit can be studied introducing auxiliary fields as done in chapter 1. In 
this case we need only two auxiliary fields (A, /x), because we relax the unit spin constraint 
in order to follow [96]^. Thus we can write 



Z = 



-«[</>] 



Y[ d^xdUxdXx exp ^ ^ (px$x+ix + ( y ('"^ ~ ^~^) "^^^ {H3$x)l^x 

X X,^ X ^ 



X 

Nr 



Nu ^ Nv ^ 



(6.2) 



Integrating the (p field in the previous expression we obtain the effective action A for the 
auxiliary fields that is the starting point for every 1/N investigation: 

y tr log 0[X] + ^{H + Hspi, 0-'[XlH + Hspt) 



Y\diixdXx exp 

X 

^ fN Nr Nu 2 Nv s\ 

+ ^ I -^Xx^ix — 2~^=^ — iT^"" — eT^"^ I 

X ^ ■ ■ / 



(6.3) 



where we have introduced the operator O that acts on scalar fields in the following manner: 



{ip,0[X],ip) = -^ipxipx+ix + ^X, 



<Px 



(6.4) 



XyfJ, 



For constant value of = A, O can be diagonalized in momenta space so that 

^[A](P) = ^(p^ + 2A-2d). (6.5) 

6.2 Gap-Equation. The = case 

We are interested to understand the crossover from Mean Field criticality to Ising criticality 
on the wing lines of the tricritical phase diagram. In order to do that it it is expected [96] 
we can put if 3 = without destroy the qualitative picture. Indeed imposing the stationarity 
of the effective action (6.3) we obtain two equations for the fields at the stationary point JI 
and A 



Jl = 2Bi(ml)+AH^B2(ml) 
X 



U _ V _n 

'+3!^+5!^ 



(6.6) 



^However the study of interaction (6.1) in the weak coupUng Hmit (chap. 2 for A/" = 3), could give the 
leading term of the expansion of rdv) and Uc{v) apart a ~ t;-factor similar to A defined in sec. 2.1.4. 

^However as we have just noticed during this work, in the large N limit unit spin constraint is usually 
irrelevant [24] [49]. 
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Figure 6.1: H'^ versus with v — 1/4, u — —1/2 and (from below) r — 5, 5.05, 5.07 and 5.15 



where we have defined 



m. 



2 _ 



= 2X-2d. 



In eq. (6.6) we have taken the notation of the previous chapters defining 

( (p2 + mg)- 



Bn{ml) 



(6.7) 



(6.8) 



From the previous equation (6.6) we can obtain H (and of course A) as a function of the 
0{N) symmetric observable Jl:^ 



^2 ^ /i-2Ei(mg) 



(6.9) 



where function of /i is obtained using eq. (6.7) and the second of eqs. (6.6). Interesting 

enough is the fact that as a function of JI in general cannot be inverted as depicted in fig. 
6.1. This is the sign of the appearance of a phase transition, the mechanism being the same 
of what presented in chapter 1 or ref. [55] for Heisenberg models. In this case the magnetic 
field H plays the role of the temperature (or spin-spin correlation length) in Heisenberg 
models. The key point is that in both model phase transitions happen at finite temperature 
(/? < cxo for the Heisenberg models) and at finite Magnetic Field {H ^ for the model 
studied in this chapter). This fact is the reason for which infrared divergences appear and 
the standard expansion fails, so that the crossover mechanism is observed instead of a one 
parameter (A^) family of critical points. 

Figure 6.1 gives evidence of the existence of a line of critical point (wing line) near which 
we have (in the following we reduce to the H > case): 



(6.10) 



■^Indeed using the equations of motions we get /x =< \(j)x\'^/N > 
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At this point the Mean Field considerations of the chapter 1 and of reference [55] can easily 
be applied with analogous results. In order to elucidate this claim, if we take as tunable 
parameters of our theory r and H {u and v remaining fixed for instance at the value reported 
in fig. 6.1) and expanding the gap-equation near the critical point He ~ 23, Tc ~ 5.07 and 
77c ~ 14 (see fig. 6.1) we get the following expansion 

H'-H', = ^a„4r-r,)"(71-7lj- 

n,m 

aoi = fori = 0,1, 2. (6.11) 

Equation (6.11) and the gap equation ior J\f — 1 Heisenberg models (1.21) are identical apart 
replacing r with p and Ji with ttIq. In particular proceeding in such a way the N = oo scaling 
fields can be obtained using eq. (1.24) and eq. (1.28). Then following the considerations 
of chapter 3 we expect to extend the N = oo analysis including also 1/N fluctuations. In 
particular this allows us to explain the crossover between Mean-Field to Ising behavior using 
universal crossover functions that have been investigated in chapter 2 and have been used in 
all the rest of this work. 

In the next section 6.3 we will obtain the Propagator and the vertices of the auxiliary 
fields 11, A, in particular we want to show the presence of the zero mode of the theory. Having 
done this it will be clear that the 1/N results obtained for the Heisenberg models (chap. 3) 
or for the Yukawa model (chap. 5) strictly apply also for Hamiltonian (6.1), confirming the 
claim that for every N finite on the wing lines one observes Ising phase transitions that cross 
-in a universal way- on Mean field phase transitions for = oo. 



6.3 Expansion and Crossover to Ising Behavior 

Now we are going to show that eq. (6.10) implies the existence of a zero mode, so that the 
picture developed in chapter 3 can be applied also in this case. Introducing fiuctuations for 
the auxiliary fields to the saddle point solutions (6.6) 

fix 



N 

Xx = A + 4^ (6.12) 
V N 



the standard 1/A^ expansion [characterized by the inverse of the propagator PAB-jiiP) 

by the effective vertices V^{^..^ .p(pi, • • -pj)] is easily obtained. In order to give an outlook 
on how a zero mode appears near a gap equation that exhibits criticality like that presented 
in fig. 6.1 we try to differentiate the gap-equation (6.3) (rewritten below) 

^J-^^A[^(jj),H{jl)]=0 (6.13) 

[^'i(p) = /i(p) and ^2(p) = '^(p)] with respect to Jl obtaining 

^ 6'S[^i-p), Hj-p)] d^B H{ji)] dH _ 

^ (5^^(0)(5*b(0) dp b'^AWa dp ' ^ ' ' 
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In the previous equation (6.14) at the critical point the second term cancels [like (/J — JIc)^] 
pointing out that is singular at the critical point confirming the starting claim. In 
particular (denoting with ^a(p) the zero mode of the theory), neglecting a normalization 
factor, from eq. (6.14) easily follows 

E^Ai(0)^B(0) ~ ijl-JI^fzAiO) 
B 

zb(0) = ^ + 0{-p-Jl^r. (6.15) 

The previous equation is similar to relation (3.29) found for Heisenberg models.^ 

To explicitly verify the previous consideration we compute detP~^. It is easy to obtain 
the following formula 



2u vfj, 1 dm, 



2 



4! 5! 4 dp 
1 

2 



P,-,^(0) = -8H^B,{ml)-2B,{ml)=(^^^ '(^-^^^^°)^) ^^"^^^ 
Using (6.9) we obtain 

detP-^(O) = -B2{ml)^ (6.17) 

confirming the presence of the zero mode at the critical point (fig. 6.1). Notice that the sign 
of dH'^/dJI is positive near the critical point [i.e. Ooa > in eq. (6.11)], this will be relevant 
to discuss the stability of the theory. 



6.4 On the stability of the Weakly-Coupled theory 

In this chapter we want to do some considerations on the stability of the weakly coupled 
theory describing the crossover between Ising to Mean Field behavior. This is in principle 
the unique think that could destroy the picture we have presented in the previous sections. 
We have just noticed in sec. 3.1 how the positivity^ of the zero mode four legs vertex is 
related to the positivity of the zero mode mass near the critical point (see eq. 3.12). 

Using equations (6.16) it is easy to compute the critical eigenvalue Aq of P~^. If we define 

m' = ^ (6.18) 



''in this case the auxihary fields are two instead of five. However relation (6.15) is general once A is taken 
to run on the number of auxiliary fields present in the theory 

^This is related to the stability of the weakly theory in the sense that if the positivity condition fails 
than one needs to consider also higher than four legs vertex that have been discarded in chapter 2 simply 
using scaling arguments. 
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we find 

so that our stability condition becomes m' < (we have just noticed that ao3 > 0), and 
using eq. (6.7) we get the following condition 

v]!^ < -lOu (6.20) 

Figure 6.1 gives numerical evidence that the critical point we have considered satisfies pre- 
vious relation (6.20). However we expect that the same result holds on all the wing hues. 



Appendix A 

Properties of xQ eq. (5.11) 



We want to proof relations (5.15). Due to the symmetry Q{—x, T) — Q{x, T) we will consider 
here only x > 0. Notice that 



1 



where we have defined: 

wM = \/p^+^ (A.2) 



Then using the previous we have: 



dw^(p) 


X 


dx 






= 1 ( 


da;2 






3 


dx^ 





(A.3) 



3 d^.(p) / / d^;,.(p) y 
(p)2 dx \ V da; / 
3 dw;a;(p) d^w;^.(p) 



Wx{p) dx dx'^ 



(A.5) 



We have x < ^'^^(p) (A.2), so that both the first and second derivative or w are positive 
while the third derivative is always negative. Then we have: 

d^(a;,p) _ dV(p)^^^^^^J_/d«ap)y ^^31^-2^ ^A.6) 



dx dx2 2T 2T\ dx J 2T 

which is always positive for the previous considerations. For the second derivative we find: 

d'^frix, p) _ d^Wx(p) ^^^-^ Wxjp) _^ 3 dwa;(p) d^w^(p) ^^^^^-2 f^x(p) 



dx^ dx^ 2T 2T dx dx"^ 2T 

/ dw^(p) y 

\ dx J 



1 ^d!^\%i,h ^ cosh-3 
2T^\ dx J 2T 2T 
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3 dw^(p) d^w^(p) / w^(p) w^{p) wjp) 
'w^{p) dx dx^ V 2T 2T 2T 



2^2 V dx J 2T 2T 



< (A.7) 
having used the fact that: 



tanha; = - — — > 2 — (^-8) 



1 sinh 2x ^ x 
2cosh^a; ~ cosh'' a; 

Immediately (A. 6) and (A.7) give us the relations (5.15). Otherwise similar results show 
that Q{x,T) 



g{x,T) = I Tdp-l-tanh^ (A.9) 
2 J Wx[p) 2T 



is a decreasing function of x. Indeed we find 

dg{x,T) 1 1 du7^(p) / , tfx(p) w^{p) .-2'Wx{p) 



tanh - cosh-' (A. 10) 



dx 2 J w^{pY dx \ 2T 2T 2T 

and the claim follows using the (A. 8). 



Appendix B 

Relations among the Yukawa model, 
medium-range models, and field 
theory 



In this appendix we want to relate the weakly coupled c/?^ theory, medium-range models, and 
the Yukawa model for d = 2. This models have been considered in this work as prototype 
for large- A/" phase transition with infrared singularities. More interesting the solution pro- 
posed (i.e. to resum the divergences introducing iV-dependent scaling variables with scaling 
relations given by a weakly coupled (p^ theory) has given us the possibility to identify the 
crossover between Ising behavior to Mean Field behavior in all such models. In this appendix 
wc want to chicidatc how one can use numerical available data on medium range model to 
obtain numerical informations on Yukawa models (or vice- versa). 

For simplicity, wc only consider the case H = 0. corresponding to M = in the Yukawa 
model. The field-theory model has been discussed in sec. 2.1, where it was shown that the 
n-point zero momentum connected correlation function XFT,n shows a scaling behavior of 
the form 

u'^~^XFT,n = /FT,n(^FT) ^FT = [t - rc{u)]/u. (B.l) 

Next, we consider systems with medium-range interactions. Consider a square lattice, Ising 
spins ax at the sites of the lattice, and the Hamiltonian 

7i = -l^J(x-y)<7.ay. (B.2) 

We assume-*^ that J(x) = 1 for < Rm, ^(x) = for \x\ > Rm- The behavior of these 
models is very similar to that observed in the Yukawa model, Rm playing the role of N . For 
any finite i?^, the system belongs to the Ising universality class, while for R^ = oo all spins 
are coupled together and one obtains mean-field behavior. In Ref. [33] it was shown that 
this model shows a crossover that interpolates between mean-field and Ising behavior. If one 
defines an effective interaction range R by 

^.^ E,.(x-y)V(x y) _ 



^One can also consider a much more general class of medium-ranged models, see Ref. [32]. 
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then for R, R„ 



oo,t = {T- T^{R))/T^{R) ^ at fixed iuR = tR^ one lias 

R^ ^"XMR.n = fMR,n(tMR): 



(B.4) 



where XMR,n is the connected zero-momentum n-points correlation function of the fields a. 
In Ref. [32] it was shown that fMR,n{x) and fFT,n{x) are closely related. Indeed, we have 



fMR,n{^) — /^l,MRA*2,MR/FT,n('^MR2^), 



(B.5) 



where Hi^uR and Amr are model-dependent constants that reflect the arbitrariness in the 
deflnitions of the flelds, of the range R, and of the scaling variable t. The constants can be 
computed using the results of Ref. [32], sec. 4.2.^ Explicitly, we have for the two-point and 
four-point zero-momentum connected correlation functions: 



Xmr,2-R ^ 

XMR,4-R ^ 



^MR 



+ 



"-MR 



In 



AlTt 



MR 



+ 871^2 + 3 



+ 0(tMR), 



+ 0(fMR)- 



(B.6) 



MR 



In the field-theory model we have instead 

1 



UX2 



1 



In ^ + 87TD2 + 3 



+ 0(r'). 



Comparing we obtain 



/^l,MR 



A^MR — A 



MR 



1 

2' 



(B.7) 
(B.8) 

(B.9) 



In sec. 5.5 we have shown a similar relation for the Yukawa model. If Xt = —goiN[T — Tc{N)], 
we find for the zero-momentum correlation functions of the field 



(B.IO) 



and 



fY,n{^) = /f^i,yy"2,y/FT,n(Aya;). (B.ll) 

In order to compute these constants we compare the one-loop expansions of the two-point 
function in field theory and in the Yukawa model. In the Yukawa model we find 

2 



d'^x(x(0)x(x)) 



NT 



NT 
oPxt 



NT 



2N \a^xtj Jp^A (27r) 



d^p y^'^p, -p, 0, 0) 



^(p) 



+ 0{x;'). (B.12) 



^Note that the function (i) defined in Ref. [32] refers to correlations of the fields (j) and not of the original 
fields (fi. However, relation (4.12) of Ref. [32] shows that in the critical crossover limit J2xi'fio'Px) ~ J2xi'l^o(t>x)- 
The same holds for xa- The expression reported here are obtained from those reported in Ref. [32] by setting 

a2 — l,a4 — —2, N — 1, Co — Co — T, and ^mr — i + cq. 
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Using the explicit expression for Tc we obtain 



1 , q;V4(0,0) 

X2 \ 1 2 



In 



3a^V4(0, 0) 

For the four-point function we have instead 



+ 87TD2 + 3 



XI 



Comparing we obtain 



a^V^iO, 0) ' 

The constants reported in sec. 5.5.3 are easily derived: 



aW4{0, 0) 



K2,FT — ^J'1,Y^^2,Y^ 



LBB 



A 



MR 



n,LBB — 



fJ'l,Y(J^2,Y 
/^1,MR/^2,Mr' 



(B.13) 



(B.14) 

(B.15) 

(B.16) 
(B.17) 



where n — 1,2. Note that Qqi is negative, so that KpT and i^^LBB are positive as expected. 

The generalization to jV = 1 Heisenberg models is straightforward and can be obtained 
repeating the previous considerations of the Yukawa model with which shares the same large 
N expansion. 
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